ABSTRACT 



Choi, Ho-Meoyng. Light-Front Quark Model Analysis of Electroweak Decays of Pseu- 
doscalar and Vector Mesons. (Under the direction of Chueng-Ryong Ji.) 

We investigate the electroweak form factors and semileptonic decay rates of pseu- 
doscalar and vector mesons using the light-front constituent quark model. Exploring 
various detailed schemes of model building, we attempt to fill the gap between the 
model wave function and the QCD-motivated Hamiltonian which includes not only 
the Coulomb plus confining potential but also the hyperfine interaction to obtain 
the correct p-vr splitting. The variational principle to this Hamiltonian significantly 
. constrains our model. For the confining potential, we use both (1) harmonic oscil- 

lator potential and (2) linear potential to compare the numerical results for these 
two cases. Furthermore, our method of analytic continuation in the Drell-Yan-West 



a^ 

> 
o 

00 



P^l {q^=q^ + g^=0) frame to obtain the weak form factors avoids the difficulty associ 



ated with the contributions from the nonvalence quark-antiquark pair creation. This 



method seems to provide a promising step forward to handle the exclusive processes 
^ ■ in the timelike region. In the exactly solvable model of scalar field theory interacting 

with gauge field, we explicitly show that our method of analytic continuation auto- 
matically yields the effect of complicated nonvalence contributions. Our calculation 
reveals that the zero mode contribution is crucial to obtaining the correct values of 
the light-front current J~ in the g^=0 frame. We quantify the zero mode contribu- 
tions and observe that the zero mode effects are very large for the light meson form 
factors even though they may be reduced for the heavy meson cases. The mass spec- 
tra of ground state pseudoscalar and vector mesons and mixing angles of uj-(j) and rj-rj' 
are predicted and various physical observables such as decay constants, charge radii, 
radiative and semileptonic decay rates, etc., are calculated. Our numerical results are 
overall in a very good agreement with the available exprimental data. 
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2.1 Three different mixing scliemes for rj and rj' and the corresponding 


spin-averaged masses. 





2.2 Radiative decay widths for the V{P) P(V)7, Ai — > 717 and P 


27 for various model parameters P (in unit of GeV) with the r]-r/ 


mixing angle Osuia) = -10° 


-23° 


. The experiment is the result from 
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, unless otherwise noted. The unit of decay width is [keV], unless 



otherwise noted. 



3.1 Best fit quark masses and the model parameters p. PL J and PLH 



are the power law(PL) models for s=2 in Eg. (3.4) with and without 
Jacobi factor, respectively. q=u and d' 



3.2 Decay constants for -k ^ iiiJ 


and V e~^e , where the mixing angle 


=' -3.3Ss used. H'iPlR 


') uses the parameters of J[PLJ) to see 


the ett'ect of the presence-absence of the Jacobi factor. The experimen- 


tal data are taken from Ref. 
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j.3 Radiative decay widths V\P 77) for rj-rj' mixing angle "9^ 



SU{Z) 



-19°. The results are obtained from Eg. (3.16|) 



3.4 Radiative decay widths for V{P) P(\^)7 transitions. We use tj-tI 



and uj-cj) mixing angles as Osuiz) = —19° and by = —3.3°, respectively 



1.1 Optimized quark masses (m^, m^) (in unit of GeV) and the Gaussian 



parameters fj (in unit of GeV) for both harmonic oscillator and lin- 
ear potentials obtained from the variational principle. 'I'he values in 



parentheses are results from the smearing function in Eg. (D.6|) instead 



of the contact term. q=u and d. 



1.2 Fit of the ground state meson masses with the parameters given in 
'I'able 4.1. Underlined masses are input data. 'Ihe masses 01 (77,77'! 



and {u, (j)) were used to determine the mixing angles of rj-rj' and cJ^ 
respectively. The values in parentheses are results from the smearing 



function in Eg. (D.6| ) instead of the contact term. 



Decay constants (in unit of MeV) and charge radii (in unit of fm"') foi 
various pseudoscalar and vector mesons. For comparison^ we use Oy = 
3.3" ± 1" for both potential cases. 'I'he experimental data are taken] 



|49|| , unless otherwise noted. Note that fx* was obtained by 



from Ref. 



4.4 Radiative decay widths for the V{P) 


P{V)'y transitions. The mix- 


ing angles, esjK?,^ = 


— 19° for rj-T]' and 


Sv 


= 3.3° ± 1° for LJ-(t), are 


used for both potential 


models, respectively. The experimental data 


are taken from Ref. 


49 





4.5 Radiative decay widths r(P — > 


77) obtained by using the axial anomaly 


plus PCAC relations (see Eq. 


|3.16)). 0su(3) = —19° for 77-77' mixing is 


used for both potential cases. 


The experimental data are taken from 


Ref. 


49 





3.1 



f'brm factors of f±{0) obtained for different zer o- momentum transfer 
limit, q'*~=0 and q^=0. The notations of a± 



used in] 



O-r. 



table are defined as a± = a±{0), ap=l + a-{0), and am=i ~ a-(0) 



respectively. J,- 



(a+) implies yj_, /f (a+) 



372 Zero-mode(Z.M.) and nonvalencefN.V.) contributions to the exact form 



factors of f±iO) for the scmilcptonic decays of K{B) n and B ^ D 



m 



1 + 1) dimensions. We distinguished the zero mode contribution at 
q^=Q from the usual nonvalence one at (/~=0 



7. 1 Optimized quark masses rrip 


GeV 


for both HO and linear potentials ob- 


taincd from the variational principle. We also include the results from 


the smearing function(yi'') 


instead of the Breit-i'ermi contact term. 


q=u and d. 





7.2 Optimized Gaussian parameters (3 


OeV 


for both HO and linear po- 


tentials obtained from the variational p 


rinciple. We also include the 


results from the smcarmg tunction(yi* ) mstcad of the Breit-i'ermi con- 


tact term. 





f''it of the ground state meson masses with the parameters given in 
I'ables 7.1| and [7.2| . Underline masses are input data. The masses of 
[uj-cf)) and (77-7/') were used to determine the mixing angles of uj-(j) and 
7/-7/', respectively 



and charge radii |fm^| for various heavy pseu- 



Decay constants 

dosca iar and vector mesons. INotc that the jp and }y are obtained by 



(0 1 ^27^75^1 1 -P) = ifpP^ and {Q\q2l^qi\V) = ^v.fy^_^; respectively, to 
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(129.8) for the HO (1 inear) parameters 
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m units of 10" s ^) is 
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I.l The LFQM description of a electroweak meson form factor: (a) the 
usual LF valence diagram and (b) the nonvalence( pair-creation) dia- 



gram. The vertical dashed hne in (b) indicates the energy-denominator 
for the nonvalence contributions. While the white blob represents the 



usual Lf*' valence wave function, the modeling of black blob has not yet 



been made. 



I.l 


a) A{Q^) testing the angular condition is shown as a function of Q^. 


'i'he solid lines are our results with jj = 0.32 and 0.36 GeV. i'or com- 




3arison, we include various choices of the wave function w'' introduced 




by Godfrey and Isgur(Gl) in Ref. 17| . 'i'he dashed, dotted and dashed- 



dotted lines correspond to (spin-independent part), wgi (OGE -|- 
linear confining term) and Wconf (linear confinging term), respectively, 
(b) The form factors of p meson with parameter (3 = 0.36 GeV. The 
sohd, dotted, and dashed lines correspond to Fc{Q^),Fm{Q'^) and 
Fq{Q'^), respectively. The dashed-dotted line is the result of non- 
relativistic limit of FclQ"^) by turning off the Melosh transformation. 20 



i.l (c) 'i'he form factors, G'l = Fi, G'2 = F2 - Fi and G'3 = Fs/2, oi 



p meson are compared with the results of QGD sum rules 50| |. The 
dotted, dashed and dashed-dotted lines correspond to G'l, G'2, and G3 
of Ref. O 21 



2.2 (a) The form factors of Ai meson with the parameter = 0.36 GeV" 



he solid, dotted, and dashed lines correspond to j^c*, -^'a/, and F^, 



respectively. The quantity A{Q^ is shown as dashed-dotted line 



(b) The form factors, d = Fi, G2 = F2 - Fi and G3 = F3/2, ol 
Ai meson are compared with the result of QCD sum rules 50|| . The 



solid, dotted, and dashed lines are our predictions of Gi, G2/2, and 
G3, respectively. The QCD sum-rule results are Gi((5^)(long dashed 
line) and G2(Q^)/2 (dashed-dotted line) 23 
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2.3 (a) The 


transition form factor of 7r^7* multiplied by Q is shown 


as a function of Q"^. The sohd hnes are our results with (3 = 0.32 and 


0.36 GeV. The dotted hne, dashed, and dashed-dotted hnes correspond 


to the results obtained using wqi in Ret. 


17 


, the BS approach of 

lG,^ = i/(i+gV^p)- 


Ref. |60| 


and the predictions of VDM model with 



(b) The normalized transition form factor of — ir^y*. Our result with 
/5 = 0.36 GeV (solid line) is compared with VDM (dashed line) with 
G^n = 1/{1 + Q^/Mp) and the pole fit of the experimental data (dotted 
line) il 



26 



2.4 The normalized form factors of K*"^ — » K ^7* (solid line) and K*^ 



-R"'^7*( dotted line) transitions with the parameter P= 0.36 GeV. The 



dotted and dashed lines correspond to charged and neutral vector kaon 
transition form factors by Ref. |62||. 'i'he dashed-dotted line corre- 



sponds to the VUiVl prediction. ~ 28 



2.5 'i'he form factors of Ai n'^'-f* transition with the parameter fj = 0.36 


GeV. GiiQ"^) and G2{Q'^) /W ot our results arc represented by the solid 


and dashed lines, 


respectively, l^'or comparison, we show also the QGU 


sum-rule results 


50 


: Gi(Q') (dotted-line) and G2(Q')/10 (dotted- 



dashed line) 31 



2.6 The decay rate tor tt"^ —>■ 7*7 transition with the parameter fj = 0.35 



GeV. Data are taken from Rets. |51|, |52||. 34 



2.7 



'I'he decay rate for rj — » 7*7 transition with the parameter fJ = 0.36 
GeV. 'I'he solid and dotted lines correspond to the t)su(3) = — 10" and 



—23° mixing schemes, respectively. Data are taken from Rcfs. |51| , . 34 



2.8 The decay rate for t]' — » 7*7 transition with the parameter (j = 0.36 
GeV. 'I'he solid and dotted lines correspond to the tlsu(3) = —10" and 



-23" mixing schemes, respectively. Data are taken from Refs. [51|, . 35 



3.1 (a) The normalized quark distribution amplitude of the pion in HO 



model. The solid, dotted, dot-dashed lines represent the results of the 
models H, H', and J, respectively. The normalization is fixed by the 



zeroth moment(the area underneath each curve) {^'^}= JI^/ d^^'^(pTr{^)=l 



with ^=Xi — X2- (b) 'I'he normalized quark distribution amplitude ol 



the pion in PL model, i'he solid, dotted, dot-dashed lines represent 



the results of the models PLH, PLH' , and PL J, respectively. The 
long-dashed line represents the Schlumpf 's result [19| 
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3.2 (a) 'I'he square of the pion charge form factor in 


HO model for low 


values of compared with the experimental data 


70 . The same line 



code as in Fig. 3.1(a) is used, (b) The square of the pion charge form 
factor in PL model. The same line code as in Fig. 3.1(b) is used. ... 47 



4.1 (a) The central potential Vo(r) versus r. Our Coulomb plus HO (solid 



line) and linear (dotted line) potentials are compared with the quasi 
relativistic potential of 1SGW2 model with k = 0.3 (long-dashed line 
and K = 0.6 (dot-dashed line) and the relativized potential of Gl model 
(short-dashed line 



Our force f 



(b) The central force fo[r) versus r. 

the linear potential is the same as that oi \m with K = 0.3 
(dotted lines). The forces of GI [^J and ISGW2 with k = 0.6 are 
the same as each other (dashed lines). Our force for the HO potential 
(solid line) is quite comparable with the other four forces up to the 
range of r < 2 fm 



oi 



12 (a) The parameters nif^ 




/^aa; f^cis, and fJss satisfying the variational 


principle given by Eq. 


(D.2 


). The solid, dotted, and dot-dashed lines 


are fixed by the sets of 


rriu, 


l3uu), ims,(3us), and ims,(3ss), respectively. 


(b) The parameters rn,. 




Pan- 0as^ and (3ss satisfying the variational 


principle given by Eq. 


(D.3 


. The same line codes are used as in (a). 





4.3 The charge form factor for the pion compared with data taken from 


Ref. 




. 'I'he solid and dotted lines correspond to the results of HO 



and linear potential cases, respectively. 



4.4 (a) Theoretical predictions of charged form factors using the pa" 



rameters of both HO (solid) and linear (dotted) potentials compared 
with a simple two-pole VMD model |8"7| (dot-dashed) 



pVDM 
^K+(KO) 



eu(d)ml/{ml + Q"^) + esm^^/iw?^ + Q"^). (b) Theoretical predictions 
of neutral form factors. The same line codes are used as in (a). . . 





4.5 The decay rate for the vr — > 7*7 transition obtained from the one-loop 


diagram. Data are taken from Refs. 


51 




52 





4.6 'I'he decay rate for the 77 7*7 transition obtained from the one-loop 



diagram. Data are taken from Rets. |51| , |52| , ^3| 





4.7 The decay rate for the 77' 7*7 transition obtained from the one- loop 


diagram. Data are taken from Refs. 


51 




52 
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5.1 The electromagnetic decays of a photon into two-body bound sates, 
i.e., 7* OT(or QQ)— » M.{qQ)M.(qQ), in scalar theory: Covariant 
representation (a), and the LF time ordered contributions to the decay 
amplitude (b) and (c).| 

5.2 Covariant triangle diagram (a) is represented as the sum of LF triangle 
diagram (b) and the LF pair-creation diagram (c). 



5.3 (a) 'I'he electromagnetic form fator of the pion in (3 + 1) dimensional 
scalar field theory for —2 < < 3 GeV^. The total, real, and imagi 



nary parts of ^^(g^) ^ are represented by thick solid, solid, and dotted 
lines, respectively, (b) The electromagnetic form fator of the pion in 



(3 + 1) dimensional scalar field theory for the entire timclikc region 
compared to the dispersion relations (data of black dots) given by 



Eg. (5.25). The same line code as in (a) is used 



5.4 The electromagnetic form fator of the kaon in (3 + 1) dimensional 



scalar field theory for —2 <■ q'^ < 5 GeV^. The same line code as in 



Fig. 5.3(a) is used, (b) The electromagnetic form fator of the kaon 
in (3 + 1) dimensional scalar field theory for the entire timelike region 



compared to the dispersion relations (data of black dots) given by 
Eg. (5.25). The same line code as in Fig. 5.3(a) is used 



aTS 'i'he electromagnetic form fator of the D meson in (3 + 1) dimensional 



scalar field theory for —10 < (/^ < 30 GeV^. The same line code as in 



Fig. 5.3(a) is used 



3.1 (a) Normalized form factor of f+jq'^) for -R" ^ tt in (1 + 1) dimension" 



The solid line is the result from the valence plus nonvalence contribu- 
tions. 'Fhe dotted line is the result from the valence contribution, (b) 



Normalized form factor of f-{q'^) for — » vr in (1 + 1) dimension. The 



same line code as in (a) is used 



3.2 (a) Normalized form factor of ,f+{q^) for i? ^ tt in (1 + 1) dimension" 



The same line code as in Fig. 6.1(a) is used, (b) Normalized form 



factor of i-{q^) for ^ tt in (1 + 1) dimension. The same line code 
as in Fig. 6.1(a) is used 



3.3 (a) Normalized form factor of f+jq'^) for B ^ D m. [1 + 1) dimension" 



The same line code as in Fig. 6.1(a) is used, (b) Normalized form 



factor of f-{q'^) for B —>■ JJ in [1 + 1) dimension. 'I'he same line code 



as in Fig. 6.1(a) is used 



7.1 The parameters m^, m^, mi,, Pgg, (3^^, etc. satisfying variational prin- 
ciple with (a) HU and (b) Linear potential models. 'I'he qq and qc etc. 
represents the sets of {mg^mg, (3aa) and {mg^mc, Poc) etc., respectively. 

7.2 The EM form factor of D meson for the spacelike = —q"' regiom 
The thick solid (light solid) and thick dotted (light dotted) lines are the 
results from the ^'"'"=0 jq^^O) frame for the HO and linear potential 
parameters, respectively. 7.3: The EM form factor of B meson for the 
spacelike = —q^ region, same line code as in f^'ig. 7.2 is used. I 



The EM form factor of pion for low Q'^ = —q^ compared with data |70 . 
7.5: The EM form factor of kaon compared with data [5^. The same 
line code as in Fig. 7.4 is used 





7.6 The form factors f+iq'^) for the K ^ n transition in timelike momen-| 




tum transfer q^ > 0. 'i'he solid and dotted lines are 


the results from the| 


(/^=0 and (/^T^O frames for the flU and linear parameters, respectively. 


The differences of the results between the two frames are the measure] 


of the nonvalence contributions from q'^^Q frame. 





rr 



'I'he decay rates dV/dq of i(e3(solid line) and i(^3(dotted line) for the 
HO parameters in q~*~=0 frame. The dot-dashed line is the contribution 



from the term proportional to /q in Eg. (7.8| ) for K^^ decay. 'I'he results 
for the linear parameters are not much different from those for the HO 
parameters 



7.8 The form factors /+( 


g^) for the K ^ n transition in spacelike momen- 


tum transfer — (/^ 


< 0. The solid and dotted lines are the results from 


the HO and linear parameters, respectively. The dot-dashed line is the 


result from Kef. 
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7.9 (a) The form factor f+{q^) for the D ^ n transition compared with a 


lattice calculation 


109 


. The same line code as in Fig. 7.2 is used for 



our results, (b) The form factor fo{q^) 
same line code is used as in fa). ... 



for the D ^ Tx transition. The 



7.10 (a) The form factor f+{q^) for the D ^ K transition compared with 


the experimental data (•) 


49^ 


and the lattice calculation 1091 


. The 



same line code as in Fig. 7.2 is used for our results, (b) The form factor 
foiq"^) for the D ^ K transition. The same line code is used as in (a). 



Ml The form factors V(g^) (solid line), ^ijq'^) (dotted line) and A2{q'^) 
( long-dashed line), where the thick (thin) lines represent the HO (lin- 



(long-dashed line), where the thick (thin) lines represent the HO 
ear) parameters, for D ^ p transition compared with the lattice cal 
culation |110|| oi Ai[q') (data). 7.12: The form factors V{q^), A^jq' , 
and A2{q'^) for D — >■ K* transition compared with the lattice calcula- 
tions ||105|| (square) and ||110|| (circle) of Ai{q^). The same line code as 



in Fig. 7. 



1 is used. 



7.13 'i'he form factor f<(q 



^ ) for the Vg V transition with the rj-rj' mixin 
llic same Imc code as m iig. 7.2 is used tor ou 



angle Usi 



results. 7.14: 'i'he form factor f+[q^} for the Vg — > r]' transition wit 



the r]-r] mixing angle 9, 



su(3y- 



T9°. The same line code as in Fig. 7. 



is used for our results 



7.15 The form factors Vjq'^), v4i(g^) and ^2(g^) fo^^ Ds (pi^s) transitioii 
The same line code as in Fig. 7.11 is used, 



7.16 (a) The form factor f+jg'^) for the — » tt transition compared witE 
the lattice calculation 108|| (data). The same hne code as in Fig. 7.2 
is used for our results, (b) The form factor /o(q'^) for the B 
transition. The same line code is used as in (a) 



IT 
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7.17 (a) 'I'he form factor V{q'^) for the B"^ p transition. 'I'he thick (thin 
solid line corresponds to our HO (linear) parameters, f'br comparison, 



we include the results of J-'A,Ay {q^) (dashed hne) given by Eg. (7.27 ) and 
lattice QCD from UKQCD ||103|| (•), APE ||105|| (square) and ELC |TTT 



(o). (b) The form factor Ai(g^) for the B^ ^ p transition. The same 
line code as in (a) is used 134 



7.17 (c) The form factor A2[q'^) for the B" ^ p transition. The same line 



coae as m a is usee 



7.18 (a) The form factor f+{q'^) for B—>-D transition. The same line code 



as m Fig. 7.2 is used, b The form factor foiq 
ihe same line coc 



for B^D transition. 
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le as in"FTg. 



IS use 



for B^^D^ transition. The same line code 



7.19 (a) The form factor f+{q"^, 

as in Fig. 7.2 is used, (b) The form factor foiq'') for Bg 



>D^ transition 



The same line code as in Fig. 7.2 is used7 



7.20 The IW function ^^^{vi ■ V2) for the HO (solid) with the common 



gaussian (3 = /3b=0.496 and linear (dotted line) parameters with the 
common (j = /i/^=0.5266, respectively. We compare our results with 
the experimental data of ARGUS |140|| (square) and CLEO |P-41|| (circle) . 141 





7.21 The dr/dq'^ distribut 


ion for B D*iv decays using the linear param- 


eters. The thick solid, dotted, and long-dashed lines are the results 
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data are taken from Ref. 
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7.22 The form factors V(g^), Ai((7^) and A2((7^) for B ^ D* transition. Th 
same line code as in Fig. 7.11 is used. 7.23: The form factors Viq"^) 



A'liq'^) and ^2(0'"') for B D* transition. The same line code as in| 
I'lg. 7.11 is used.] 144 

7.24 The two form-factor ratios, Rijw) (solid line) and R2{w) (dotted line"] 
for the HO parameters compared with the experimental data" 149|| of 



for the HO parameters compared with the experimental data "14911 of 
-Ri(l) (circle) and -R2(l) (square). 7.25: Our predictions of the ratio 
of ^^^{w)/E,^^ {w) for both HO (solid line) and linear (dotted line) 
parameters compared with the experimental data [59| 145 



1 



Chapter 1 
Introduction 

Quantum Chromodynamics (QCD) which describes the strong interactions of 
quarks and glouns, has been a very active field of research in nuclear and parti- 
cle physics. To calculate the hadronic properties from QCD, however, is such a 
formidable task that progress has been slow and many problems are still not solved 
today. The equations to be solved are highly non-linear and one is dealing with for- 
mulas for infinitely many relativistic particles, i.e., a quantum field theory requiring 
careful renormalization. Thus QCD presents one of the most challenging problems of 
theoretical physics. 

Presently, however, this situation is changing. The continuous slow progress in com- 
putational and analytic techniques has now reached a point where real QCD calcu- 
lations have become feasible for a rapidly increasing number of problems in hadron 
physics. Together with the much improved experimental situation this should trig- 
ger an enormous boost for hadron physics, e.g., various static and non-static hadron 
properties such as form factors, structure functions and cross sections etc., in the 
years ahead. These results are thought to be very important as they will help us in 
understanding the role of QCD or the quark substructure and their interactions inside 
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the nuclear matter. For best comparison with the experiment, rehab le and versatile 
theoretical calculations have to be at hand. To avoid the extreme complexity of real 
QCD, many phenomenological "QCD-inspired" models have been developed, often 
with astonishing success. 

The main purpose of this thesis is to develope "QCD-inspired" model describing the 
mass spectrum and hadron properties of low-lying hadrons at small momentum trans- 
fer region(nonperturbative QCD). Even though more fundamental nonperturbative 
QCD methods such as QCD sum-rule techniques [Q, ^, Q and lattice QCD calcula- 
tions 0, ^ are available, there is still growing interest in using simple relativistic 
quark models |, |, [ig, |l7l |ll l^, ||, |2| to describe 



hadron properties. Our relativistic quark model [|1^, [T^, |2T], |22[ is based upon the 
Fock state decomposition of hadronic state which arises naturally in the "light-front 
quantization" of QCDQ. 



1.1 Light-Front Degrees of Freedom 

The foundations of light-front (LF) quantization date back to Dirac ||2^ who showed 
that there are remarkable advantages in quantizing relativistic field theories at a 
particular value of LF time t = t + z/ c rather than at a particular ordinary time t. 
In this formalism, a hadron is characterized by a set of Fock state wave functions, the 
probability amplitudes for finding different combinations of bare quarks and gluons 
in the hadron at a given LF time r. For example, a meson can be described by 

1^) = 5Z i^qq/M + '4^qgg\qq9)'^qqg/M + " " " • (1-1) 

qq qqg 

^The interested readers are referred to extensive review by Brodsky, Pauli and Pinsky [p3| 
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These wave functions provide the essential hnk between hadronic phenomena at short 
distances(perturbative) and at long distances(non-perturbative) [^]. The quantum 
field theory at equal-r is quite different from the one at equal-t. The distinguished 
features in the LF approach are the dynamical property of rotation operators and 



the simplicity of the vacuum except the zero modes |2^, |2^. Let us begin by 

illustrating these two distinguished features. 

The boost operators of the ten Poincare generators contain interactions (i.e. dynam- 
ical operators) changing particle numbers in equal t. Thus, solving the relativistic 
scattering and bound state problems in the cannonical quantization at the equal t, 
depends on the reference frame due to the boost non-invariance when the Fock-space 
is truncated for practical calculations. The problem of boost operators changing par- 
ticle numbers at equal t can be cured by the LF quantization since the quantization 
surface r = is invariant under both longitudinal and transverse boosts (i.e. kine- 
matical operators) defined at equal r. In return, however, the rotational invariance 
is violated at equal r when the Fock space is truncated for practical calculations [^ . 
The transverse angular momentum whose direction is perpendicular to the direction 
of the quantization axis z in equal r involves interactions (i.e. dynamical operators) 
and thus it is not easy to specify the total angular momentum of a particular hadronic 
state. Also r is not invariant under parity |31|. The number of problems seems to 
conserve even though the boost problem at equal t is now replaced by the rotation 
problem at equal r. However, we point out that the rotation problem is much easier 
to deal with compared to the boost problem because the rotation is compact, i.e., 
closed and periodic [Bl]. Moreover, the problem of assigning quantum numbers J^^ 



to hadrons can be circumvented by the Melosh transformation [^. This ought to be 
regarded as an advantage of the LF quantization over the ordinary equal t quanti- 
zation. The restoration of the rotational symmetry in the LF quantization has also 
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been discussed [031. 



The characteristics of vacuum at equal r has also a dramatic difference compared 
to the vacuum properties at equal t. Suppose that a particle has mass m and four- 
momentum k = {k^, k^, k"^, k^), then the relativistic energy-momentum relation of the 
particle at equal-r is given by 

r ^ !i±^, (1.2) 

where the LF energy conjugate to r is given by k~ = k^ — k^ and the LF momenta 
fc"*" = k^+k^ and k_L = {k^, k"^) are orthogonal to k~ and form the LF three-momentum 
k = (A;^,k_|_). The rational relation given by Eq. ( |1.2| ) is drastically different from 
the irrational energy-momentum relation at equal-t given by 



= Vk2 + m2, (1.3) 

where the energy k^ is conjugate to t and the three-momentum vector k is given by 
k = {k^, k"^, k^). The important point here is that the signs of k^ and k~ are correlated 
and thus the momentum fc+ is always positive because only the positive energy k" 
makes the system evolve to the future direction (i.e. positive r), while at equal-t the 
signs of k^ and k are not correlated and thus the momentum k^ corresponding to k~^ 
of equal-r can be either positive or negative. This provides a remarkable feature to 
the LF vacuum; i.e., the trivial vacuum of the free LF theory is an eigenstate of the 
full Hamiltonian, viz., the true vacuum |^4|. This can be proved by showing that the 
full LF Hamiltonian annihilates the trivial perturbative vacuum [^. For example, 
in QED, the application of the interaction Hip = J d^xjip'-y^ipA^ to the perturbative 
vacuum |0) results in a sum of terms {ki)a\k2)d) {k:^)\()) . While the conservation of 
the LF momentum requires Yl^=i = 0, the massive fermions with finite k~ cannot 
have kf = due to Eq. (|1.2|) . Thus, Hip annihilates the trivial vacuum |0) and so 
does the full Hamiltonian H^p = Hip + H^fj since |0) is annihilated by the free 
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Hamiltonian H^p by definition. This feature is drastically different from the equal-t 
quantization where the state i^|0) is a highly complex composite of pair fiuctuations. 
This greatly complicates the interpretation of the hadronic wave functions at equal-t 
quantization. 

This apparent simplicity of the LF vacuum may yield a problem to understand the 
novel phenomena such as the spontaneous symmetry breaking, Higgs mechanism, 
chiral symmetry breaking, axial anomaly, ^^-vacuua, etc., because these were known as 
the direct consequencies of the nontrivial vacuum structures of various field theories. 
Thus, the question of how one can realize these nontrivial vacuum phenomena from 



the trivial LF vacuum arises ||3^. In fact, some efforts have been made to show that 
nontrivial vacuum phenomena can still be realized in the LF quantization approach 
if one takes into account the nontrivial zero-mode (A;^ = 0) contributions 37]. As 



an example, it was shown that the axial anomaly in the Schwinger model can be 



obtained in the LF quantization approach by carefully analyzing the contributions 
from zero modes. 



1.2 Construction of Light-Front Quark Model 

The distinguished features in the LF approach can now be summarized as (1) the 
dynamical property of the rotation operators and (2) the simplicity of the vacuum 
except the zero modes. Thus, in the LF quantization approach, one can take advan- 
tage of the rational energy-momentum dispersion relation and build a clean Fock state 
expansion of hadronic wave functions based on a simple vacuum by decoupling the 
complicated nontrivial zero modes. The decoupling of zero modes can be achieved 
in the light-front quark model (LFQM) since the constituent quark and antiquark 
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(a) 



(b) 



Figure 1.1: The LFQM description of a electroweak meson form factor: (a) the usual 
LP valence diagram and (b) the nonvalence(pair-creation) diagram. The vertical 
dashed line in (b) indicates the energy-denominator for the nonvalence contributions. 
While the white blob represents the usual LF valence wave function, the modeling of 
black blob has not yet been made. 



acquire appreciable constituent masses. Furthermore, recent lattice QCD results 38 



indicated that the mass difference between v( and pseudoscalar octet mesons due to 
the complicated nontrivial vacuum effect increases (or decreases) as the quark mass 
rUq decreases (or increases); i.e., the effect of the topological charge contribution 
should be small as increases. This supports us in building the constituent quark 
model in the LF quantization approach because the complicated nontrivial vacuum 
effect in QCD can be traded off by the rather large constituent quark masses. 



One can also provide a well-established formulation of various form factor calculations 
in the LF quantization method using the well-known Drell-Yan-West (g"*" = + = 
0) frame |3^, ^ , which provided an effective formulation for the calculation of various 
form factors in the space-like momentum transfer region = —Q^ < 0. In g+ = 
frame, only parton-number-conserving Fock state (valence) contribution is needed 
when the "good" components of the current, J+ and Jj_ = {J^, Jy), are used [^]. For 
example, only the valence diagram shown in Fig. 1.1(a) is used in the LFQM analysis 
of spacelike meson form factors. 
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The key approximation of the LFQM is the mock-hadron approximation P2| to trun- 



cate the expansion by retaining only the lowest Fock state and treat the lowest Fock 
state as a free state as far as the spin-orbit part is concerned while the radial part 
is given by the ground state of the harmonic oscillator wave function. Then, the 
assignment of the quantum numbers such as angular momentum, parity and charge 
conjugation can be given to the LF wave functions by the Melosh transformation | p2|| . 
For example, the meson bound state \M) is represented by 

\M) = ^%\QQ), (1.4) 

where Q and Q are the effective dressed quark and antiquark. The model wave 
function in momentum space is given by0 



^QQ = = ^^<l>ixM)niXixM), (1.5) 

where 0(x, k_|_) is the radial wave function and dkn/dx is a Jacobi factor, and 7^^^\_(a;, 
kx) is the spin-orbit wave function obtained by the interaction-independent Melosh 
transformation from the ordinary equal-time static spin-orbit wave function assigned 
by the quantum numbers J^'-' . When the longitudinal component fc„ is defined by 
kn = {x — 1/2) Mq + (m? — m^)/2Mo, the Jacobian of the variable transformation 
{x, k±} ^ k = {kn, kj_) is given by 



dkn Mo ( r(m2-m?^i2 



Mi 



dx 4x(l — x) 
The explicit form of the spin-orbit wave function is given by 



;i.6) 



,i,,i 



X {\,\nl{l-x,-k^,m,)\\'^)q^\'^\JJs), (1.7) 

■^The wave function in Eq.( |l.5| ) is represented by the Lorentz- invariant variables Xi ~ p^/P^, 
kit = p±^ - x,P± and A,, where P'^ = (P+,p-,P^) ^ (po + P^ (Mg^ + Pl)/P+,Pj_) is the 
momentum of the meson M, and and Xi are the momentum and the helicity of constituent 
quarks, respectively. 
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where 



7?.M(a;,k| , m) 



m + xMq — icr ■ (fi X k) 



1.8) 



(m + xMo)2 + ki 

is the Melosh transformation operator with n= (0,0,1) being a unit vector in the z 
direction. 

Ahhough the spin-orbit wave function is in principle uniquely determined by the 
Melosh transformation given by Eq. (|1.8| ), a couple of different schemes for handling 
the meson mass Mq in Eq. ( |1.8| ) have appeared in the literature [§, |^, [1^, [TT], 

While in the invariant meson mass (IM) 



g, jig 0, |T5|, |I7|, y, ig, 

10| , m , |12| , [l^, |14|, |15|, |17[ |18|, |l^ , the meson mass square Mq is given by 



scheme P, 



k^ + ^ k^ + m? 



1.9) 



X 1 — X 

in the spin-averaged meson mass (SM) scheme [^U|, |2T|, Mq was taken as the 
average of physical masses with appropriate weighting factors from the spin degrees 
of freedom. 

This thesis is divided mainly into two parts, i.e., the analyses of the spacelike (g^ < 0) 
and the timelike (g^ > 0) form factors. In Chapter 2 of this thesis [22|, we first study 
the radiative decays of pseudoscalar (vr, K, r], t]'), vector (p, K*, u, 0) and axial vector 
{Ai) mesons using the SM scheme, which was originally argued by Dziembowsky and 
Mankiewicz |^0| and later developed by Ji and Cotanch [^]. The purpose of this work 



is to investigate more observables with the same model as Refs. |]2y, ^ . Our overall 
predictions of pseudoscalar, vector and axial-vector meson radiative decay processes 
are in remarkably good agreement with the experimental data. 



In Chapter 3 |]T2|, we discuss the relations among the LP quark models [0, |TT], |18|, |T9| 
which are based on the IM scheme. We noticed that there are a couple of differences 
within the same IM sheme, i.e., (a) the presence-absence of the Jacobi factor dkn/dx 
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in choosing harmonic oscillator (HO) wave function and (b) the difference in the 
choice of radial wave function, e.g., HO versus power- law wave functions. We found 
that the difference from the Jacobi factor is substantially reduced in the numerical 
predictions for the physical observables once the best fit parameters are chosen. The 
difference in the choice of radial wave function is, however, still appreciable even 
though the best fit parameters are used. Comparing the two different meson mass 
schemes, i.e., SM and IM schemes, we also found that once the best fit parameters 
were used both schemes provided the predictions that were not only pretty similar 
with each other but also remarkably good compared to the available experimental 
data for form factors, decay constants, charge radii, etc., of various light pseudoscalar 
and vector mesons as well as their radiative decay widths. Through the analyses in 
Chapters 2 and 3, we model the radial wave function rather than a potential. 



In Chapter 4 ||13[, we develop LFQM attempting to fill the gap between the model 
wave function and the QCD-motivated potential, which includes not only the Coulomb 
plus confining potential but also the hyperfine interaction, to obtain the correct p- 
TT splitting. In our LFQM we analyzed both the mass spectra and the wave 
functions of the light pseudoscalar and vector mesons and predicted mixing angles of 
u!-(f) and rj-i]' and various physical observables such as decay constants, charge radii, 
radiative decay rates, etc. 

The timelike (g^ > 0) form factor analysis in the LFQM is more subtle than the 
analysis of the spacelike form factors since the = frame is defined only in the 
spacelike region (g^ < 0). While the ^ frame can be used in principle to com- 
pute the timelike form factors, it is inevitable (if q'^ ^ 0) to encounter the nonvalence 
diagram arising from the quark-antiquark pair creation (so called "Z-graph"). For 
example, the nonvalence diagram in the case of semileptonic meson decays is shown 
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in Fig. 1.1(b). The main source of the difficulty, however, in calculating the nonva- 
lence diagram is the lack of information on the black blob which should contrast with 
the white blob representing the usual LF valence wave function. In fact, we recently 



found 1^ in the analysis of kaon electorweak (^^3) decays that the omission of nonva- 
lence contribution leads to a large deviation from the full results. The timelike form 
factors associated with the hadron pair productions in e^e~ annihilations also involve 
the nonvalence contributions. Therefore, it would be very useful to avoid encounter- 
ing the nonvalence diagram and still be able to generate the results of timelike form 
factors. This can be done by the analytic continuation from the spacelike form factor 
calculated in the Drell-Yan-West {q^ = 0) frame to the timelike region. 



Before applying our LFQM ||T3|| to the electroweak form factors and semileptonic 
decay rates of pseudoscalar and vector mesons, we start from the heuristic model 
calculations to compare the results of the timelike form factors obtained from both 
g"*" 7^ (i.e. valence + nonvalence contributions) and q'^ = (i.e. valence contribution 
only) and to show the analytic continuation method is correct. The choice of the 
components of the LF current is especially important when one needs to calculate 
timelike processes in g'^ = frame such as the semileptonic decays of pseudoscalar 
mesons. We also come to this point later in this thesis. 



In Chapter 5 of this thesis |^ , we investigate the form factors of qQ bound states both 
in spacelike and timelike region using an exactly solvable model of (3 + 1) dimensional 
scalar field theory interacting with gauge fields. Based on the LF quantization, the 
Drell-Yan-West {q~^=0) frame as well as the purely longitudinal momentum (g"*" 7^ 
and q_L=0) frame were used for the calculations of the M — > 7* + M form factors. 
We then analytically continue the form factors in the spacelike region to the timelike 
region and compare those with the direct results of the timelike 7* — M + M form 
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factors. Analytically continued results coincide exactly with the direct results and 
it verifies that the method of analytic continuation is capable of yielding the effect 
of complicate nonvalence contributions. The meson peaks analogous to the vector 
meson dominance (VMD) phenomena are also generated at the usual VMD positions. 



In Chapter 6 we discuss the zero mode(g^ = mode of a continuum theory) 
contribution, which is crucial to obtain the correct values of the LF current J~ in the 
Drell-Yan-West(g"'' = 0) frame. In the exactly solvable model of (l+l)-dimensional 
scalar field theory interacting with gauge fields, we quantify the zero mode contribu- 
tion and observe that the zero mode effects are very large for the light meson form 
factors even though they are substantially reduced for the heavy meson cases. 

Finally, we present in Chapter 7 [O, H4| the analysis of exclusive O^^O^(l^) semilep- 



tonic meson decays using our LFQM [13|. Our method of analytic continuation to 



obtain the weak form factors avoids the difficulty associated with the contribution 
from the nonvalence quark-antiquark pair creation. Our numerical results for the 
decay rates are in a good agreement with the available experimental data and the 
lattice QCD results. In addition, our model predicts the two unmeasured mass spec- 
tra of ^So{hh) and '^Si{hs) systems as M;,5= 9295 (9657) MeV and Mbj= 5471 (5424) 
MeV, for the HO (liner) confining potential, respectively. Conclusions and discus- 
sions of this thesis follow in Chapter 8. In Appendix A, the spin-orbit wave functions 
7l'l'^\_{x, kx) of pseudoscalar and vector mesons are presented for the IM scheme. In 
Appendix B, the details of how to obtain the spin-averaged masses of rj, rj', u and 
are presented. In Appendix C, we present the derivation of the formulas used for the 
electromagnetic decay widths. In Appendices D and E, we show how to fix the model 
parameters and determine the mixing angles of uj-(f) and rj-r]', respectively, from our 
QCD-motivated quark potential model discussed in Chapter 4. The derivations of 
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the matrix element of the form factors for — > semileptonic decays in both ^+=0 
and ^''''7^0 frames follow in Appendix F. 
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Chapter 2 

LFQM Predictions for Radiative 
Meson Decays: SM Scheme 



As we discussed briefly in the Introduction, the main idea of the spin-averaged meson 
mass (SM) scheme |2^, which was originally argued by Dziembowski and 



Mankiewicz |^0|, is to take Mq in the spin-orbit wave function lZ^-^{x, k_|_) given by 
Eq. ( |1.7| ) as the average of physical masses with appropriate weighting factors from 
the spin degrees of freedom. For distinction from the invariant meson mass (IM) 
scheme (see Appendix A), we use the following covariant form of the spin-orbit wave 
function in the SM scheme[]: 



where the operators TM,^l are given by 



J 



PC 



0" 

r 

1^ 



V,fJ, 



a,fi 



{mp+ /')75, 



ma 



75- 



(2.1) 



(2.2) 
(2.3) 

(2.4) 



^See the Appendix of Ref. for more detailed derivations, where the Lepage-Brodsky conven- 
tion p^, EH] was used to obtain the Dirac spinors u and v given by Eq. (2.1). 
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Here mp{^)^a is the spin-averaged meson mass of the pseudoscalar (vector), axial vector 
mesons and the space components e{^) of the polarization four- vectors e{^) in the 
rest frame have the components e(±) = =f(1? ±^5 0)/v^, ^(O) = (0, 0, 1). 

In addition to the above treatment of the dynamical character of the angular mo- 
mentum operator in LF dynamics, we assume the radial wave function of the ground 



state mesons are described by the harmonic-oscillator (HO) wave functions [p5[] : 

2 k2, + 771,2 



(x, kj_) = Aexp 



E^^^"^/8/5' , (2.5) 

This Brodsky-Huang-Lepage (BHL) oscillator prescription [^, ^ of Eq. (|2.5| ) can 
be connected to the HO wave function in equal-t formulation by equating the off-shell 
propagator E = — (J2^=ikf) in the two frames |T^. As explained in Ref. p6| , 
the equal-t HO wave functions are known to give a reasonable first-approximation 
description of the static properties in a scheme without short-range hyperfine interac- 
tions. Therefore, it is reasonable to assume the spin-averaged mass as the meson mass 
in Eqs. (2.2)-(2.4). For example, the n and p masses are equal and both approximated 
by the spin-averaged value mu = {\^-k + |"^p)exp. 



As the SM scheme in LFQM P0| , |21[| provided a remarkably good description of the 
static properties for the pion and K mesons and reproduced the basic features of 
the lattice QCD and the QCD sum-rule results for vr, K, p, and Ai mesons, it is our 
intention to investigate more observables with the same model. In this work [p2| , 
we present a comprehensive study of the radiative decays of pseudoscalar (vr, K, t], t]'), 
vector {p, K* ,uj,(j)) and axial vector(Ai) mesons. We shall also contrast our results 
with the results based on a different treatment of meson mass such as the invariant 



mass (IM) scheme |jT^, in the calculation of the magnetic and quadrupole moments 
of the p meson, and the transition form factors of p ^ TC'y* and Ai — > 717*. The spin- 
averaged masses of it, K,p, and Ai mesons are given by 777^ = 777p = 0.612 GeV, 
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Table 2.1: Three different mixing schemes for rj and rj' and the corresponding spin- 
averaged masses. 



^5(7(3) 


Xri 








m^[MeV] 


m^/ [MeV] 


0° 










842 


885 


-10° 










843 


884 


-23° 


0.85 


0.53 


0.53 


0.85 


834 


873 



rriK = rriK* = 0.793 GeV and m^^ = 1.120 GeV Since we now consider r],r]',uj 
and mesons also in this work, we present the details of how to obtain the spin- 
averaged values of 1], rj', uj and mesons in Appendix B. The flavor assignment of r] 
and r]' mesons in the quark and antiquark basis are as follows: 



{uu + dd) _ 

V = 

, (uu + dd) _ 

V =X^' ^ + Yr,'SS, 



(2.6) 
(2.7) 



where = F^, = - sin 6p and = X^/ = cos 6p with 6p = 6su{3) - Oideai ~ 6su{3) - 



35°. Of particular interest are the values of mixing angle 6 



SU{3) 



-10° (so called 



"perfect mixing" ) ^ and —23° that are used in our analysis. The spin-averaged 
masses of 1] and rj' for each scheme are given in Table |2.1| . For uj and (p mesons, we 
use the scheme of ideal mixing |H9i and obtain the spin-averaged masses = 928 



MeV and = 799 MeV. Once the spin-averaged masses are fixed, then besides the 
well-known constituent quark masses of {u,d,s) quarks, i.e., = m^^ = 330 MeV 
and rus = 450 MeV, the only papameter in this model is the Gaussian parameter 
(3 which determines the broadness (or sharpness) of radial wave function. We will 
present our numerical results for a typical (3 value of /3 = 360 MeV throughout this 



section, unless stated otherwise. 



This Chapter is organized as follows: In Section 2.1, we calculate the form factors, 
charge radii, magnetic and quadrupole moments of the p and Ai mesons and compare 
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with the results of QCD sum rules ^0|. In Section 2.2, the transition form factors 
and the decay widths of the transitions V — > P7*, P — > V-y* (V = p,K*,u,(f) and 
P = IT, K, rj, rj') and Ai 717* are presented including the comparison with other 
theoretical results as well as the experimental data. In Section 2.3, we present the 
calculation of the transition form factors of the 7r° 7*7, rj 7*7 and f]' — 7*7 
transitions and compare our results with the recent experimental data [0, p3 
Summary and discussions of our major results follow in Section 2.4. 



2.1 The Form Factors of the p and Al Mesons 



Our analysis will be carried out using the Drell-Yan-West (g+ = 0) frame [^: 

P={P\p-,P^) = iP\—,0^), 

Q = {(t^q ,q±) = (0, -^,q±), 



(2i 



where M is the meson mass and the photon momentum is transverse to the direction 
of the incident spin-one system, with q\ = Q'^ = —<f- 



The Lorentz invariant electromagnetic form factors = 1, 2, 3) for a spin 1 particle 



are defined |54| by the matrix elements of the current operator between the initial 
|P, A) and the final |P', A') eigenstate of momentum P and helicity A as follows: 



{P',\'\r\PA) = step 



g"^(P + P'YF^iQ^) + (g'^^g'^ - g'^°g^)F2(Q' 



(2.9) 



+ q''q^{P + P'YF^{Q^)/{2M^) 

where q = P' — P and the polarization vectors of the initial and final mesons e = ex 
and e' = ey, respectively, are defined by 

5(0) 



5(±1) = ^(0,0,1,±^), e'(±l) = ^i(0,^,l,±0. (2.10) 
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Also, the Lorentz invariant form factors Fi{Q'^) are related to the charge, magnetic 
and quadrupole form factors of a meson pil as follows: 



Fc = F^+ '^kFq , 

Fm = F2 , (2.11) 
= Fi - F2 + (1 + k)F3 , 

where k = Q'^/AM'^ is a kinematic factor. At zero momentum transfer,these form 
factors are proportional to the usual static quantities of charge e, magnetic moment 
yUi, and quadrupole moment Qi. 

Fc(0) = l, eFA/(0) = 2M/ii, eFgiO) = M^Q^. (2.12) 

In the LFQM, the matrix element can be calculated by the convolution of initial and 
final LF wave functions of a meson: 

{P',\'\J^\P,\) = eJ rfa; /— 40(a;,kx + (l-x)q^)</>(a;,kx) 

JO J ibTi^ 

^ IIxl;,A^(a;,kx + (l-x)qx)r^XAi,A2(a;,k^) 

A, A' 

+ 62(1 ^ 2 of the first term), (2.13) 

where the spin-covariant functions xi^i^i) ^ire given by Eqs. (2.2)-(2.4) and the 
vertex is obtained from the expression \u{p') / {p'^Y/'^]'y^[u{p) / {p'^Y^'^] |^. Note 
that the jacobi factor dkn/dx given by Eq. (|l]^) is absorbed into the normalization 
constant A in Eq. ( p.5| ) throughout the Chapter 2. 

The relationship in Eq. (p.9|) between the covariant form factors and current matrix 
elements can be applied [BE, in principle, to any choice of Lorentz frame. As 



discussed by Brodsky and Hiller p%, in the standard LF frame [39|, i.e., = 



Qy = 0, and = Q, the three form factors can be obtained from the '+' component 
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of three helicity matrix elements: 



Fc-- 
Fm 

Fn -- 



2P+(2k+ 1) 
2 

2P+(2k + 1) 
1 

2P+(2k + 1) 



3 



(2k- 



2k 



-FoUTT(2/^-l)i^;- 



00 



00 



K+1 



(2.14) 



v/2k K 

where we defined {P' , \'\J^\P, X) = Fyy After a tedious but straightforward calcu- 
lation, we find the following expressions for the helicity form factors of the p meson: 



00 



dx 



x{l — x) 



exp 



+ ml + x{m2 — ml^'' 



x{l — x) 



2x^(1 - xf + x(l - x){aj + al- 2^^) + (aiaa)^ 
aj + al- 4aia2)^^ + ^'^ 



^+0 



Q 



,r dx 



2^2/3 

x[aia2 + - a;(l - x)], 



xfl — x) 



(ai - 02) 



Fl_ 



Q 



A(3 



1 (1 — x)dx 
I 1 aia2exp 







^ + ml + x{ifi\ — m\ 
xil — x) 



(2.15) 

(2.16) 
(2.17) 



where 



(i-x)^g^ 
^ 16/52 ' 

{xmp + mi) 



mi 



m,; 



a2 



2(3' 

((1 - x)mp + m2) 



(2.18) 



213 ' 2p 
and mi, m2 are the constituent masses of the quark and antiquark. The normalization 
constant Np = {AApj3^ /tcP'^)'^ is fixed by the definition of charge, Fc{0) = 1. For the 
systems of spin 1 or greater, in addition to the parity and time-reversal invariance of 
the current operator J^(0) an additional constraint on the current operator 



comes from the rotational covariance requirement [^, |58|, The angular condition 
for the spin-1 system is given by O, pBI, p3: 



[ 1 + 2k) F+ + F+ - VS^F 



— — n 
+0 -'^00 ~ 



(2.19) 
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As a matter of fact, the expressions of the right-hand side in Eq. ( |2.14| ) are not unique 



because of the angular condition in Eq. (|2.19|) . As pointed out in Refs. [[T^, p9 
unless the exact Poincare covariant current operator beyond one-body sector is used, 
the angular condition is in general violated (i.e., A(Q^) 7^ 0) and the calculation 



of the form factors Fi is dependent on the expressions on the r.h.s. of Eq. (|2l|). 



Examples of different choices of current combinations can be found in Ref. for 
the calculation of p meson form factors. As shown in Fig. 2.1(a), our result obtained 
for A((5^) is comparable with other choices given in Ref. IT] 



The magnetic (in unit of e/2M) and quadrupole moments (in unit of e/M^) of the p 
meson are obtained by Eq. (p.l2|) , 

pi=2.1, Qi = 0.41. (2.20) 

These values are not much different from the values of other model predictions based 
on the IM scheme presented in Ref. {pi = 2.3, Qi = 0.45) and Ref. ^ {pi = 2.26, 
Qi = 0.37). We also calculated the electromagnetic radii associated with the form 
factors, F^, F^j and as (4^)= 14 GeV-^, {rU= 22 GeV^^ (r|.^)= -5 
GeV-2, respectively. The results of Fc{Q^), Fm{Q^) and Fq{Q^) for < < 
5 GeV^ are shown in Fig. 2.1(b). To see the effect of the Melosh transformation 
(the measure of relativistic effects), we calculated the charge form factor by turning 
the Melosh rotations off and included the result in Fig. 2.1(b). As one can see in 
this figure, the breakdown of rotational covariance occurs around = 1 — 2 GeV^. 
Therefore, the model calculations of the form factor may not be reliable beyond 1 — 2 
GeV^ range. However, for the comparison with other model calculations, we have 
displayed our results beyond this region. The charge radius from this nonrelativistic 
form fator is ('"^)non-rel ~ ^'^'^ GeV"^ which is about 30% smaller than that of 
relativistic charge radius. In Fig. 2.1(c), we also compared our results with the 
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2.0 



1.5 - 




Figure 2.1: (a) A(Q^) testing the angular condition is shown as a function of . 
The sohd hues are our results with /3 = 0.32 and 0.36 GeV. For comparison, we in- 
clude various choices of the wave function introduced by Godfrey and Isgur(GI) 
in Ref. |T^. The dashed, dotted and dashed-dotted lines correspond to (spin- 
independent part), wgi (OGE + linear confining term) and Wcon/ (linear confinging 
term), respectively, (b) The form factors of p meson with parameter /? = 0.36 GeV. 
The solid, dotted, and dashed lines correspond to Fc{Q'^),Fm{Q'^) and Fq^Q"^), re- 
spectively. The dashed-dotted line is the result of non-relativistic limit of FciQ'^) by 
turning off the Melosh transformation. 
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Figure 2.1: (c) The form factors, Gi = Fi, G2 = F2 - Fi and G3 = F3/2, of p 



meson are compared with the results of QCD sum rules . The dotted, dashed and 
dashed-dotted lines correspond to Gi, G2, and G3 of Ref. . 



previous calculations of these form factors which were made in the framework of 
QCD sum rules by loffe and Smilga [pO| .^ 



The electromagnetic form factors of the Al meson are defined by Eq. (|2.9| ) as in the 
case of the p meson. Using a similar method taken in the p meson case, we find the 
following helicitiy component of the Al form factors analogous to Eqs. (2.15)-(2.17): 



00 



N 



Al 



dx 



xil — X, 



exp 



x(l — x) 



X 



6x\l - xY + 2x\l - xYlaf + ai- + x(l - x)[(aia2) + T 



-^\ciia2f - {a\ + al + 4aia2)^^ + 



(2.21) 



+0 



Q /-i dx 
T^Nai / —exp 
2p JO X 



x{l — x) 



^2 - aij 



^ The definition of the form factors(G'i) by QCD sum-rule |Q and the definition in this paper 
are related as follows: Gi = F\, G2 = F2 — Fi and G3 = F3/2. 
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2x^(1 - xf + x(l - x)(4^2 - aias) + - aiOa^^ 



exp 



.1 {l-x)dx 

W ""'Jo ~~x 

x(a2 - aif[x{l - x) 



_2 ^2 



x(l — x) 



(2.22) 



(2.23) 



As one can see in Fig. 2.2(a), the behavior of each form factor looks similar to the p 
meson case. In Fig. 2.2(a), we also included |A((5^)| for Ai. The amount of deviation 
from the rotational covariance is almost the same as that of the p meson case. 

The electromagnetic radii associated with F^^ and are obtained as {rp^) = 
10 GeV-2, {r^pj= 20 GeV-^ and (4^)= -10 GeV-^, respectively. The magnetic 
and quadrupole moments of the Ai meson are also obtained as 



/xi = 2.16, Qi = 1.34. 



(2.24) 



While the magnetic moment of Ai does not differ from that of the p meson, the 
quadrupole moment of Ai is about 4 times greater than that of the p meson in 
accordance with the fact that the Ai meson is a bound state with nonzero orbital 



angular momentum / = 1[14]. In Fig. 2.2(b), we compare our results with those of 
QCD sum rules] 



2.2 The Form Factors for V{P) P{V)Y and Ai 
7r7* Transitions 

The transition form factors oi A ^ B^* {{A, B) = (p, tt), (p, rj), {u, it), {uj, t]), {K*, K), 
{r]', p), {r]', u), (0, rj), (0, r]')) and Ai nj* are defined by 

{BiP')\r\AiP,X)) = eGAB(Q')e'^^"^£.(P, A)P^P;3, (2.25) 
(7r(P')l^1^i(^, A)) = ^\iV- qg>^' - V'-qnG.iQ') 
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Figure 2.2: (a) The form factors of Ai meson with the parameter P = 0.36 GeV. 
The sohd, dotted, and dashed hnes correspond to Fq, Fm, and Fq, respectively. The 
quantity |A((5^)| is shown as dashed-dotted hne. (b) The form factors, Gi = Fi, 
G2 = F2 — Fi and Gs = -F3/2, of Ai meson are compared with the result of QCD 
sum rules 0]. The solid, dotted, and dashed lines are our predictions of Gi, 6*2/2, 
and G3, respectively. The QCD sum-rule results are G'i((5^)(long dashed line) and 
G'2(Q2)/2 (dashed-dotted line) 
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+^{V ■ qq^ - q^V^)q''G,m\e,{P, A), (2.26) 

where e(P, A) denotes the polarization vector of the initial particles and V = P + P'. 
We used the same q^ = frame for the calculation of transition form factors as 
defined in Eq. ( p.8| ). As shown in Appendix C, the width of the decay A B'-f is 
given by 

where a is the fine-structure constant, Sa is the spin of the initial particle and Ma{b) 
is the mass| of the meson A (B). In the case of Ai — > tt'j transition, the decay width 
is expressed in terms of (^1(0) (see Appendix C), 



4a 

T{A, - vr7) = ^ 



Gi(0) 



2 /Ml. -M2\3 



2Ma, 



(2.28) 



Ma, 

In the calculations of the transition form factors V -P7*, we used the '+' component 
of the matrix elements of the current operator J'^. Since both sides of Eq. (p.25|) are 



vanishing for the longitudinal polarization Ap(^) = for any Q value, we considered 
Ap(a;) = 1 to calculate the transition form factor. Even though the matrix element J+ 
of this component (Ap(t^) = 1) also vanishes for real-photon(Q2 = Q) limit, this factor 
of = is nothing to do with the form factor itself because the r.h.s. of Eq. ( p.25| ) 



also vanishes in the limit of = 0. In other words, both sides of Eq. ( |2.25| ) in the 
case of Xp{Lo) = 1 have the same kinematic factor which does not vanish for any 
except = 0. Therefore, one can extract the transition form factor. We thus find 
the following expressions for the form factors Gp^iQ^) and Gujtt{Q'^)- 

Gpn{Q^) = {eu + ed)Ipvri™M,mg,P), 

G^.„{Q^) = {cu - ed)Ipvri^M,mg,P), (2.29) 



■'This must be the physical mass rather than the spin-averaged mass because the phase factor 
is nothing to do with the modeL Spin-averaged masses are used only for the calculation of form 
factors. 
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^PvAmM,m,,P) = ^^^l -exp[ . 



X 



a\ + a\ + + [a;(l — x) — i^^] + a\a}2{a\ + al^^ 



(2.30) 



and a}^^^ = (xm^^.^jj + mg)/2/? with i and f meaning incident and outgoing mesons. 
The normahzation constants of the p and vr-mesons are given by N.,^ = ^Np = 
2(2A^(,)/3V7rP+)2. 

We obtained the following prediction of decay widths of the p^(770) — > 77^7 and 
c<j(782) — > TT'-f transitions: 

r(p± ^ vr±7) = 69 keV (r;?_^,±^ = 68 ± 8 keV), 

r(c<; ^ 7r7) = 708 keV {VJZ^^ = 717 ± 51 keV). (2.31) 

Our results for the decay widths Tpj, and F^^ are in a very good agreement with the 
experimental data. The electromagnetic radii of these form factors are obtained as 
(^G'p^)= 7 GeV^^ and {rQ^J= 20 GeV^^, respectively. In Fig. 2.3(a), we present 
the transition form factor of p — > Tcy* 

for < g2 < 8 GeV^ and compare with 



other model predictions of Refs. [0, |6^. In Fig. 2.3(b), our prediction(solid line) 
of ^ 717* transition form factor in the spacelike region is compared with the 
VDM (dashed line) with F^^ = 1/(1 + Q^/M^) and the pole fit (dotted line) of the 
experimental data in the timelike region |^T|, |62| . 

Similarly, we calculated all other radiative decay processes between vector (1 ) and 
pseudoscalar mesons(0~"'") using the one loop integral formula /pv'^*(mM, rUg, [3) given 
by Eq. ( [2.30 ). We needed to change only the spin- averaged meson masses tum and 



constituent quark masses mj(i = 1,2), accordingly. Thus, for the p rj'j* and 
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Figure 2.3: (a) The transition form factor of — > 7r"'"7* multiplied by Q/" is shown 
as a function of . The solid lines are our results with /3 = 0.32 and 0.36 GeV. The 
dotted line, dashed, and dashed-dotted lines correspond to the results obtained using 



wgi in Ref. ||T^, the BS approach of Ref. |6y] and the predictions of VDM model 
with GpTT = 1/(1 + Q^/M^). (b) The normalized transition form factor oi lo ^ Tc-f*. 
Our result with f3 = 0.36 GeV (solid line) is compared with VDM (dashed line) with 
Guin = 1/(1 + Q'^/Mp) and the pole fit of the experimental data (dotted line) |]6l| . 
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uj — s> ?77* decays, we obtain the following transition form factors: 

Gm{Q'^) = X^{eu - ed)Ipvri'^M,mq,p), 

G^niQ'^) = Xri{eu + ed)Ipvri'^M,mg,p). (2.32) 

Using the two different mixing schemes, i.e., 6su{3) = —10° and —23°, we obtain the 
decay widths of the transitions p ^ rj^y and c<j — 777 as ^p}^^^ = 56 keV, T~}^^^ = 
6.4 keV and ^p^^-y = 65 keV, Tj^"^ = 7.4 keV, respectively. Both schemes are 
in excellent agreement with the experimental data of r°^^^ = 58 ± 10 keV and 
r^^^^ = 7.0 ±1.8 keV. The electromagnetic charge radii of p — > 777* transition are 
also predicted as (r^ ) = 16 GeV-^ for -10° and 18 GeV^^ _2^o 

mixing angle. 

The charge radii for u — 777* transition are (tq^ ) = 5 GeV-^ for -10° and 6 GeV"^ 
for —23°, respectively. 

The form factors of the 1]' —>■ p'j* and rj' — > uj'j* are given by 

Gri'piQ^) =X^/(e„ - ed)Ipvri^M,mq,P), 

Gri'ujiQ'^) = Xr^>{eu + ed)Ipvri^M,mq,P). (2.33) 

Our predictions of the decay widths are given by r~,l^J^^ = 117 keV, r~,l^J^^ = 9.7 keV 
and = 72 keV, r-,^^^^ = 6.0 keV. The experimental data of T^^^^^ = 61±8 keV 

and r°f^^^ = 5.9±0.9 keV. It is interesting to note that in case of transitions involving 
7]', the result of —23° mixing scheme is much better than that of —10° mixing scheme. 
The electromagnetic charge radii for —23° mixing angle are predicted as {tq , )= 12 
GeV^^ and {r^ , )= 4 GeV^^, respectively. 

The transitions of the K*^ — ^ K^'y* and — > K^'-f* in which the constituent quarks 
have unequal masses are also interesting processes to test our model predictions. As 
shown in Ref. the predictions for the kaon charge radius (r^Y^"^, the kaon form 
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Figure 2.4: The normalized form factors of K*^ ^^+7* (solid line) and K*^ 
i^°7* (dotted line) transitions with the parameter (3= 0.36 GeV. The dotted and 
dashed lines correspond to charged and neutral vector kaon transition form factors 
by Ref. The dashed-dotted line corresponds to the VDM prediction. 



factor Fk, and the decay constant fx are consistent with available experimental data. 
The transition form factors of the charged and netural K* decays are given by 

Gk*±k±{Q'^) = ±(e« + es)Ipvy*{mM,mg,p), 

Gk*oko{Q'^) = (cd + e,)Jpy^*(mM,mg,/?). (2.34) 

Using the spin-averaged masses of mK*±= 792.1 MeV and mK*o= 796.5 MeV and 
the constituent quark mass ms= 0.45 GeV [^], we obtain the decay widths for these 
charged and neutral vector kaon decay processes as TK*±^K±-y= 53 keV, TK*o-,KO'y= 
122 keV, respectively. The experimental data are r^?±^^±^= 50 ± 5 keV, r^?o^^o^= 
117 ± 10 keV. The electromagnetic transition charge radii are predicted as (r^ ,± ± ) = 
11 GeV~^ and ('^g^.o^q) = 10 GeV~^. In Fig. 2.4, we report the Q^-dependence 
of transition form factor of charged and neutral vector kaon, i.e., Gx*+k+{Q'^) and 
Gx-'OKoiQ'^), respectively, for < < 10 GeV^. Even though we showed the re- 
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suit of /5 = 0.36 GeV only, we note that our results become much closer to those of 



Bethe-Salpeter(BS) quark model prediction with smaller values of /5, i.e., (3 ~ 0.3 
GeV. 

For the decay processes of — > rjij]')'-)*, the transition form factors of — 777* and 
Vj'^* are given by 

G^n'iQ^) = Yn>2esIpvr i^M, mq, (3). (2.35) 

As shown in Appendix B, we obtain the spin-averaged mass of 0-meson as = 0.799 
GeV. Our predictions of the decay widths of the transitions — 777 and Vj^ 
with Qsu{Z) = -10° and -23° are given by r^^°^ = 61 keV, T^%^ = 0.28 keV and 
r^^^^ = 45 keV, r^^^,^ = 0.45 keV, respectively. The current experimental data are 

pexp _ e-fi o -U 9 O V^\r ^r^A F'^^P 



56.9 ± 2.9 keV and F^J;^,^ < 1.8 keV. It will be very interesting to compare 



our results, especially for r]''y, with the more precise measurements envisioned at 
TJNAF. The electromagnetic charge radii for —10° mixing angle are also predicted 



as (rf. ) = 5 GeV"' and (rj. ) = 5.5 GeV"'. 

In the case of transition Ai 717*, the kinematical factors in front of the form factor 
G2 in Eq. (2.26) yields zero for a certain in the standard Drell-Yan frame(Px = 
0, = 0) as discussed in Ref. [|T^]. This leads to a technical difficulty in extracting 



the form factors numerically from Eq. (2.26). Thus, as pointed out in Ref. fT^, we 
use the following symmetric coordinate frame for the calculation of of Ai n'-f* 
transition form factors: 



p+ ' 2^^;' V p+ '2" 

g=(o,^^i^,q.V (2.36) 



P+ 



30 



Then, we find the following expressions for the Ai form factors: 



P+(ml 



[mz 



Q 



2ml 







[nil — m. 



(2.37) 
(2.38) 



where = {7t{P')\J'^\Ai{P, A)) and its matrix elements are given by 

+ TTil + x(m2 — mf 



G'r 







exp 



2/3 V 2 Jo X 
x(ai -4)k(l -x) +^^] 
"1 (ix 



2iV^,iV, 



x(l — x) 



exp 



x(l — x) 

+ + x{m2 — ml^'' 



x{l — x) 



Q Na,K f^dx 

2 exp 



P 



x^ 



^ + mf + x{m\ 



7^2 

2 



x(l — x) 



(2.39) 
(2.40) 

(2.41) 



with a common facfor in the longitudinal component of G: 
Qo = 2x^(1 - xf[{a{ + 4) - (4 + 4)] 

+x(i — x)[4ct2(4 + 4) ~ 44 ('^1 + '^2)] ^ '^'^[(4 + 4) + ('^1 + '^2)] 
+^^[44(4 + 4) + 44(4 + 4)]- (2-42) 

In Fig. 2.5, we show the result of the Ai transition form factors, GilQ"^) (solid line) and 
G2{Q'^) (dashed line). While the form factor G2{Q'^) agrees with the predictions of the 
QCD sum-rules(dot-dashed line) in the region 

1 < g2 < 3 GeV^ the form factor 
GiiQ"^) seems to be quite different from the QCD sum-rule result (dotted line) p]| . 
However, there are other QCD sum rule calculations of Gi and G2 [|63l , |6^ and the 
results are rather different from each other. The authors of QCD sum-rules also 
pointed out that their predictions for the transition Ai TC'y* are of semiqualitative 
results. The electromagnetic radii corresponding to these form factors are evaluated 
as (r-Gi) = 15 GeV~^ and (r^^) = -30 GeV~^, respectively. 



The decay width of ^4^(1260) — > 7r"'"7 is obtained from Eq. (|2.28| ) as 

T{Ai 7r7) = 705 keV. (2.43) 
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Figure 2.5: The form factors of Af 7r+7* transition with the parameter f3 = 0.36 
GeV. G'i((5^) and G2(Q^)/10 of our results are represented by the sohd and dashed 
hues, respectively. For comparison, we show also the QCD sum-rule results |pO|| : 
Gi(Q2) (dotted-line) and G2(<5^)/10 (dotted-dashed line) 



The prediction of VDM is given by F^fit^ = 1000 - 1500 keV and the 

experimental value has been reported |6^ as T{Ai —>■ 7r+7) = 640 ±246 keV using 
the measurment of Primakoff production of the Ai resonance. Thus, our predicted 
width is quite consistent with the corresponding experimental data. In Refs. |T^, |T5|, 
the same decay width was calculated using the invariant mass of Ai meson instead of 
the spin-averaged mass. Their results F(y4^ — > 7r+7) = 319 keV [|T?| and 250 keV ||T3| 



seem to have rather large discrepancy from the experimental data. 
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2.3 The Form Factors for (77^,77,77') 7*7 Transi- 
tions 

The transition form factor of P ^ l*l{P = tt'', r], and r]') is defined from the matrix 
element of electromagnetic current = {■y{P + q)\J^\P{P)) as follows: 

= ie^Gp,{Q')e^,,^P''6Pq'', (2.44) 

where P and q are the momenta of the incident pseudoscalar meson and virtual 
photon, respectively, and e is the transverse polarization vector of the final (on-shell) 
photon satisfying e_L ■ q_L = 0. We again used the standard q~^ = frame for the 
calculation of the decays involving two photons, 

P=iP+,p-,P^) = {l,ml,,0^), g = (0,g2-m^,qx), (2.45) 

where P~^ is arbitrary but for simplicity we choose P~^ = 1 and thus we have = 
—q'j_ = —Q^. Each component of the final-state real photon is uniquely determined in 
this reference frame by the conservation of momentum and the on-mass shell condition 
of the real photon in the final state. 

12 



The Q -dependent decay rate for P — > 7*7 is given by [^] 

Tp_,,.,{Q') = ^a'M'pGl^iQ'), (2.46) 

where Mp is the physical mass of P = tt^,i], and rj'. Here, the decay width is given 
by Tp^^y at = 0. If we choose the '+' component of the current, the vertex factor 
in Eq. ( p.44 ) is given by 



vx'{x2,'ki_) ,ux{xi,k± + q_i_) ux{xi,k± + q±) ,uxixi,k±) 

■7 



X2 y/Xi y^Xl y/Xl 

1 



~ + q±y + m?]/xi - {k\ + m?)/x2 



+ (1^2) 



(2.47) 
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Here, only anti-parallel helicities of constituents contribute to the integrations and 
our model wave function \E'p(x, k^) for anti-parallel helicities is given by 



TT A^p («i«2 - k^ ) / k^ 

•A r, ) 



exp 



(2.48) 



2/?3^ 2 ' x{l-x) """'V 8x{l -x)l3\ 
A straightforward calculation for the transition P gives the following result: 



])Ip^*^{m^,mu(d), P), 



(2.49) 



where 



^ /-i dx 
Ur—— / —exp 
irp Jo X 



~ 1 
m 



X 



X 



2x{l-x) 

^2 



.2x{l — X 
4^2 + 



_|_ fjri^ 

2x{l — x) — exp ( -^^—j- ^ ) (aia2 + 'rh'^ + 4^^^ 



(it exp 



t /t 



2x{l-x) 

Using Eqs. ( p.46| )-( p30| ), the decay widths for P — 77 are obtained as 

r(7r° ^ 77) = 6.50 eV, r(r/ 77) = 0.47(0.65) keV, 
Y(r]' 77) = 7.9(5.6) keV, 



(2.50) 



(2.51) 



where the values for r/ — > 77 and rj' —>■ 77 are obtained for the —10° (—23°) mixing 



scheme. The experimental data|^9| are given by F^^p = 7.8±0.5 eV, = 0.47±0.05 
keV and = 4.3 ± 0.6 keV. The agreement of our results with the experimental 
data is not unreasonable. In Figs. 2.6-2.8, the Q^-dependence of the decay rate 
Fp^((5^) for P = 7T^,ri, and rj' are shown and compared with the recent experimental 



data |5T|, Our predictions for all of these processes are overall in a good 

agreement with the experimental data up to a rather large Q^. 
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Figure 2.7: The decay rate for rj 7*7 transition with the parameter f3 = 0.36 GeV. 
The sohd and dotted hnes correspond to the 9su{3) = —10° and —23° mixing schemes, 
respectively. Data are taken from Refs. pTl E2l . 
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Figure 2.8: The decay rate for rj' 7*7 transition with the parameter /3 = 0.36 GeV. 



The sohd and dotted hnes correspond to the 9 



5(7(3) 



— 10° and —23° mixing schemes, 



respectively. Data are taken from Refs. pTl E2l R3 



2.4 Summary and Discussion 



In this work, we have investigated the radiative decays of pseudoscalar(7r, i^, 77, 77'), 
vector(p, K*,uj, 0) and axial vector(Ai) mesons as well as the form factors of p and Ai 
mesons using a simple relativistic constituent quark model. We summarized all of our 
predictions on the meson decay widths for various values of (3 {(3 = 0.32, 0.34, 0.36, 0.38 
GeV) in Table Remarkably, most of our predictions with the parameter j3 in this 
region are within the experimental errors. We have also investigated the sensitivity 
of our results by varying quark masses. For ±10% variation of the non-strange quark 
masses, the decay width F and the transition charge radius (r^)^''^ of the radiative 
meson decays change by 3-4% and 2-4%, respectively. Changing strange quark mass 
by ±10% for the process of K*^ i^^7 yields 13-15% and 3-4% difference in the 
decay width and the transition charge radius, respectively. 
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Table 2.2: Radiative decay widths for the V{P) P{V)j, Ai and P ^ 2j 

for various model parameters j3 (in unit of GeV) with the rj-rj' mixing angle 6su(3) = 
— 10° [—23°]. The experiment is the result from unless otherwise noted. The 

unit of decay width is [keV], unless otherwise noted. 



Widths (T) 


13 = 0.32 


0.34 


0.36 


0.38 


Expt. 


H'^ — > TT'^'^V 


78 


73 


69 


64 


60 ± 8 


UJ — > TT'J 


775 


742 


708 


674 


717 ± 51 




60 


57 


53 


50 


50 ±5 




134 


128 


122 


116 


117 ± 10 




66 [77] 


60 [70] 


56 [65] 


51 [60] 


58± 10 


UJ ^ T]^ 


7.4 [8.5] 


6.9 [7.9] 


6.4 [7.4] 


6.0 [6.8] 


7.0 ± 1.8 


7]' p7 


137 [89] 


126 [80] 


117 [72] 


108 [66] 


61 ±8 


r]' — > 


11.2 [7.3] 


10.4 [6.6] 


9.7 [6.0] 


9.1 [5.6] 


5.9 ±0.9 




54 [40] 


58 [42] 


61 [45] 


65 [47] 


56.9 ±2.9 




0.26 [0.43] 


0.27 [0.44] 


0.28 [0.45] 


0.29 [0.46] 


< 1.8 


A\ — > 7r7 


620 


664 


705 


742 


640 ± 2460' 


TT^ ^ 27 


7.58 


7.06 


6.50 


5.91 


7.8 ± 0.5 [eV] 


77 ^ 27 


0.61 [0.78] 


0.53 [0.71] 


0.47 [0.65] 


0.42 [0.58] 


0.47 ±0.05 


7?' ^ 27 


8.8 [6.5] 


8.3 [6.1] 


7.9 [5.6] 


7.3 [5.1] 


4.3 ±0.6 



In our point of view, the success of this model hinges upon the simplicity of the light- 



front vacuum. The recent lattice QCD results [^ indicate that the mass difference 
between 77' and pseudoscalar octet mesons due to the complicated nontrivial vacuum 
effect increases(or decreases) as the quark mass decreases(or increases), i.e., the 
effect of the topological charge contribution should be small as increases. This 
makes us believe that the complicated nontrivial vacuum effect can be traded off 
by the constituent quark masses. This may mean that there is a suppression of 
complicate zero-mode contribution from the light-front vacuum in our model due 
to the rather large constituent quark masses. 



Our approach in this work was to model the wave fuction rather than to model the 
potential. However, we have attempted to compare our results with various other 
available theoretical results including the potential models and the QCD sum rules. 
At the very least, our results seem to be quite comparable with the results from 
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modelling the potential. The results on the angular condition are also not drastically 
different from the result of the potential models (See Fig. 2.1(a)). Furthermore, 
the agreement with the QCD sum rule results is not unreasonable. Our model has a 
predictive power and more experimentally measurable quantities should be calculated 
and compared with data. A particularly interesting prediction from our model is the 
branching ratio of — > rj^ estimated as 10~"^ for the 77-77' mixing angle of —23°. For 
the mixing angle —10°, the estimation for this branching ratio is reduced by 60% 
from that of —23°. Thus, it will be very interesting to compare our results with the 
precise measurements envisioned at TJNAF. 
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Chapter 3 

Relations Among the Light-Front 
Quark Models and the Model 
Prediction of 77 — 77^ Mixing Angle: 
IM Scheme 



One of the popular quark models in the LF formalism is the invariant meson mass 



(IM) scheme |10, 11, 17, 18, 19 in which the IM square Mq is given by 



E 



(3.1) 



The corresponding spin-orbit wave functions TZ'l'^'\^{x, k^) of pseudoscalar and vector 
mesons are obtained by the Melosh transformation and their expressions are given in 
the Appendix A. 



While the expressions of the spin-orbit wave functions for pseudoscalar and vector 
mesons are unique, the models differ in choosing the radial wave function. For exam- 
ples, Huang, Ma, and Shen adopted the Brodsky- Huang- Lepage (BHL) oscillator 
prescription [^, ^ (model H): 



(j)^{k^) = Aexp 



i = l 



Xi 



(3.2) 
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On the other hand, the authors of Ref. (model C) and Ref. (model J) used 
the following radial wave functions: 

0^(^)(fc2) = N^^j^ exp(-kV2/?^), (3.3) 
where k = (A;„, k^) is the three momentum and the normalization constants Nc{j) are 



Nc = {4:/ y/nP^y/'^ and Nj = ny2/3Nc, respectively. In addition to the harmonic 
oscillator (HO) wave functions, Schlumpf |]19[ used a power-law (PL) wave function 



given by 

0^^(A;2) = iVpower(l + kV/5')-^ (3.4) 

where the power s may be chosen typically as s=2 (model PL). Furthermore, more 
realistic model wave function from the Godfrey-Isgur potential was investigated 



in Ref. 17 . 



Regarding the differences in choosing HO radial wave functions, Huang, Ma, and 



Shen [18| pointed out the ambiguity of introducing a factor, such as J Mq/Ax{1 — x) 



in Ref. [rT| or ^l/2x{l — x) in Ref. [^, to the BHL wave function as a consequence 
of the jacobian relating the instant-form momentum to the LF momentum. This 
ambiguity which we will discuss in detail in the following section distinguishes the 
model H from the models C and J. 

The purpose of this work is to (1) show explicitly the difference of the model H from 
the other two HO models {C and J), (2) point out nevertheless that the numerical 
predictions from all of these HO models {H,J and C) are almost equivalent once the 
best fit parameters(i.e., constituent quark masses and Gaussian parameters (3) are 
chosen, (3) apply these models to the two photon decay processes and the radiative 
decays between pseudoscalar (P) and vector (V^) mesons, and (4) obtain the best value 
of rj-f]' mixing angle. In addition, we change the form of the radial wave function from 
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HO to PL (model PL) and investigate the presence-absence effect of the Jacobi factor 
to further identify the differences in the model wave functions. 

This Chapter is organized as follows: In Secion 3.1, we analyze the relations among 
the three HO model wave functions and show that they are basically equivalent with 
each other modulo a Jacobi factor. In Section 3.2, we list several constraints to fix 
the free parameters of the relativistic quark model in comparison with the relavant 
experimental data for the decay constants. In our analysis we consider both HO 
and PL wave functions and restrict ourselves on the light-meson sector {u-, d- and 
s-quarks) with equal quark and anti-quark masses(mg = mq). In Section 3.3, we 
analyze the two-photon decays of tt,-/], and rj' and the transitions of V{P) —>■ P(\^)7 
searching for the best value of rj-r]' mixing angles. Summary and discussion are 
followed in Section 3.4. 



3.1 Relations among the Light-Front Quark Mod- 
els 



Since the spin-orbit wave functions 7l'l'^\^{x, k^) obtained from the Melosh rotation 
are independent from the radial wave functions as far as we consider the same IM 
scheme, the relations among the three HO models {H, J and C) are obtained by 
comparing only the radial wave functions in Eqs. (p. 21) and 



The relation between (p {k ) ||TT| and (f) {k ) [|T^ can be obtained by the definition of 



the normalization of the wave functions 



^ I d'kl^p^i^)]' = (^/ d'k\<p\e)\' = 1, (3.5) 
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which leads to 



<j>'{k')=nJ-<P^{e), (3.6) 



where nc{= 3) in Eq. ( |3.5|) is the color factor. Eq. ( |3.6|) implies that the two 
models C and J are equivalent to each other. Comparing (f)^(k'^) and (f)'^^'^\k'^), 
we use the jacobian J' = dkn/dx = Mo/4x(l — x) of the variable transformation 
{x,k^} ^ k = (A;„,k^) in which the longitudinal component kn is defined as a 
function of x and k^ by A;„ = (x — |)Mo ||Tl[. Then the phase space [cPk] in Eq. ( |3.5| ) 
is transformed into (a;, k^) variables as follows: 

d^k = dx£k^ ^ (3.7) 
4x(l — X) 

The wave function (j)^ {k"^) can then be compared directly with (p^^k"^) (or (j)''{k'^)) in 
(x, k^)-space by the definition of the normalization of the wave functions 



From Eqs. ( p^ ) and (|37q ), we obtain 



Here, we should note that in order to satisfy the above equivalence relation the nor- 
malization constant A in Eq. ( ^.2| ) is no longer a constant but a function of the LF 
variables (x, k^): 



/ x{l — x) ' 









A{x,k^) = vrW I _^ )'^%xp(mV2/3^), (3.10) 



for {/)'^(fc^). However, in the model H [jT8[, A is not a function of {x, k_|_) given by Eq. 
( p.lOj ) but a constant. Therefore, the predictions of the model H are different from 
those of the models C and J. 
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With the analysis of the wave functions of each model, we have now shown that the 
two models C and J are completely equivalent to each other and the model H is 
different from C and J unless A satisfies Eq. ( ^.lUI ). If we were to use A{x, k^) given 
by Eq. (|3.10|) , then we can explicitly show the equivalence of all the observables 



among the three models. Similarly, one can consider the presence-absence effect of 
the Jacobi factor for the PL model wave function. In the next section, we will present 
the constraints to fix the free parameters in both HO and PL wave functions. 

3.2 Constraints on the Wave Functions of + and 
1 Mesons 

In order to fix the free parameters(mq and jSgq) of the models H and PL, we use the 
method adopted by Ref . |T^ (model J) . Considering the mesons with equal quark and 
anti-quark masses, we need four parameters in the u—,d— and s-quark sector, i.e., 
nT'u(d),^s, and Pqq{q = u{d),s). We fix these four parameters by fitting the pion decay 
constant /,r and the vector meson decay constants fviV = P,^^ and 0). The decay 
constant fp of a pseudoscalar meson P{qiq2) defined by (0|g27^75Q'i|-P) = iV^fpP^ 
can be evaluated as follows[| 



where ^ = (1 — x)mi + x'm2. Likewise, the decay constant fv of a vector meson is 
defined by {0\q2'j^qi\V) = Myfv^^- For the decay of p° = (uu — dd)/\^, we obtain 



/p = /y(mi,m2,/3)(e„ - ed)/V2, (3.12) 



^For convenience, we use the wave function of the model C ||Tl|, i.e., iff [k"^) = 
(4/V^/33)i/2 exp(-/fcV2/32). 
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Table 3.1: Best fit quark masses and the model parameters (3. PLJ and PLH are 
the power law(PL) models for s=2 in Eq. (3.4) with and without Jacobi factor, 
respectively. q=u and d. 





H 


J 


PLH 


PLJ 


mg[GeV] 


0.25 


0.25 


0.28 


0.28 


m,[GeV] 


0.37 


0.37 


0.37 


0.37 


/?,,-[GeV] 


0.36 


0.3194 


0.40 


0.307 


/3.i[GeV] 


0.38 


0.3478 


0.415 


0.328 



where0 

/y(mi,m2,/3) 



V6 

47]-2 



J dx J d'^kj_] 



A + 



2k2 



Mq + mi + m2 



(3.13) 



9x /42 + ki 

While the vr and p decay const ant s(/7r, fp) are used to fix the parameters m = rriu = 
rrid and /? = Pun = Pdd: the and decay constants(/a;, Z^) are used to fix and 
(3ss- Since the flavor mixing of u and </> in the quark basis is given by 



</) = — sin Sv {uu + dd) / \p2 — cos byss, 
uj = COS 6v{uu + dd) / V2 — sin Syss, 



(3.14) 



where Sy = Osu{3) — 35.26°, we use the vector mixing angle as Sy = —3.3° for the 
direct comparison with the results^ of the decay constants and in Ref. [|l^ . 



The constituent quark masses and the model parameters (3uu = Pdd ^ind (3ss fixed 



by the available experimental data for the decay constants f^^, fp, and are 
summarized in Table 3.1 for both HO and PL models. The PL models with and 
without Jacobi factor are denoted by PLJ and PLH, respectively. In Table 3.2, we 
summarize the results of the model H{PLH) for the decay constants and compare 

these results with those of the model J{PLJ) and the experimental data. In order 

^We have given the general result for unequal quark and antiquark masses since we shall use it 
to calculate, e.g., the decay constants fx*, fo*, and fs* etc. in the following chapters. 

■^The frf, and /^^ are obtained as f^=— sin Syly {mi, m2, 0){eu + ed)/V2 — cos 6yly{mi, m2, /3)es 
and fuj= cos 6yly{mi, m2, P){eu + ed)/\/2 — sin 6yly {mi, m2, [3)es, respectively. 
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Table 3.2: Decay constants for n ^ fiu and V —>■ e^e", where the mixing angle 
6v = —3.3° is used. H\PLH') uses the parameters of J{PLJ) to see the effect of 
the presence-absence of the Jacobi factor. The experimental data are taken from 
Ref. [E9l. 



fpiv) 


H' 


H 


J 


PLH' 


PLH 


PLJ 


Expt.[MeV] 


fn 


86.5 


92.1 


92.4 


79.2 


92.8 


92.5 


92.4±0.25 


fp 


138.5 


156.1 


151.9 


120.1 


155.5 


154.7 


152.8± 3.6 


fuj 


42.0 


47.5 


46.1 


43.4 


47.3 


47.1 


45.9± 0.7 


f<t> 


73.0 


80.2 


79.7 


62.3 


79.4 


79.6 


79.1± 1.3 



to see the effect of the presence-absence of the Jacobi factor given by Eq. (pT^), we 
present the results of the model H{PLH) using the same parameters of the model 
J {PLJ), which we call H'{PLH'), as well as the results of the model H{PLH) using 
the best fit parameters. As one can see from the comparison between H'{PLH') and 
J{PLJ) in Table 3.2, the effect of the Jacobi factor is not negligible. However, once 
the best fit parameters for the model H{PLH) are chosen as shown in Table 3.1, the 
numerical results of the decay constants are almost equivalent between the models 
H{PLH) and J{PLJ) regardless of the Jacobi factor. 



The quark distribution amplitude defined by 4>M,\{,Xi,Q)=^'°^[cPki]^! 



i.e., the probability amplitude for finding quarks in the Lz=Q (s-wave) projection 
of the wave function collinear up to the scale Q pT| , and the form factor of 
the pion from the models H'{PLH'), H{PLH) and J {PLJ) are plotted in Figs. 
3.1(a)[3.1(b)] and 3.2(a)[3.2(b)], respectively. As shown in Fig. 3.1(a)[3.1(b)], a 
rather substantial difference between and (j)'l{(j)^^^^is shown due to the 

presence- absence of Jacobi factor. Because the quark distribution amplitude depicts 
just functional dependence of the model wave function, the deviation in the quark 
distribution amplitude due to the presence-absence of Jacobi factor does not diminish 



''Because of the presence of the damping factor in Eqs. (3_^)-(3^), we extended the integral Hmit 
to infinity without loss of accuracy. 



45 



substantially even though the best fit parameters are used as one can see from the 
comparison between and in Fig. 3.1(a)[3.1(b)]. However, the 

deviation between model H{PLH) and model J{PLJ) is reduced more substantially 
in the physical observable such as the form factor of pion in Fig. 3.2(a) [3.2(b)] once the 
best fit parameters are used. This observation is more pronounced in the calculation 
of decay widths as we present in the next section. 



3.3 Decay Widths for P ^ 77 and V{P) P(y)j 
Transitions 



Applying these models to calculate the decay widths for P lli^P = T^,rj,rj') and 
V{P) P{V)'y(y = p, u, (p) transitions, the decay width for P ^ 77 is given by 



4 



2„2 



a gp^^rup 



(3.15) 



where the coupling constants gp^^ are expressed in terms of three pseudoscalar decay 
constants obtained by the axial-vector anomaly plus PCAC (partial conservation of 



axial- vector current) [|T0], |7^, |74 

1 



5'-7r77 
fl'»777 



gv' 



1 
1 



cos 6 



2V2 



SU{3) 



77 



fo 

— sm 9su{3) + —r- 
Js Jo 



sin 9 



SUi3) 



(3.16) 



cos 6' 



5(7(3) 



Here, the predictions of the decay constants fs and fo of the model H{PLH) are 
given by 

fs/f^ = 1.148(1.096), /o//. = 1.074(1.048). (3.17) 

These results are very close to the predictions of the model J{PLJ) |]lOl which are 
given by /g//^ = 1.156(1.094) and /o//^ = 1.078(1.047). The predictions of the 
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Figure 3.1: (a) The normalized quark distribution amplitude of the pion in HO model. 
The solid, dotted, dot-dashed lines represent the results of the models H, H', and J, 
respectively. The normalization is fixed by the zeroth moment (the area underneath 
each curve) (^°)= J^i d^^^(j).„{^)=l with ^=xi—X2- (b) The normalized quark distribu- 
tion amplitude of the pion in PL model. The solid, dotted, dot-dashed lines represent 
the results of the models PLH, PLH', and PLJ, respectively. The long-dashed line 



represents the Schlumpf's result ||T9 
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0.6 - 



0.5 



0.00 0.02 0.04 0.06 0.08 0.10 0.12 
Q'[GeV'] 




0.5 



0.00 0.02 0.04 0.06 0.08 0.10 0.12 
Q'[GeV'] 



Figure 3.2: (a) The square of the pion charge form factor in HO model for low values 
of compared with the experimental data The same line code as in Fig. 3.1(a) 
is used, (b) The square of the pion charge form factor in PL model. The same line 
code as in Fig. 3.1(b) is used. 
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Table 3.3: Radiative decay widths r(P —>■ '-ff) for rj-rj' mixing angle 6 



SU{3) 



-19" 



Widths 


H 


J 


PLH 


PLJ 


Experiment 


r(7r 77) 


7.79 


7.73 


7.67 


7.72 


7.8±0.5[eV] 


T{ri^ 77) 


0.49 


0.485 


0.52 


0.52 


0.47 ± 0.05 [keV] 


Y{r]' 77) 


4.51 


4.45 


4.64 


4.68 


4.3±0.6[keV] 



chiral perturbation theory in Ref. ||7^ has been reported as /s/Ztt = 1.25 and /o/Ztt = 
1.04 ± 0.04. We have used the rj-i]' mixing (see Eqs. (2.6)-(2.7)) analogous to uj-(j) 
mixing in Eq. ( |3.14| ) and observed that 9su{3) ~ — 19° gives a good agreement with 
the experimental data of decay rates. The decay rate and the transition form factors 
for P —>■ '~f*'~f {P=7r,rj, and f]') are given by Eqs. (2.46) and (2.49), respectively. 
However, the integral formula l2{mM,fnq, (3) in Eq. (2.50) should be changed into 
the following form: 



dx I d k\ 



xinir, 



4772 Jo J " dx , I ml + k2 k'^ + m; 



(3.18) 



where k'j_ = k_L + (1 — x)q_L. For Osuiz) = —19°, the numerical results of the decay 
widths are summarized in Table 3.3. As shown in Table 3.3, the predictions of decay 
widths for P ^ 77 are not much dependent on the models once the best fit parameters 
are obtained and the model predictions with 6su{3) = —19° are in excellent agreement 
with the experimental data. 



To further justify our observation on the Jacobi factor and the model predictions of 
rj-rj' mixing angle, we have calculated also the decay widths for V{P) P{V)'~f. We 
can also use the same formula for the decay width given by Eq. (2.27) but the integral 
formula Jpyy (mM, ""^g, /?) in Eq. (2.30) should be changed into the following form^: 

fdh\ |0^(A;')|' 



(mq,/?) 



2tt Jo 



d'ki 



\ dx J x[Mq — (mi — 7712)^] 



Even though wc considered the equal quark and antiquark masses in this work, nevertheless, 
we have given the general results for the unequal quark and antiquark masses since we shall use it 
to calculate, e.g., the decay widths for K*'^ — > K'^'^ and K*'^ — > K'^'j in the next chapter. 
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Table 3.4: Radiative decay widths for V{P) — > P{V)'y transitions. We use 77-77' and 
u}-(f) mixing angles as 6su{3) = —19° and Sy = —3.3°, respectively. 



Widths 


H 


J 


PLH 


PL, J 


Exi)eiinient [keV] 


r(p^ — * n^^) 


75 


76 


89 


97 


68 ±8 


r{uj TT'-f) 


712 


730 


855 


938 


717 ±51 


r(0 ^ 7r7) 


5.5 


5.6 


6.6 


7.2 


5.8 ±0.6 


r(p ^ 777) 


59 


59 


69 


76 


58 ±10 


T{uj^ 777) 


8.6 


8.7 


10.2 


11.1 


7.0 ± 1.8 


r(0 ^ 777) 


55.9 


55.3 


72.4 


74.2 


56.9 ±2.9 




66.1 


67.5 


79.1 


86.8 


61 ±8 


r(V ^ ^7) 


4.7 


4.8 


5.5 


6.1 


5.9 ±0.9 


T{(j) ri'^y) 


0.56 


0.57 


0.71 


0.76 


< 1.8 



X 



A + 



(3.19) 



Mo + mi+ 777,2- 

The results are summarized in Table 3.4, where the a;-0 mixing angle Sy — —3.3° is 
also taken into account. Again, the predictions of the model H{PLH) with 9su{3) — 
—19° are similar to those of the model J {PL J). Thus, we confirm again that the 
numerical results of the model LL[PLLL) and the model J [PL J) are not much different 
from each other regardless of the Jacobi factor once the best fit parameters are used. 
However, the difference in the choice of radial wave function(e.g. HO wave function 
versus PL wave function) is still appreciable even though the best fit parameters are 
used. The experimental data of decay widths both for P — > 77 and V{P) P{V)^ 
are quite consistent with the HO model predictions for 9su{z) ~ —19°. 



3.4 Summary and Discussion 



In this work, we have first shown that the HO models C and J are actually equivalent, 
while the model H is not exactly same with the other models C and J due to the 
Jacobi factor sJMq/^x{1- x). The effect of the presence-absence of the Jacobi factor 
is in principle not negligible as we have shown in Table 3.2 and Figs. 3.1 and 3.2 
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[comparison between H'{PLH') and J {PL J)]. However, once the best fit parameters 
are used, the numerical resuhs of the physical observables from the model H{PLH) 
are almost equivalent to those of the model J{PLJ) regardless of the presence-absence 
of the Jacobi factor as shown in the calculations of the decay const ant s(/7r, fp, fu, and 
fcf,), the pion form factor and the decay rates of V{P) P(V)7 and P — > 77 
transitions. 

In the case of best fit, the effect of the presence-absence of the Jacobi factor amounts 
to the difference in the best fit parameters as shown in Table 3.1. Thus, the effect from 
the Jacobi factor is there no matter what we do. However, what we have shown in 
this work indicates that it is rather difficult to pin down a better model in the present 
phenomenology if the two models (e.g., H and J) are differ only by the Jacobi factor. 
Nevertheless, the difference between the HO and PL results are substantial enough to 
state that the overall agreement with the data of various radiative decay widths (see 
Table 3.4) is clearly better in the HO models than in the PL models. We also found 
an excellent agreement of all three HO mode\{H, C and J) predictions on the P — 77 



and V{P) Piy)^ processes with the experimental data, if we use Eq. ( |3.16|) in 
two photon decay processes obtained by the axial anomaly and PCAC relations and 
choose the rj — t]' mixing angle 6su(3) ~ —19°. 
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Chapter 4 

Mixing Angles and 
Electromagnetic Properties of 
Ground State Pseudoscalar and 
Vector Meson Nonets in LFQM 



It has been realized that relativistic effects are crucial to describe the low-lying 



hadrons made of u, d, and s quarks and antiquarks |4^. The LFQM takes advan- 
tage of the equal LF time {T=t+z/c) quantization and includes important relativistic 
effects in the hadronic wave functions. The distinct features of the LF equal-r quan- 
tization compared to the ordinary equal-t quantization may be summarized as the 
suppression of vacuum fluctuations with the decoupling of complicated zero modes 
and the conversion of the dynamical problem from boost to rotation. Taking advan- 
tage of this LF approach, we have shown in Chapters 2 and 3 that once the best fit 
parameters were used, both SM and IM schemes in the LFQM provided a remark- 
ably good agreement with the available experimental data for form factors, decay 
constants and charge radii etc. of various light pseudoscalar and vector mesons as 
well as their radiative decay widths. 



52 



However, the radial function so far has been taken as a model wave function rather 
than as a solution of the QCD-motivated dynamical equation. Because of this, we 
could not calculate the mass spectra of mesons directly from the model but assume 
a couple of schemes such as SM or IM schemes. Even though the authors of Ref. |17 



adopted the quark potential model developed by Godfrey and Isgur to reproduce 
the meson mass spectra, their model predictions included neither the mixing angles of 
uj-(f) and 77-77' nor the form factors for various radiative decay processes of pseudoscalar 
and vector mesons. 

In this work, we are not taking exactly the same quark potential developed by Godfrey 
and Isgur . However, we attempt to fill this gap between the model wave function 
and the QCD-motivated potential, which includes not only the Coulomb plus confin- 
ing potential but also the hyperfine interaction, to obtain the correct p-vr splitting. 
For the confining potential, we take a (1) harmonic oscillator (HO) potential and (2) 
linear potential and compare the numerical results for these two cases. We use the 
variational principle to solve the equation of motion. Accordingly, our analysis covers 
the mass spectra of light pseudoscalar {n, K,!],!]') and vector {p, K* ,uj, (j)) mesons 
and the mixing angles of uj-(f) and rj-r]' as well as other observables such as charge 
radii, decay constants, radiative decay widths, etc. We exploit the invariant meson 
mass scheme in this model. We also adopt the parametrization to incorporate the 
quark-annihilation diagrams |7^, mediated by gluon exchanges and the SU(3) 
symmetry breaking, i.e., 7^ mg, in the determination of meson mixing angles. 

This Chapter is organized as follows: In Section 4.1, we set up a simple QCD mo- 
tivated effective Hamiltonian and use the Gaussian radial wave function as a trial 
function of the variational principle. We find the optimum values of the model pa- 
rameters, quark masses {mu(d),^s) and Gaussian parameters ( Puu=Pud=Pddj Pus, Pss) 
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for the two cases of confining potentials (1) and (2). We also analyze the meson mass 
spectra and predict the mixing angles of uj-cf) and rj-rj' . We adopt a formulation to 
incorporate the quark- annihilation diagrams and the effect of SU(3) symmetry break- 
ing on the meson mixing angles. In Section 4.2, we calculate the decay constants, 
charge radii, form factors, and radiative decay rates of various light pseudoscalar and 
vector mesons and discuss the numerical results of the two confining potentials (1) 
and (2) in comparison with the available experimental data. A summary and discus- 
sions follow in Section 4.3. The details of fixing the model parameters and the mixing 
angle formulations are presented in Appendices D and E, respectively. 



4.1 Model Description 



The QCD-motivated effective Hamiltonian for a description of the meson mass spectra 
is given by Jl^ ^ 



ml + k^ + J ml + fc2 + V,, 



I nlml 



qq 



"q > < •'qq 

\^ nlml ^^■^qq\^ nlml 1 



(4.1) 



where M^q is the mass of the meson, A;^ = + /c^, and Y^nim) the meson wave 
function given in Eq. (1.5) 

In this work, we use the two interaction potentials Vqq for the pseudoscalar (0""*") and 
vector (1 ) mesons: (1) Coulomb plus HO, and (2) Coulomb plus linear confining 
potentials. In addition, the hyperfine interaction, which is essential to distinguish 
vector from pseudoscalar mesons, is included for both cases, viz., 

4k 2S ■ S- 

Vqq = Voir) + V^y,{r) = a + V,onf - ^ + ^^^VVcoui, (4.2) 



where Vconf = hrir"^) for the linear (HO) potential and {Sq ■ S'g)=l/4 (—3/4) for the 
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vector (pseudoscalar) meson. Even though more reahstic solution of Eq. ( [4.1| ) can be 
obtained by expanding the radial function 0n,z=o(^^) onto a truncated set of HO basis 
states [T^, i.e., YJn!=i Cn(t>n,o{k'^) , our intention in this work is to explore only the 
and 1 ground state meson properties. Therefore, we use the IS" state harmonic 
wave function (pioik"^) as a trial function of the variational principle 

1 \l/2 



no{x, kx) = (^-^ ) exp{-e/2f5'), (4.3) 
where 0(a;, k_|_) is normalized according to 



2 



'^dx I rfV(^)|0io(x,k^)p = 1. (4.4) 

Because of this rather simple trial function, our results could be regarded as crude 
approximations. However, we note that this choice is consistent with the LFQM wave 
function which has been quite successful in describing various meson properties |jlO 



TT], [T2|, Furthermore, Eq. (|4.3| ) takes the same form as the ground state 

solution of the HO potential even though it is not the exact solution for the linear 
potential case. As we show in Appendix D, after fixing the parameters a, b, and k, 
the Coulomb plus HO potential Vo(r) in Eq. ( |4.2| ) turns out to be very similar in the 
relevant range of potential (r < 2 fm) to the Coulomb plus linear confining potentials 
[see Figs. 4.1(a) and 4.1(b)] which are frequently used in the literature fl^ 



|78| , |79| , pO| , pl| . The details of fixing the parameters of our model, i.e., quark masses 
{fnu{d),fns), Gaussian parameters {Pud: Pus, Pss), and potential parameters {a,b,K) in 
Vqq given by Eq. ( ^.2| ), are summarized in Appendix D. 



Following the procedure listed in Appendix D, our optimized model parameters 
are given in Table 4.1. In fixing all of these parameters, the variational principle 
[Eq. ( p.l| )] plays the crucial role for ud, us, and ss meson systems to share the same 



55 



> 

CD 



- 









, , , , , 


/ ^' 






H.O. 


//' 






Linear 


'; ' 






ISGW2(c=0.3| 


J/ 






ISGW2(c=0.6| 


\\\ 






Gl 


'1 

ill, , , 









0.0 



0.5 



1.0 
r[fm] 



1.5 



2.0 



0.0 




H.O. 

Linear=ISGW2(K=0.3) 
GI=ISGW2(k=0.6) 



1.0 
r[fm] 



2.0 



Figure 4.1: (a) The central potential Vo(r) versus r. Our Coulomb plus HO (solid 
line) and linear (dotted line) potentials are compared with the quasi-relativistic po- 
tential of ISGW2 model with n = 0.3 (long-dashed line) and k = 0.6 (dot-dashed line) 
and the relativized potential of GI model (short-dashed line), (b) The central force 
/o(r) versus r. Our force for the linear potential is the same as that of ISGW2 |^ 
with K = 0.3 (dotted lines). The forces of GI H and ISGW2 with k = 0.6 are the 



same as each other (dashed lines). Our force for the HO potential (solid line) is quite 
comparable with the other four forces up to the range of r < 2 fm. 
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Table 4.1: Optimized quark masses [mq^ms) (in unit of GeV) and the Gaussian 
parameters (3 (in unit of GeV) for both harmonic oscillator and linear potentials 
obtained from the variational principle. The values in parentheses are results from 



Potential 






(3qq 






HO 


0.25 


0.48 


0.3194 


0.3681 (0.3703) 


0.3419 (0.3428) 


Linear 


0.22 


0.45 


0.3659 


0.4128 (0.4132) 


0.3886 (0.3887) 



potential parameters (a, 6, n) regardless of their quark-antiquark contents [see Figs. 
4.2(a) and 4.2(b)]. 

We also determine the mixing angles from the mass spectra of {uj ■,(!)) and {r],r]'). 
Identifying (/i, f2) = {4>, uj) and (r/, r/') for vector and pseudoscalar nonets, the physical 
meson states /i and /2 are given by 



l/i) = — sinS\nn) — cos(5|ss), 
1/2) = cos 6\nn) — sin 6\ss), 



(4.5) 



where \nn) = l/\^\uu + dd) and 5 = 9su{3) ~ 35.26° is the mixing angle. Taking 
into account SU(3) symmetry breaking and using the parametrization for the (mass)^ 
matrix suggested by Scadron |76[, we obtain 



tan b 



^,_{Ml-Ml^{Ml-MX, 



(4.6) 



which is the mo del- independent equation for any meson qq nonets. The details of 
obtaining meson mixing angles using quark-annihilation diagrams are summarized in 
Appendix E. In order to predict the and 17-77' mixing angles, we use the physical 
masses |^ of Mj^ = {ra^.m^ and M^^ = (m(^,m^/) as well as the masses of M]^= 
996 (952) MeV and Mj^-= 732 (734) MeV obtained from the expectation value of H^s 
in Eq. ( [4.1|) for the HO (linear) potential case (see Appendix D for more details). 
Our predictions for uj-(j) and rj-i]' mixing angles for the HO (linear) potential are 
\Sv\ ~ 4.2°(7.8°) and Osuis) ~ -19.3°(-19.6.°), respectively. The mass spectra of 
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Figure 4.2: (a) The parameters niq, nis, Pqq, Pqs, and Pss satisfying the variational 
principle given by Eq. (D.2). The solid, dotted, and dot-dashed lines are fixed by the 
sets of {mu,Puu), {'ITT' s, Pus), and {ms,(3ss), respectively, (b) The parameters m, 



13' 



(3qq, (3qs, and Pss satisfying the variational principle given by Eq. (D.3). The same 
line codes are used as in (a). 
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Table 4.2: Fit of the ground state meson masses with the parameters given in Table 
4.1. Underlined masses are input data. The masses of (77,77') and (a;, 0) were used to 
determine the mixing angles of 77-77' and uj-cf), respectively. The values in parentheses 
are results from the smearing function in Eq. (|D.(j|) instead of the contact term. 



'So 


Expt.[MeV] 


HO 


Linear 




Expt. 


HO 


Linear 


TT 


135 


135 


135 


P 


770 


770 


770 


K 


494 


470 (469) 


478 (478) 


K* 


892 


875 (875) 


850 (850) 


V 


547 


547 


547 


UJ 


782 


782 


782 


v' 


958 


958 


958 





1020 


1020 


1020 



light pseudoscalar and vector mesons used are summarized in Table 4.2. Since the 
signs of 5v for cj-^ mixing are not yet definite [0, [76], ^ in the analysis 

of the quark-annihilation diagram (see Appendix E), we will keep both signs of 5v 
when we compare various physical observables in the next section. 



4.2 Application 

In this section, we now use the optimum model parameters presented in the previ- 
ous section and calculate various physical observables: (1) decay constants of light 
pseudoscalar and vector mesons, (2) charge radii of pion and kaon, (3) form factors 
of neutral and charged kaons, and (4) radiative decay widths for the V{P) P(y)'^ 
and P — > 77 transitions. These observables are calculated for the two potentials (HO 
and hnear) to gauge the sensitivity of our results. 

Our calculation is carried out using the standard LP frame ( = 0) with ql = = 
—q^. We think that this is a distinct advantage in the LFQM because various form 
factor formulations are well established in the LP quantization method using this well- 
known Drell-Yan-West frame {q^ = 0). The charge form factor of the pseudoscalar 
meson can be expressed for the "+" component of the current J'^ as follows: 

F{Q^) = egI{Q^,mg,mg) + eqI{Q'^,mg,mg), (4.7) 
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Figure 4.3: The charge form factor for the pion compared with data taken from 
Ref. ||8^. The sohd and dotted hnes correspond to the results of HO and hnear 
potential cases, respectively. 



where eq{eg) is the charge of quark (antiquark) and 

f-i 



I{Q ,mq,mq 



X 



^ dx J d'k^^^^^<P{x,k^W{x,k', 
A^ + k^- k'x 



+ kl^A^ + k'i 
with the definition of A and k'_|_ given by 

A = xniq + (1 — x)mq, k' ± = k^ + (1 — x)q±. 



(4. 



(4.9) 



The charge radius of the meson can be calculated by r'^=—6dF{Q^)/dQ'^\Q2=Q. Since 
all other formulas for the physical observables such as pseudoscalar and vector meson 
decay constants, decay rates for the V{P) P(y)'y and P ^ 77 transitions, have 
already been given in Chapter 3, we do not list them here again. 

In Fig. 4.3, we show our numerical results of the pion form factor for the HO (solid 
line) and linear (dotted line) cases and compare with the available experimental 
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Table 4.3: Decay constants (in unit of MeV) and charge radii (in unit of fm^) for 
various pseudoscalar and vector mesons. For comparison, we use \6v\ = 3.3° ± 1° for 
both potential cases. The experimental data are taken from Ref . , unless otherwise 
noted. Note that fx* was obtained by {0\sYu\K*) = 2^/^MK*fK'e'"- 



Observables 


Sy = -3.3° ± 1° 


Sy = +3.3° ± 1° 


Experiment 


HO 


Linear 


HO 


Linear 


u 


92.4 


91.8 


92.4 


91.8 


92.4±0.25 


Ik 


109.3 


114.1 


109.3 


114.1 


113.4± 1.1 


fp 


151.9 


173.9 


151.9 


173.9 


152.8± 3.6 


Ik* 


157.6 


180.8 


157.6 


180.8 






45.9±1.4 


52.6±1.6 


55.1±1.3 


63.1±1.5 


45.9± 0.7 


u 


82.6^ 0.8 


94.3^0.9 


76.7+ 1.0 


87.6+1.1 


79.1± 1.3 




0.449 


0.425 


0.449 


0.425 


0.432 ± 0.016 [B8] 




0.384 


0.354 


0.384 


0.354 


0.34± 0.05 [H 


/vi2 


-0.091 


-0.082 


-0.091 


-0.082 


-0.054± 0.101 m] 



data |8^ up to the ~ 8 GeV^ region. Since our model parameters of mu= 0.25 
GeV and Puu= 0.3194 GeV for the HO case are the same as the ones used in Refs. [jlO[ 
and |19[ , our numerical result of the pion form factor is identical with the Fig. 2 (solid 
line) in Ref. [0. In Figs. 4.4(a) and 4.4(b), we show our numerical results for the 
form factors of the charged and neutral kaons and compare with the results of vector 
model dominance (VMD) [^, where a simple two-pole model of the kaon form fac- 
tors was assumed, i.e., Fk+{ko){Q'^) = eu{d)ml/{ml + Q^) + esm^^/im^^ + Q"^). From 
Figs. 4.4(a) and 4.4(b), we can see that the neutral kaon form factors using the model 
parameters obtained from HO and linear potentials are not much different from each 
other even though the charged ones are somewhat different. 



The decay constants and charge radii of various pseudoscalar and vector mesons 
for the two potential cases are given in Table 4.3 and compared with experimental 
data|^, |8^. While our optimal prediction of 5v was l^v^l = 4.2°(7.8°) for HO (linear) 
potential model, we displayed our results for the common 5v value with a small 
variation (i.e., \5v\ = 3.3° ± 1°) in Table 4.3 to show the sensitivity. The results for 
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Figure 4.4: (a) Theoretical predictions of charged form factors using the parame- 
ters of both HO (sohd) and hnear (dotted) potentials compared with a simple two-pole 
VMD model [0] (dot-dashed), F^^f^o) = e„(rf)m2 /(mj + Q^) + e-sm\l{m\ + Q^). (b) 
Theoretical predictions of neutral form factors. The same line codes are used as 
in (a). 
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Table 4.4: Radiative decay widths for the V{P) P{V)'~f transitions. The mixing 
angles, Osu(3) = —19° for rj-r]' and \6v\ = 3.3° ± 1° for u-cj), are used for both potential 



Widths 


= -3.3° ± 1° 


Sv = +3.3° ± 1° 


Expt.[keV] 


HO 


Linear 


HO 


Linear 


1 [p — > 7r^7) 


76 


69 


76 


69 


68 ± 8 


r{uj ^ TT'-f) 


730±1.3 


667±1.3 


730+1.3 


667+1.3 


717 + 51 


r(0 ^ 7r7) 






5 fi+y-y 

3-0-2.9 


e: 1+3.6 

'-'•-'--2.6 


5.8 + 0.6 


r(p ^ ^7) 


59 


54 


59 


54 


58 + 10 


r(^ 777) 


8.7^ 0.3 


7.9^0.3 


6.9+ 0.3 


6.3+0.3 


7.0 + 1.8 


r(0 r/7) 


38. 7± 1.6 


37.8± 1.5 


49. 2± 1.6 


47.6± 1.5 


55.8 + 3.3 


r(r/' ^ p7) 


68 


62 


68 


62 


61 +8 


T{rj' u-f) 


4.9± 0.4 


4.5± 0.4 


7.6± 0.4 


7.0± 0.4 


6.1 + 1.1 


r(0 ^ r]'-f) 


0.41^0.01 


0.39^0.01 


0.36+ 0.01 


0.34+0.01 


< 1.8 


T{K*° 7^07) 


124.5 


116.6 


124.5 


116.6 


117+ 10 


r(ir*+ i^+7) 


79.5 


71.4 


79.5 


71.4 


50 + 5 



both potentials are not much different from each other and both results are quite 
comparable with the experimental data. The decay widths of the V{P) —>■ P(y)7 
transitions are also given for the two different potential models in Table 4.4. Although 
it is not easy to see which sign of 6v for the HO potential model is more favorable to 
the experimental data, the positive sign of 6v looks a little better than the negative 
one for the processes of uj{(j)) rj'-f and 77' w-f transitions. Especially, the overall 
predictions of the HO potential model with positive 6v seem to be in good agreement 
with the experimental data. However, more observables should be compared with the 
data in order to give a more definite answer for this sign issue of uj-(f) mixing angle. 
The overall predictions of the linear potential model are also comparable with the 
experimental data even though the large variation of the mixing angle 6v should be 
taken into account in this case. 



In Table 4.4, we show the results of P{= 'iT,ri,i]') 77 decay widths obtained from 
our two potential models with the axial anomaly plus partial conservation of the axial 
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Table 4.5: Radiative decay widths r(P —>■ 77) obtained by using the axial anomaly 
plus PCAC relations (see Eq. (3.16)). 6*5^(3) = —19° for 77-77' mixing is used for both 
potential cases. The experimental data are taken from Ref. 



Widths 


HO 


Linear 


Experiment 


r(vr 77) 


7.73 


7.84 


7.8 ± 0.5 [eV] 


Y{r] 77) 


0.42 


0.42 


0.47 ± 0.05 [keV] 


Y{r]' 77) 


4.1 


3.9 


4.3 ± 0.6 [keV] 



vector current (PCAC) relations given by Eq. (|3.16|) . The predictions of T]{ri') 77 
decay widths using PCAC are in a good agreement with the experimental data for 
both the HO and linear potential models with rj-rj' mixing angle, 6su(3) = —19°. 
The predictions of the decay constants for the octet and singlet mesons, i.e., rjg and 
T]o, are fs/U = 1.254(1.324) and fo/U = 1.127(1.162) MeV for the HO (linear) 
potential model, respectively. Our predictions of fs and /o are not much different 
from the predictions of chiral perturbation theory |^ reported as fs/fn = 1.25 and 
/o/Ztt = 1.04 ± 0.04, respectively. Another important mixing-independent quantity 
related to /§ and /o is the i? ratio defined by 

.3 1 / /2 £2^ 



R = 



Y(n -> 77) T{ri' 77) 



rri'i 



m: 



Tin 



_ I J TT _|_ Q J TT 



(4.10) 



Our predictions, R = 2.31 and 2.17 for the HO and linear potential model cases, 
respectively, are quite comparable to the available experimental data |^ , -Rexpt = 
2.5±0.5(stat)±0.5(syst). Also, the Q^-dependent decay rates Tp^iQ"^) are calculated 
from the usual one-loop diagram [|l^, ^ and the results are shown in Figs. 4.5-4.7. 
Our results for both potential models are not only very similar to each other but 



53| up to 



also in remarkably good agreement with the experimental data |51|, [ 
the ~ 10 GeV^ region. We think that the reason why our model is so successful 
for P — * 7*7 transition form factors is because the dependence (~ 1/Q^) is due 
to the off-shell quark propagator in the one-loop diagram and there is no angular 



condition [B^ associated with the pseudoscalar meson. 
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Figure 4.5: The decay rate for the vr 7*7 transition obtained from the one-loop 
diagram. Data are taken from Refs. 




2 4 6 8 10 

Q'[GeV'] 



Figure 4.6: The decay rate for the 77 7*7 transition obtained from the one-loop 
diagram. Data are taken from Refs. |5^ . 
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Figure 4.7: The decay rate for the Vj 7*7 transition obtained from the one-loop 



diagram. Data are taken from Refs. [51, 52, 53 



4.3 Summary and Discussion 



In the LFQM approach, we have investigated the mass spectra, mixing angles, and 
other physical observables of light pseudoscalar and vector mesons using QCD-motivated 
potentials given by Eq. ( |4.2| ). The variational principle for the effective Hamiltonian is 
crucial to find the optimum values of our model parameters. As shown in Figs. 4.1(a) 
and 4.1(b), we noticed that both central potentials in Eq. ( ^.2| ) are not only very 
similar to each other but also quite close to the Isgur-Scora-Grinstein-Wise model 
2 (ISGW2) |81[ potentials. In Figs. 4.1(a) and 4.1(b), we have also included the 
Godfrey-Isgur (GI) potential for comparison. Using the physical masses of (cj, 0) and 
(77, r/'), we were able to predict the and 77-77' mixing angles as ~ 4.2°(7.8°) and 
^51/(3) ~ — 19.3°(— 19.6°) for the HO (linear) potential model, respectively. We also 
have checked that the sensitivity of the mass spectra of (t^, 0) to ~ 1°(5°) variation 
of i-e., from by = 4.2°(7.8°) to 3.3° for the HO (linear) potential case, is within 
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the 1%(5%) level. 

Then, we applied our models to compute the observables such as charge radii, decay 
constants, and radiative decays of P{V) V{P)'~f* and P 11* ■ As summarized 
in Tables 4.3, 4.4, and 4.5, our numerical results for these observables in the two 
cases (HO and linear) are overall not much different from each other and are in 



a rather good agreement with the available experimental data |^9|. Furthermore, 
our results of the R ratio presented in Eq. ( [4.101 ) are in a good agreement with the 
experimental data |]5D|, |5T|. The Q"^ dependences of P — 77* processes were also 
compared with the experimental data up to ~ 8 GeV^. The Q"^ dependence 
for these processes is basically given by the off-shell quark propagator in the one- 
loop diagrams. As shown in Figs. 4.5-4.7, our results are in an excellent agreement 



with the experimental data [|5T1 , M i |53| . Both the pion and kaon form factors were 
also predicted in Figs. 4.3 and 4.4, respectively. We believe that the success of 
LFQM hinges upon the suppression of complicated zero-mode contributions from the 
LF vacuum due to the rather large constituent quark masses. The well-established 
formulation of form factors in the Drell-Yan-West frame also plays an important role 
for our model to provide comparable result with the experimental data. Because 
of these successful applications of our variational effective Hamiltonian method, the 
extension to the heavy {h and c quark sector) pseudoscalar and vector mesons will be 
shown in Chapter 7 and the 0^^ scalar mesons is currently under consideration. 

While there have been previous LFQM results on the observables that we calculated 
in this work, they were based on the approach of modeling the wave function rather 
than modeling the potential. Our approach in this work attempting to fill the gap 
between the model wave function and the QCD-motivated potential has not yet been 
explored to cover as many observables as we did in this work. Nevertheless, it is not 
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yet clear which sign of uj-(j) mixing angle should be taken, even though the overall 
agreement between our HO potential model with the positive sign, i.e., 5y ~ 3.3°, 
and the available experimental data seem to be quite good. If we were to choose the 
sign of X as X > in Eq. (E.4), then the fact that the mass difference — nip 
is positive is correlated with the sign of the uo-cf) mixing angle In other words, 

> implies 5y > from Eqs. (E.3)-(E.5). Perhaps, the precision measurement 
of — > ?7'7 envisioned in the future at TJNAF experiments might be helpful to give 
a more stringent test of 5y. In any case, more observables should be compared with 
the experimental data to give more definite assessment of this sign issue. 
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Chapter 5 

Exploring the Timelike Region For 
the Elastic Form Factor in a Scalar 
Field Theory 



The Drell-Yan-West (g'^=g°+g^=0) frame in the LF quantization provided an effective 
formulation for the calculation of various form factors in the spacelike momentum 
transfer region q^=—Q^< [40|. In g+=0 frame, only parton-number-conserving Fock 



state (valence) contribution is needed when the "good" components of the current, 
and 3 _L={Jx) Jy), are used For example, only the valence diagram shown in 



Fig. 1.1(a) is used in the LFQM analysis of spacelike meson form factors. Successful 



LFQM description of various hadron form factors can be found in the literatures ||T0 

However, the timelike (g^> 0) form factor analysis in the LFQM has been hindered by 
the fact that g+=0 frame is defined only in the spacelike region {q'^=q^q^ ~ q']_< 0). 
While the q^^^ frame can be used in principle to compute the timelike form factors, 
it is inevitable (if q^^^) to encounter the nonvalence diagram arising from the quark- 
antiquark pair creation (so called "Z-graph"). For example, the nonvalence diagram 
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in the case of semileptonic meson decays is shown in Fig. 1.1(b). The main source 
of the difficulty, however, in calculating the nonvalence diagram (see Fig. 1.1(b)) 
is the lack of information on the black blob which should contrast with the white 
blob representing the usual LF valence wave function. In fact, it was reported |^ 



that the omission of nonvalence contribution leads to a large deviation from the full 
results. The timelike form factors associated with the hadron pair productions in 
e^e~ annihilations also involve the nonvalence contributions. Therefore, it would be 
very useful to avoid encountering the nonvalence diagram and still be able to generate 
the results of timelike form factors. 

In this work, we show an explicit example of generating the exact result of the timelike 
form factor without encountering the nonvalence diagram. This can be done by the 
analytic continuation from the spacelike form factor calculated in the Drell-Yan-West 
(g^=0) frame to the timelike region. To explicitly show it, we use an exactly solvable 
model of (3 + 1) dimensional scalar field theory interacting with gauge fields. Our 
model is essentially the (3 + 1) dimensional extension of Mankiewicz and Sawicki's 
(1 + 1) dimensional quantum field theory model |^3[, which was later reinvestigated 



by several others ||2^, |2^, 0, ^ ^ . The starting model wave function is the solution 



of covariant Bethe-Salpeter (BS) equation in the ladder approximation with a rela- 



tivistic version of the contact interaction |93|. Here, we do not take the Hamiltonian 



approach. The covariant model wave function is a product of two free single particle 
propagators, the overall momentum-conservation Dirac delta, and a constant vertex 
function. Consequently, all our form factor calculations are nothing but various ways 
of evaluating the Feynman perturbation-theory triangle diagram in scalar field theory. 
In this model, we calculate: (A) the timelike process of 7* M + M transition in 
g+T^O (g^>0) frame, (B) the spacelike process of M — > 7* + M in g^T^O (g^<0) frame, 
and (C) the spacelike process of M — ^ 7* + M in g+=0 frame. Using the analytic 
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continuation from g^<0 to q^>0, we show that the result in (C), which is obtained 
without encountering the nonvalence contributions at all, exactly reproduces the re- 
sult in (A). In fact, all three results (A), (B), and (C) coincide with each other in the 
entire range. We also confirm that our results are consistent with the dispersion 
relations ^ ^ |100|] . We consider not only for the equal quark/antiquark mass 



case such as the pion but also for the unequal mass cases such as K and D. 

This Chapter is organized as follows: In Section 5.1, we derive the timelike electro- 
magnetic (EM) form factor of 7* M+M process in the frame (A) and discuss 
the singularities occuring from the on-energy shell of quark-antiquark pair creation. 
In Section 5.2, the spacelike form factor of M — > 7* + M process is calculated both 
in the q~^^0 (B) and q~^=0 (C) frames. We then analytically continue the spacelike 
form factors to the timelike region. The singularities occured in the timelike region 
are also discussed. In Section 5.3, for the numerical calculation of the EM (vr, K, 
and D) meson form factors for three different cases (A), (B) and (C), we use the con- 
stituent quark and antiquark masses {mu=md=0.25 GeV, ms=0.48 GeV, and mc=1.8 
GeV) [|13l, Q and show that the form factors obtained from those three different 
cases are indeed equal to each other for the entire region. The meson peaks anal- 
ogous to the vector meson dominance (VMD) are also obtained. The conclusion and 
discussion follows in Section 5.4. 



5.1 Form Factors in the Timelike Region 

The EM local current J^(0) responsible for a virtual photon decay into two qQ bound 
states in the scalar field theory can be calculated using the diagrams shown in Fig. 
5.1. The covariant diagram shown in Fig. 5.1(a) is equivalent to the sum of two LF 
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b b c c 




Figure 5.1: The electromagnetic decays of a photon into two-body bound sates, i.e., 
7* — >• gg(or QQ)^ Jli{qQ)Jli{qQ), in scalar theory: Covariant representation (a), 
and the LF time ordered contributions to the decay amplitude (b) and (c). 

time-ordered diagrams in Figs. 5.1(b) and 5.1(c). The EM current J^{0) obtained 
from the covariant diagram of Fig. 5.1(a) is given by 

J^{0) = ejd'k- {q-2kr- — — 1 ^— - 

J [q ~ ky — + le {Q — k — P2) — + le 

^ 72 2 ~ ~'~ ^q("^9 ^ "^<9 ^^^^ term), (5.1) 

k 

where ?Tig(Q) and eg(Q) are the constituent quark (antiquark) mass and charge, respec- 
tively. The corresponding form factor F{q^) of the qQ bound state in timelike(g^ > 0) 
region is defined by 

J^(0) = (Pi-P2)'^i^(g'), (5.2) 

where q = + P2, = P| = and M is the mass of a qQ bound state scalar 
particle. Using the Cauchy integration over k~ in Eq. ( |5.1D , we can find each time- 
ordered contribution (Figs. 5.1(b) and 5.1(c)) to the timelike form factor -F(g^ > 
4M^) in Eq. ( [5.2[ ). This procedure allows us to analyze the singularity structure of 
each LF time-ordered diagram as well. For the calculation of each LF time-ordered 
contribution, we take the purely longitudinal momentum frame, i.e., q^^O, q±=0 and 
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Pi^=P2±=0. Accordingly, g^(= q^q ) > 4M^ is given by 

g2 = M2(l + a)Va, (5.3) 

where a=P'^ / P^=q~^ / — 1 is the longitudinal momentum fraction and the two 
solutions for a are given by 

Note that both aj-=l correspond to the threshold, g^=4M^. The EM form factor 
F(g^) in Eq. (|5^ ) is independent of the subscript sign of a. Thus, one can take either 
or a_ to calculate F{q^). Here, for convenience, we use a=a_ which ranges from 
to 1 for the physical momentum transfer region, i.e., a_ — * as ^ oo and 1 as 
q^ — > 4M^. Of course, we can use a=a+ equally well and verify that the two results 
{a+ and aS) are exactly same for the calculation of F{q^). 

Since > Pi > P2 > for a=a_, the Cauchy integration over k~ in Eq. ( |5.1|) has 
two nonzero contributions to the residue calculations, one coming from the interval 
(i) < < Pi (see Fig. 5.1(b)) and the other from (ii) Pi < < q^ (see Fig. 
5.1(c)). The internal momentum k~^ is defined by k^=xPi', where x is the Lorentz 
invariant longitudinal momentum variable. The "good" -current J~^{0) is used in our 
computation of the two LF diagrams Figs. 5.1(b) and 5.1(c). In the following, for 
simplicity, we won't explicitly write either the obvious second term in Eq. ( |5.1|) nor 
the charge factor (e, or eq). 

In the region of < A;+ < P^^, the residue is at the pole of k~ = ['m'^ + A;^ — ie]/k~^, 
which is placed in the lower half of complex-^" plane. Thus, the Cauchy integration 
of in Eq. (|5.1| ) over k~ in this region gives 

/■P+ fy+ _ 2k^ 
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X 



X 



{ml + kl)/k+ - {ml + fc2 )/(gH 



(5.5) 



q- - - {ml + kl)/k+ - (m| + kl)/{P+ - k+) 

where the subscript h in J{j''(0) imphes the current in Fig. 5.1(b) and M{=—m) is the 
normahzation constant. The last two factors in Eq. (|5.5|) correspond to the two and 
three particle intermediate states. We represent the energy denominators of the two 
and three particle intermediate states as E\ and in Fig. 5.1(b), respectively. 



To analyze the singularities of Eq. ( p.5| ), we further integrate over and obtain 

'£}> 

dx— j:^ ^ln( 



t4-/^n ^r 7 x(l + a - 2x)/(l + a) , ,,^2 



(5.6) 



where £^2=^(1 + ct — x)M'^/a — m^ and £^|=x(l — x)M^ — [xm^ + (1 



x]m 



While 



is not zero {S^ 7^ 0) in general for the entire physical region, can be zero when 
> Am'^^Qy The singular structure of S!^ — £2 term in Eq. ( [5.6| ) depends on whether 
a qQ bound state scalar particle is strongly bounded (M^ < m^ + m^) or weakly 
bounded (M^ > + wi^). As we will show in our numerical calculations (Section 
5.4), anomalous threshold appears for > m^ + m^ while only the normal threshold 
of bound state exists for < + m^. 

In the region of Pi < k^ < q^, the residue is at the pole of k^=q^ — [nig + (qj_ — 
kx)^ — «e]/(g+ — /c"*"), which is placed in the upper half of complex-/c~ plane. Thus, 
the Cauchy integration of J^{0) in Eq. ( |5.1| ) over k~ in this region yields the result 



+ 



2k^ 



-U dk+d^k^—^ ^ , ,^ 

Jp+ k+ {q+ - k+){P^ - A;^ 

1 



X 



{ml + kl)/k+ - {ml + kl)l{q+ - k+) 



(5.7) 



q- - P{ + (m| + kl)/{P+ - k+) - {ml + kl)/{q+ - k^ 
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where the subscript c in J^{0) means the current in Fig. 5.1(c). After the integration 
over the in Eq. ( |5.7|) , we obtain 

1 ^^aX{l + a- 2aX)/{l + a) ,^ fS'^ 



where x=l + a{l - X), £I=X{1 - X)M'^ - [Xm| + (1 - X)m2] and £^=aX{l + a- 
aX)M'^ /a — rn?. The pole structure in Eq. ( |5.8|) is equivalent to that of Eq. ( ^.6|) . 



Consequently, the timelike form factor in Eq. ( |5.2|) is given by 



ax{l + a — 2ax)^ / £^2 



.3 02 \<-3 

where £2, £3, £2 £3 are defined in Eqs. (|5.6|) and (^]8| 



5.2 Form Factors in Spacelike Region and the An- 
alytic Continuation to the Timelike Region 

In this section, we calculate the EM form factor in spacelike momentum transfer 
region and then analytically continue to the timelike region to compare the result 
with the timelike form factor (i.e. Eq. ( p.9| )) that we obtained in the previous section. 
The EM current of a qQ bound state in spacelike momentum transfer region is defined 
by the local current j^(0); 

r{o) = {p, + P2rHq'), (5.10) 

where q=Pi — P2, q^<0 and J-'{q^) is the spacelike form factor. The EM current j^(0) 
obtained from the covariant triangle diagram of Fig. 5.2(a) is given by 

(5.11) 
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(b) 



(c) 



Figure 5.2: Covariant triangle diagram (a) is represented as the sum of LF triangle 
diagram (b) and the LF pair-creation diagram (c). 

As in the case of the timelike form factor in Section 5.2, the Cauchy integration of 
k~ in Eq. ( |5.11D has also two contributions to the residue calculations, one coming 
from the interval < < P2 (see Fig. 5.2(b)) and the other from P2 < < 
(see Fig. 5.2(c)). Again, only the "good" -current j'^{0) in Eq. (|5.11|) is used to obtain 
the contributions from Figs. 5.2(b) and 5.2(c). 



In the region of < /c+ < P2, the residue is at the pole of k~ = [mQ + k\ — ie]/k'^ , 
which is placed in the lower half of complex-/c~ plane. Thus, the Cauchy integration 
of in Eq. (|5.11|) over k~ in this region yields 



(Pi + P2 - 2A;)- 



■^io '^^^'^'^^lP,-ky{P,-k-qyk^ 

1 



X 



X 



Pf - [ml + kl]/{P+ - k+) - [ml + kl]/k^ 

1 



(5.12) 



Pf - g- - [ml + (kx + q±)2]/(Pi+ - k+) - [m| + kl]/k^ 
where the subscript h in implies the current in Fig. 5.2(b). 

In the region of P2 < k^ < P^, the residue is at the pole of k^=P{ — [m^ + k\ — 
it]/ {Pi — k^), which is placed in the upper half of complex-A;~ plane. Thus, the 
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Cauchy integration of in Eq. ( p.ll| ) over k in this region becomes 
'■P^ - (Pi + P2- 2k)+ 



M / dk+<fki 



(Pi - A;)+(Pi -k- q)+k^ 
1 



X 



Pf - K + kl]/{P^ - k+) - [ml + kl]/k+ 



(5.13) 



p- _ p- + [ml + A;i]/(P+ - k+) - [ml + (q^ + kx)2]/(P2+ - A;+ 

As one can see from Eq. ( |5.13| ), the nonvalence contribution (Fig. 5.2(c)) vanishes 
only in q'^=0 frame. In the following, we investigate the spacelike form factor J^{q^) 
given in Eq. ( p.l(j| ) using both g^T^O and g+=0 frames. We then analytically continue 
to the timelike region in order to compare the result with the direct calculation of 
the timelike form factor F{q^) presented in the previous section. 



5.2.1 The 7^ and q_L = frame 



In the purely longitudinal momentum frame g^T^O, q_L=0, and Pi_l=P2±=0, the 
momentum transfer q^=q^q^ can be written in terms of the longitudinal momentum 
fraction p=P+/P+ = 1 - q+/P^- 

q"^ = -M\l- (3)^(3 <0, (5.14) 

where the two solutions of (3 are given by 




(5.15) 



The form factor J-'{q^) in Eq. (|5.10| ) is also independent of the subscript sign of /?. 
However, < P^ was used in obtaining Eqs. ( [5. 121 ) and ( |5.13D and thus here we 
use P=P- (0 < /9 < 1) in spacelike region. As shown in Eqs. ( |5.11| )-( pT^ ), the sum 
of valence (Fig. 5.2(b)) and nonvalence (Fig. 5.2(c)) diagrams is equivalent to the 
covariant triangle diagram in Fig. 5.2(a). 
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For the analysis of singularity structures, we integrate over k_L and obtain from the 
valence contribution(Fig. 5.2(b)); 

jm = ^^l ci.^^-^ln(|), (5.16) 

where £|=x(l— x)M^ — [xmg + (l— x)m|] and S2=Px{l—(3x)M'^ — [(3xm'^+{l—(3x)mQ]. 
It turns out that Eq. (|5.16|) has no singularities because S!^j^O, ^^It^O, and — ^^It^O 



for the entire region. On the other hand, the k±_ integration for the current j^(0) 
in Eq. ( |5.13| ) yields 

,?(0)^.Arr^^ '^-t-^'f-^ ..f|Y (5,17) 



"0 S2 ^£2' 

where x=l - (1 - p)X, £I={1 - P)X[l - (1 - I3)X]M'^ - [(1 - /9)A'm| + (1 - (1 - 
P)X)ml] and 4'=-(l - - ^)MV/3 - m^. While 4' corresponding to the 

energy denominator of the two particle intermediate state does not vanish both in 
spacelike and timelike region, of the three-particle intermediate state is not zero 
only for the spacelike momentum transfer region. For the timelike region, can 
be zero so that the singularities start at g^^jjj=4m^(-Q^ for Yl^iYQQ) vertex. The 
singularity structure of S2 — in Eq. ( |5.17| ) is the same as in the case of timelike 



form factor (Section 5.4), following the condition of a qQ bound state. 

The EM form factor ^{q"^) in Eq. ( |5.10| ) of a qQ bound state in q^^O frame is then 
obtained by 



where £^25^35^2 £2 are defined in Eqs. (|5.16| ) and (|5.17|) . Here, /3 is a function 



of q^. According to the analytic continuation, the sign of q^ in Eq. (|5l8D must be 
changed from — to + for the timelike region. 
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5.2.2 The Drell-Yan-West {q+ = 0) frame 



In Drell-Yan-West frame, = 0, = —q'j_, and Pi_l=P2±=0, the '+' component of 
the current has only the valence contribution, i.e., jfc'(O) in Eq. ( |5.12|) . The current 
j^(0) in g+=0 frame is given by 



rl 2x(l — x) 



(5.19) 



where i±=k± + xq±/2, A=x{l - x)M^ - [xm^ + (1 - a;)m|], and ^^=x^g^/4. The 
angle (0 < < 27r) is defined by i± ■ q±= |^±||q±| cos0. 



Integrating Eq. (|5.19|) over and we obtain 



dx- 



:i-x) 



^gi(4M2 + ql) y^(x - x+)(x - x_) 

X tanh"^ 



4M2 + ql J[x-x+){x-x.) 



where 



2(M2 -ml + nA) 
x± = ^ — — ± 



4(M2 -m2 + m|)2 



(5.20) 



(5.21) 



4M2 + g2 (4M2 + g2)2 (4M2 + g2)- 

The analytic continuation from spacelike to timelike region in the g"^=0 frame requires 
the change of gj_ to iq± in Eqs. (|5.20| ) and ( |5.21|) . We note from Eqs. ( |5.20| ) and ( ^.21| ) 
that the result of the timelike region exhibits the same singularity structure as the 
direct analyses in q^^O frame, i.e., Eqs. (|5.6|) and ( p.8|) , even though the nonvalence 
contribution in Fig. 5.2(c) is absent here. 

After some manipulation, we obtain the EM form factor of a qQ bound state in the 
q~^=0 frame as follows 



^(g^g+ = 0) 



arcsin{i^) 



X tanh ^ 



g2 a + b sin 6 



d6{l — a — 6 sin I 



g2 - 4M2 ib cos 6 



(5.22) 
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where a = (x+ + x_)/2, b = (x+ — x_)/2, and q'^=—q'j_- While the representations in 
Eqs. ( 5.18|) and ( ^.221) look apparently different, it is amazing to realize that the two 
formulas (Eqs. ( |5.18| ) and ( |5.22| )) turn out to be actually identical. As we will show 
explicitly in the next section of numerical calculations, all three results of Eqs. (|5.9|) , 
( ^.18| ) and (|5.22|) indeed coincide exactly in the entire range. 



5.3 Numerical Results 



For our numerical analysis of vr, K, and D meson form factors, we use the physical 
meson masses together with the following constituent quark and antiquark masses: 
mu=md=0.25 GeV, ms=0.48 GeV, and mc=1.8 GeV ||13|, ^ Since our numerical 
results of the EM form factors obtained from Eqs. ( p. 91 ), ( |5.18| ) and ( |5.22D turn out to 
be exactly same with each other for the entire region, only a single line is depicted 
in Figs. 5.3, 5.4 and 5.5 for the form factor calculations of vr, K, and D mesons, 
respectively. 

In should be noted from our constituent masses that < + rriQ for vr and K 
and > + tuq for D meson cases. As discussed in Ref. for the analysis of 
the one-particle matrix element of a scalar current, the sigularity for > + itlq 
case starts at 

9min = -^K(Q) - - ^QmMM + - ml^Q^], (5.23) 

for 'y*qq{'y*QQ) vertex, while the singularity for < + itlq case starts on the 
positive g^-axis at the threshold point g^;j^=4m^^g^ for 'y*qq{l*QQ) vertex. Our nu- 
merical results exhibit all of these threshold behaviors coming from the normal (vr, K) 
and anomalous (D) cases. As a consistency check, we also compare our numerical 
results of the form factor F{q^) = Re F{q^) +i Im F{q^) with the dispersion relations 
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given by 



ReF(,^) ^ 1P^^^^.'^ (5.24) 



= 

TT J-oo q'^ — q^ 



ImF(g^) = -'-Pf^^^^dq-, (5.25) 
TT J— oo q — q 



where P indicates the Cauchy principal value. 

In Fig. 5.3(a), we show the EM form factor of the pion for —2 GeV^ < < 3 GeV^. 
The imaginary part (the dotted line) of the form factor starts at g^jj^=4m^^^=0.25 
GeV, which is consistent with the condition for < + case. It is interesting 
to note that the square of the total form factor |F^(g^)p (thick solid line) produces 
a p meson-type peak near ~ M^. However, it is not yet clear if this model indeed 
reproduces all the features of the vector meson dominance (VMD) phenomena. Even 
though the generated position of peak is consistent with VMD, the final state inter- 
action is not included in this simple model calculation. We believe that much more 
complex mechanisms may be necessary to reproduce the realistic VMD phenomena. 
More detailed analysis along this line is under consideration. Nevertheless, it is re- 
markable that this simple model is capable of generating the peak and the position 
of peak is quite consistent with the VMD. 

In Fig. 5.3(b), we show the timelike form factor of the pion for the entire > region 
and compare the imaginary part of our direct calculations (dotted line) obtained from 
Eqs. ( [5.9| ), ( |5.18| ), and ( |5.22| ) with the result (data of black dots) obtained from the 



dispersion relations given by Eq. (5.25). Our direct calculation is in an excellent 
agreement with the solution of the dispersion relations. Our result for the real part 
are also confirmed to be in complete agreement with the dispersion relations. For 
high q"^ region, the imaginary part of the form factor is dominant over the real part 
(thin solid line). 
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Figure 5.3: (a) The electromagnetic form fator of the pion in (3 + 1) dimensional 
scalar field theory for —2 < < 3 GeV^. The total, real, and imaginary parts of 
|F7r(g^)P are represented by thick solid, solid, and dotted lines, respectively, (b) The 
electromagnetic form fator of the pion in (3 + 1) dimensional scalar field theory for 
the entire timelike region compared to the dispersion relations (data of black dots) 
given by Eq. (5.25). The same fine code as in (a) is used. 
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Figure 5.4: The electromagnetic form fator of the kaon in (3 + 1) dimensional scalar 
field theory for — 2 < 5^ < 5 GeV^. The same hne code as in Fig. 5.3(a) is used, (b) 

The electromagnetic form fator of the kaon in (3 + 1) dimensional scalar field theory 
for the entire timelike region compared to the dispersion relations (data of black dots) 
given by Eq. (5.25). The same line code as in Fig. 5.3(a) is used. 
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Figure 5.5: The electromagnetic form fator of the D meson in (3 + 1) dimensional 
scalar field theory for —10 < < 30 GeV^. The same line code as in Fig. 5.3(a) is 
used. 

In Fig. 5.4(a), we show the kaon form factor for —2 GeV^ < 5'^ < 5 GeV^. The kaon 
also has the normal singularity. However, it has two thresholds for the imaginary 
parts; one is g^in=4m^ and the other is g^in=4mg. These lead to the humped shape 
(dotted line) of the imaginary part shown in Fig. 5.4(a). While we have in principle 
two vector-meson-type peaks (i.e. p and (p), one can see in Fig. 5.4(a) only (p meson- 
type peak for the timelike kaon EM form factor above the physical threshold at 
g^jj^=4M^. We also show in Fig. 5.4(b) the imaginary part from our direct calculation 
is in an excellent agreement with the result (data of black dots) from the dispersion 
relations for the entire timelike region. Again, the imaginary part is predominant 
for high region. 

In Fig. 5.5, we show the D meson form factor for -10 GeV^ < < 30 GeVl Unlike 
the normal threshold of vr and K form factor calculations, the D meson form factor 
shows anomalous thresholds according to Eq. ( ^.23| ), i.e., g^^^ ~ 0.24 GeV^ (compared 
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to 4m^=0.25 GeV^ for normal case) and g^j^ ~ 12.4 GeV^ (compared to 4m^=12.96 
GeV^ for normal case) for the 7* — c? and 7* — c vertices, respectively. Similar to the 
kaon case in Fig. 5.4, we also have two unphysical peaks, i.e., p and J/tpllS) meson 
type peaks due to d and c quarks, respectively. However, the timelike form factor 
of D meson has no pole structure for the the physical > 4M|, region. In all of 
these figures (Figs. 5.3-5.5), it is astonishing that the numerical result of Eq. ( p.22| ) 
obtained from g+=0 frame without encountering the nonvalence diagram coincides 
exactly with the numerical results of Eqs. ( [5. 91 ) and ( |5.18| ) obtained from q^^O frame. 



5.4 Summary and Discussion 



In this chapter, we presented the investigation on the EM form factors of various 
mesons both in spacelike and timelike regions using an exactly solvable model of 
(3 + 1) dimensional scalar field theory interacting with gauge fields. Our calculations 
demonstrated that one can compute the timelike form factor without encountering the 
nonvalence contributions. We calculated the form factor in spacelike region using the 
Drell-Yan-West {q^=0) frame and showed that its analytic continuation to the time- 
like region reproduces exactly the direct result of timelike form factor obtained in the 
longitudinal momentum {q'^y^O and q±=0) frame. It is remarkable that the analytic 
continuation of the result in Drell-Yan-West frame to the timelike region automati- 
cally generate the effect of the nonvalence contributions to the timelike form factor. 
Another interesting result in our model calculations is that the peaks analogous to the 
VMD were generated and the position of peaks were indeed quite consistent with the 
VMD. Even though much more detailed analyses including the final state interaction 
may be necessary to reproduce the entire feature of VMD, our results seem pretty 
encouraging for further investigations. Using the dispersion relations, we have also 
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confirmed that our numerical results of imaginary parts start at ?min— and the 
normal thresholds appear for tt and K (M^ < 171^+ itlq) systems while the anomalous 
threshold exists for D (M^ > m'^ + ttiq) system. Thus, it is hopeful that one can use 
the same technique of analytic continuation and calculate the timelike form factors 
in more realistic models. Detailed anlysis along this line is underway. Applications 
to the semileptonic decay processes in (3 + 1) dimensional scalar field model are also 
in progress. 
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Chapter 6 

Nonvanishing Zero Modes in the 
Light-Front Current 



One of the distinguishing features in LF quantization is the rational energy-momentum 
dispersion relation which gives a sign correlation between the LF energy(P~) and the 
LF longitudinal momentum(P+). In the old-fashioned time-ordered perturbation the- 
ory [KIT], this sign correlation allows one to remove the so-called "Z-graphs" such as 



the diagram of particle-antiparticle pair creation(annihilation) from(to) the vacuum. 
As an example, in the theory of scalar fields interacting with gauge fields [0, the 



covariant triangle diagram shown in Fig. 5.2(a) corresponds to only two LF time- 
ordered diagrams shown in Figs. 5.2(b) and 5.2(c), while in the ordinary time-ordered 
perturbation theory. Fig. 5.2(a) would correspond to the six time-ordered diagrams 
including the "Z-graphs". Furthermore, the Drell-Yan-West (^^=0) frame may even 
allow one to remove the diagram shown in Fig. 5.2(c) because of the same reason- 
ing from the energy-momentum dispersion relation and the conservation of the LF 
longitudinal momenta at the vertex of the gauge field and the two scalar fields. 

Based on this idea, the Drell-Yan-West {q~^=0) frame is frequently used for the bound- 
state form factor calculations. Taking advantage of q^=0 frame, one may need to 
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consider only the valence diagram shown in Fig. 5.2(b), where the three-point scalar 
vertices should be replaced by the LF bound-state wavefunction. 

In this work, however, we point out that even at q'^=0 frame one should not over- 
look the possibility of non-zero contribution from the nonvalence( pair creation or 
annihilation) diagram shown in Fig. 5.2(c). As we will show explicitly in the simple 
(1 -|- 1) -dimensional scalar field theory interacting with gauge fields, the current J" is 
not immune to the zero mode contribution shown in Fig. 5.2(c) at g^=0. While the 
current J"*" does not have any zero mode contribution from Fig. 5.2(c), the processes 
that involve more than one form factor, e.g., semileptonic decay processes, require 
the calculations of more components of the current other than J"*" in order to find 
all the necessary form factors in q~^—0 frame. For instance, in the analysis of the 
semileptonic decays between two pseudoscalar mesons, two form factors, f±{q^), are 
involved and one has to use not only J"*" but also J~ (or J-*- in 3 -|- 1 dimensions) to ob- 
tain both form factors in q~^=0 frame. Thus, the zero mode contribution is crucial to 
obtaining the correct results of electroweak form factors. Only a brief exactly solvable 
model calculation is provided here. A full, detailed treatment of (3 + 1) dimensional 
semileptonic decay processes such as K ^ n, B ^ tt, B ^ D etc. will be presented 
in a separate communication. 

This Chapter is organized as follows: In Section 6.1, we describe the general formalism 
of the semileptonic decay form factors for non-zero momentum transfer in (1-1-1)- 
dimensions and then discuss the zero mode problem in the limiting cases of the form 
factors as g"*" 0. The numerical results of the zero mode contributions for K ^ ir, 
S — > TT, and B ^ D transitions are given in Section 6.2. We also briefly discuss 
in Section 6.2 the zero modes in the electromagnectic form factor calculations. The 
summary and discussion follow in Section 6.3. 
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6.1 Zero modes in ^0 Transitions 



The semileptonic decay of a Qiq bound state into another Q2q bound state is governed 
by the weak current, viz., 

no) = {P2\Q2rQi\Pi) = + + (^1 - ^2)^ (e.i) 

where P2—P1 — q and the non-zero momentum transfer square q^—q~^q~ is timehke, 
i.e., q'^—[0, (Ml — M2)^]. One can easily obtain q'^ in terms of the fraction a as follows 

q'^{l-a){M',-^), (6.2) 
a 

where a—P2/Pi—l — q'^/Pi- Accordingly, the two solutions for a are given by 



M2 
Ml 



Ml + M| - g2 hMl + M| - g2x 2 



±W ^„ ^ " -1 



(6.3) 



2M1M2 V V 2M1M2 

The +(— ) sign in Eq. (6.3) corresponds to the daughter meson recoiling in the pos- 
itive(negative) ^-direction relative to the parent meson. At zero recoil (g^^g^g^) and 
maximum recoil(g^=0), a± are given by 

a^O) = 1, a_(0) = (^j . (6.4) 

In order to obtain the form factors f±{q^) which are independent of a±, we can define 

{P2\Q2rQi\Pi)\a=a^ = 2P+i7+(a±) for /i = +, (6.5) 

= 2i^-^)H-{a^) for = (6.6) 

and obtain from Eq. (6.1) 
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where (5±—a^{Q)/a±. 

Now, the current J'*(0) obtained from the covariant triangle diagram of Fig. 5.2(a) is 
given by 

1 1 



P2 — ky — ml + iek'^ — m? + ie 

(6.9) 



Prom this, we obtain for the "±" -components of the current J^{0) as 

J±(0) = -27ri(/± + /±), (6.10) 

where and I2 corresponding to diagrams Figs. 5.2(b) and 5.2(c), respectively, are 
given by 

^ V X J \ a a—x X 



/■I ^ — It A- n 
I}(a) = / dx ' 7° , ; ; ^, (6.12) 

and 



ra M^ + Ml/a-2myx 
L [a) — I ax 7 ; rvS — 

^» .(l-.)(„-.)(Mf-^-I?)(^-^ 



(6.13) 



.(1 - .)(a - .)(Mf - - ^) (^ + - ^ - Mf) 
Note that at zero momentum transfer limit, q^—q^q~ — > 0, the contributions of 
I^{a) come from either limq+^Ql^(a)=I^{a^{0)) or limq-^Ql^{a)=I^{a^{0)). It 
is crucial to note in q~^—0 frame that while /2'(q;+(0)) vanishes, /2~(q;+(0)) does not 
vanish because the integrand has a singularity even though the region of integration 
shrinks to zero. Its nonvanishing term is thus given by 
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Table 6.1: Form factors of /±(0) obtained for different zero- momentum transfer limit, 
g"'"=0 and q~=0. The notations of a±, ap, and am used in table are defined as a± = 
a±{0), ap=l + a_(0), and 0^=1 — a-(0), respectively. lf{a±) implies ^21=1 lt{oi±)- 



Form factor 


q+ = 


g~ = 


/+(0) 


/r(a+)/2 


/r(a_)/2M2 


/-(O) 


[/rK)/M2-apJ+(a+)/2]/a„ 


[/+(«_) - apIr{a.)/2Mf]/am 



Table 6.2: Zero-mode(Z.M.) and nonvalence(N.V.) contributions to the exact form 
factors of /±(0) for the semileptonic decays of K{B) —>■ vr and i? ^ D in (1 + 1) 
dimensions. We distinguished the zero mode contribution at g+=0 from the usual 
nonvalence one at q~=0. 



Frame 




N.V.(Z.M.) factor 


K 


B ^ir 


B ^ D 


g+ = 


^..M.(0)//r(0) 


None 











j^.M.^o)/r^"(o) 


(X J2-(a+(0)) 


6.9 


0.03 


0.1 


q- = 


/^■Mo)//r(o) 


oc /2"(a-(0)) 


2.8 


1.3 


0.05 


jN.vi.^o)//i""(o)^"J 


oc /2+(a_(0)) 


3.8 


3.8 


0.6 


JlN.V..^O)/f""(0)^"J 


oc/2-(a-(0)) 


-11.1 


-4.0 


-1.1 



H We show the separate contributions of the nonvalence terms proportional to 
/^(a_(0)) and /^(a_(0)) to the exact form factor of /-(O) at q~=0. 

This nonvanishing term is ascribed to the term proportional to k^=Pi —m\/ (P^—k^) 
in Eq. (6.14), which prevents Eq. (6.14) from vanishing in the limit, a — 1. This 
is precisely the contribution from "zero mode" at g"^=0 frame. l2{<y+{0)) should be 
distingushed from the other nonvanishing pair-creation diagrams at q^=0 frame, i.e., 
/^(a-_(0)). Some relevant but different applications of zero modes were discussed in 
the literature ^ . 



In Table 6.1, we summarized the form factors /±(0) obtained from both currents, 
and J~ , for different zero momentum transfer limit, i.e., g^=0 or q~=0. As shown in 
Table 6.1, the nonvalence contributions, l^{a±{0)), are separated from the valence 
contributions, l^{a±{0)). Of special interest, we observed that the form factor /-(O) 
at q^=0 is no longer free from the zero mode, /2"(a;_|_(0)). 
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6.2 Numerical Results 

To give some quantitative idea how much these nonvalence contributions /;^(q;±(0)) 
are for a few different decay processes, we performed model calculations for X — > tt, 
B ^ TT, and B ^ D transitions in (1 + 1) dimensions using rather widely used con- 
stituent quark masses, mu((i)=0.25 GeV, mc=1.8 GeV, and mb=5.2 GeV. Numerically, 
we first verified that the form factors, /+(0) and /-(O), obtained from the frame 
are in fact exactly the same with /+(0) and /-(O) obtained from the q~—Q frame, 
respectively, once the nonvalence contributions (including zero mode) are added. The 
nonvalence contributions to the form factors of /±(0) at q'—Q are also shown in Ta- 
ble 6.2. In Figs. 6.1a(b)-6.3a(b), the effects of pair-creation (nonvalence) diagram to 
the exact form factors are shown for the non-zero momentum transfer region for the 
above three decay processes. Especially, the zero mode contributions /^(q;+(0)) to 
the exact solutions for the /_(0) at g+=0, i.e., /^■^•(0)//!"ii(0), are estimated as 6.9 
for X ^ TT, 0.03 for B ^ and 0.12 ior B ^ D decays. The zero mode contribu- 
tions on /_(0) at q~^=0 frame are drastically reduced from the light-to-light meson 
transition to the heavy-to-light and heavy-to-heavy ones. This qualitative feature of 
zero mode effects on different initial and final states are expected to remain same 
even in (3 + 1) dimensional case, even though the actual quantitative values must be 
different from (1 + 1) dimensional case. 

Furthermore, we have found the effect of zero mode to the EM form factor: 

J'^(0)-(2Pi-g)'^FM(g'). (6.16) 

The EM form factor at q^=Q using J~(0) current is obtained by 

Fm{0) = n\ / dx '-^ -2 + 1/ml , (6.17) 

2:(l-x)2(M2^^-^) ^ 
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2.0 



1.8 - 




1.2 - 




0.00 0.05 0.10 

q'[GeV'] 

Figure 6.1: (a) Normalized form factor of f+iq"^) ior K ^ tt in (1 + 1) dimension. The 
solid line is the result from the valence plus nonvalence contributions. The dotted 
line is the result from the valence contribution, (b) Normalized form factor of f-{(f) 
for X — > TT in (1 + 1) dimension. The same hne code as in (a) is used. 
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Figure 6.2: (a) Normalized form factor of f+{(f) for S — > tt in (1 + 1) dimension. 
The same line code as in Fig. 6.1(a) is used, (b) Normalized form factor of f-{q^) 
for S — > TT in (1 + 1) dimension. The same line code as in Fig. 6.1(a) is used. 
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Figure 6.3: (a) Normalized form factor of for S — > D in (1 + 1) dimension. 

The same line code as in Fig. 6.1(a) is used, (b) Normalized form factor of f-{q^) 
for S — > in (1 + 1) dimension. The same line code as in Fig. 6.1(a) is used. 
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where is the normahzation constant and the l/rn^q in Eq. ( |6.17| ) is the "zero mode" 
term. Numerically, using the previous quark masses, the effects of zero modes on 
the form factors of F^(0) and Fb(0), i.e., F^^-^-(0)/F^"^i(0) and Ff^-(0)/Ff ^^(O), 
are estimated as 16.9 and 0.75, respectively. Again, the zero mode contribution is 
drastically reduced for the heavy meson form factor. However, it gives a very large 
effect on the light meson form factors. The similar observation on the EM form factor 
was made in the Breit frame recently ||2^. In (3 + 1) dimensions, however, we note 
that the relation between the Breit frame and the Drell-Yan-West frame involves the 
transverse rotation in addition to the boost and therefore the results obtained from 
the Breit frame cannot be taken as the same with those obtained from the Drell-Yan- 
West frame or vice versa. 



6.3 Summary and Discussion 



In conclusion, we investigated the zero mode effects on the form factors of semilep- 
tonic decays as well as the electromagnetic transition in the exactly solvable model. 
Our main observation was the nonvanishing zero mode contribution to the J~ current 
and our results are directly applicable to the real (3 + 1) dimensional calculations. 
The effect of zero mode to the /-(O) form factor is especially important in the ap- 
plication for the physical semileptonic decays in the Drell-Yan-West (q'^=0) frame. 
To the extent that the zero modes have a significant contribution to some physical 
observables as shown in this work, one may even conjecture that the condensation 
of zero modes could lead to the nontrivial realization of chiral symmetry breaking in 
the LP quantization approach. The work along this line is in progress. 
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Chapter 7 

LFQM Analysis of Exclusive 
0~^0~(1~) Semileptonic Meson 
Decays 



In recent years, the exclusive semileptonic decay processes generated a great excite- 
ment not only in measuring the most accurate values of the Cabbibo-Kobayashi- 
Maskawa (CKM) matrix elements but also in testing diverse theoretical approaches 
to describe the internal structure of hadrons. Especially, due to the anticipated abun- 
dance of accurate experimental data from the i?-factories (e.g. HERA-B at HERA, 
BaBar at SLAG and Belle at KEK), the heavy-to- heavy and heavy-to-light meson 
decays such as B^D{D*), B^-n{p), D^-n{p) etc. become invaluable processes de- 
serving thorough analysis. While the available experimental data of heavy meson 
branching ratios have still rather large uncertainties , various theoretical methods 
have been applied to calculate the weak decay processes, e.g., lattice QCD [ |102| , |103| , 
|T04i |105| , |T06| , [T07| , [T08| , |109|, |lTg, |lTl|] , QCD sum rules |T|, |TI|, |lT5|, |rT|], Heavy 
quark effective theory [ |117| ], and quark models [^, 



rg, |n|, pii . 



122 , 123 , |124| , 125 |. In particular, the weak transition form factors determined by 



the lattice QCD ||109|| provided a useful guidance for the model building of hadrons. 
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making definitive tests on existing models, even though the current error bars in the 
lattice data are yet too large to pin down the best phenomenological model of hadrons. 
These weak form factors, however, are the essential informations of the strongly inter- 
acting quark/gluon structure inside hadrons and thus it is very important to analyze 
these processes with the viable model that has been very successful in analyzing other 
processes. 

In this work, we report the analysis of exclusive semileptonic decays of 0~ — > 0~ and 
0~ ^ 1^ heavy meson decays using our LFQM which has been quite successful 
in the analysis of EM form factors and radiative decays. In addition to these heavy 
meson semileptonic decays, the light-to-light weak form factor analysis such as K^^ 
decays will be disscussed by comparing with the experiment as well as many 



other theoretical models, e.g., chiral perturbation theory(CPT) ||126| , |127|| , the effec 



tive chiral Lagrangian approach ||128|| , vector meson dominance ||129|| , the extended 
Nambu-Jona-Lasino model ||130|| , Dyson-Schwinger approach ||131|| and other quark 



models [BIJ . 



The LFQM takes advantage of the equal LP time (r=t + z/c) quantization p3[ and 
includes important relativistic effects in the hadronic wave functions. The distin- 
guished feature of the LF equal-r quantization compared to the ordinary equal-t 
quantization is the rational energy-momentum dispersion relation which leads 
to the suppression of vacuum fluctuations with the decoupling of complicated zero 



modes ||26|, ^ |2^ and the conversion of the dynamical problem from boost to 



rotation |^ . Moreover, one of the most distinctive advantages in the LFQM has been 
the utility of the well-established Drell-Yan-West (^7^=0) frame for the calculation of 
various form factors By taking the "good" components of the current [J^ and 
J_l), one can get rid of the zero mode problem and compute the full theoretical 
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prediction for the spacelike form factors in q~^=0 frame. The weak transition form 
factors that we are considering, however, are the timehke > observables. Our 
method is to rely on the analytic continuation from the spacelike region to the time- 



like region calculating the "good" components of the current in the q~^=0 frame 
If we were to take the q~^^0 frame, then we must take into account the higher Fock- 
state (nonvalence) contributions arising from quark- ant iquark pair creation (so called 
"Z-graph") as well as the valence configuations. In fact, we notice that a few pre- 
vious analyses ||119|| were performed in the q^y^O frame without taking into account 



the nonvalence contributions. We find that such omission leads to a large deviation 



from the full results [HTI]. Our method is to rely on the analytic continuation from 
the spacelike region to the timelike region calculating the "good" components of the 
current in the g"'"=0 frame. 



The key idea in our LFQM for mesons is to treat the radial wave function as a trial 
function for the variational principle to the QCD-motivated Hamiltonian saturating 
the Fock state expansion by the constituent quark and antiquark. The spin-orbit wave 
function is uniquely determined by the Melosh transformation (see Appendix A). We 
take the same QCD-motivated effective Hamiltonian given by Eqs. (4.1) and (4.2) and 
the Gaussian radial wave function 0(A;^) given by Eq. (4.3) as our trial wave function to 
minimize the central Hamiltonian |T^. The model parameters of heavy-quark sector 
(c and b) such as rric, mb, Puc, Pub, etc. are then uniquely determined by the same 
procedure as the light-quark analysis [|13| discussed in Chapter 4. The procedure of 



determining model parameters constrained by the variational principle (see Eqs (D.l)- 
(D.3) in Appendix D) is shown in Fig. 7.1, where the lines of qq and qc {q=u and d) 
etc. represent the sets of {mg,mg, Pgg} and {mg,mc, Pqc}, respectively, etc. Because 
all the lines in Fig. 7.1 should go through the same point of (6=0.18 GeV^,/t = 0.313), 
the parameters of rric, mb, Puc, Pud, etc. are all automatically determined without 
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Figure 7.1: The parameters m^, rric, mj,, (3qs, l3qc, etc. satisfying variational principle 
with (a) HO and (b) Linear potential models. The qq and qc etc. represents the sets 
of {rngjing, (3gq) and {rnqjUic, (3gc) etc., respectively. 
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any adjustment. Our model parameters obtained by the variational principle are 



summarized in Tables 7.1 and 7.2 



Our predictions of the ground state meson mass spectra and the decay constants are 
summarized in Tables |7.3| and |7.4] , respectively, and compared with the experimental 
data and the lattice QCD results [|102| , |107|| . Our predictions of ground state me- 



son mass spectra agree with the experimental data |^ within 6% error. Furthermore, 
our model predicts the two unmeasured mass spectra of ^So{bb) and ^Si{bs) systems 
as M65=9295 (9657) MeV and Mfe,-=5471 (5424) MeV for the HO (linear) potential, 
respectively. Our predictions for the decay constants of heavy mesons are quite com- 
parable with the lattice QCD ||102| , |107|| anticipating future accurate experimental 
data. 

This Chapter is organized as follows: In Section 7.1, we calculate the EM form factors 
of D and B mesons as well as vr and K mesons in the Drell-Yan-West {q^=0) frame 
and compare the results with those obtained from the valence contributions in q~^^0 
frame. Comparing the form factors between the two reference frames, we quantify the 
nonvalence contributions from q~^y^O frame. In Section 7.2, we discuss the formalism 
of the form factors for various 0^ 0^ and 0^ 1^ semileptonic decays. We 
calculate the weak form factors in g"'"=0 frame (i.e. spacelike < region) and 
then analytically continue the form factors to the timelike region by changing qj_ 
to iq± in the form factors. In Section 7.3, our numerical results of the observables 
for various 0^ 0^ and 0^ — 1^ semileptonic decays are presented and compared 
with the available experimental data as well as other theoretical results. Of special 
interest, we compare our analytic solutions of the form factors with the simple pole 
approximation motivated by the vector dominance model. The monople-type form 
factors turn out to be good approximations to our analytic solutions for most decay 
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Table 7.1: Optimized quark masses mq[GeV] for both HO and linear potentials ob- 
tained from the variational principle. We also include the results from the smearing 
function(SF) instead of the Breit-Fermi contact term. q=u and d. 



Potential 


ruq 


TTLs 


rUc 




HO 


0.25 


0.48 


1.8 


5.2 


SF 


0.25 


0.48 


1.8 


5.2 


Linear 


0.22 


0.45 


1.8 


5.2 


SF 


0.22 


0.45 


1.8 


5.2 



Table 7.2: Optimized Gaussian parameters /?[GeV] for both HO and linear potentials 
obtained from the variational principle. We also include the results from the smearing 
function(SF) instead of the Breit-Fermi contact term. 



Potential 










Ps-c 


Pc-c 


Pqb 


Psb 


Pbb 


HO 


0.3194 


0.3681 


0.3419 


0.4216 


0.4686 


0.6998 


0.4960 


0.5740 


1.8025 


SF 


0.3194 


0.3703 


0.3428 


0.4280 


0.4782 


0.7278 


0.5122 


0.5980 


1.9101 


Linear 


0.3659 


0.4128 


0.3886 


0.4679 


0.5016 


0.6509 


0.5266 


0.5712 


1.1452 


SF 


0.3659 


0.4132 


0.3887 


0.4697 


0.5042 


0.6548 


0.5307 


0.5767 


1.1806 



Table 7.3: Fit of the ground state meson masses with the parameters given in Ta- 
bles |7.1| and |7.2| . Underline masses are input data. The masses of {oo-cj)) and {rj-rj') 
were used to determine the mixing angles of uo-cj) and 77-77', respectively. 



^^0 


Expt. [MeV] 


HO (SF) 


Linear (SF) 




Expt. [MeV] 


HO (SF) 


Linear (SF) 


TT 


135±0.00035 


135 (135) 


135 (135) 


P 


770± 0.8 


770 (770) 


770 (770) 


K 


494± 0.016 


470 (469) 


478 (478) 


K* 


892± 0.24 


875 (875) 


850 (850) 


V 


547± 0.19 


547 (547) 


547 (547) 


LO 


782± 0.12 


782 (782) 


782 (782) 


v' 


958±0.14 


958 (958) 


958 (958) 




1020±0.008 


1020 (1020) 


1020 (1020) 


D 


1869±0.5 


1821 (1821) 


1836 (1840) 


D* 


2010± 0.5 


2024 (2026) 


1998 (1997) 


Ds 


1969±0.6 


2005 (2004) 


2011 (2014) 




2112±0.7 


2150 (2150) 


2109 (2108) 


Vc 


2980± 2.1 


3128 (3123) 


3171 (3173) 


J/V' 


3097±0.04 


3257 (3244) 


3225 (3221) 


B 


5279± 1.8 


5235 (5231) 


5235 (5237) 


B* 


5325± 1.8 


5349 (5349) 


5315 (5314) 




5369±2.0 


5378 (5372) 


5375 (5376) 


(bs) 




5471 (5466) 


5424 (5423) 


m 




9295 (9353) 


9657 (9651) 


T 


9460± 0.21 


9558 (9459) 


9691 (9675) 
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Table 7.4: Decay constants [MeV] and charge radii [fm^] for various heavy 
pseudoscalar and vector mesons. Note that the fp and fy are obtained by 
(0|g27'^75Q'i|-P) = ifpP'^ and (0|g27'^5'i|^) = Myfv^^, respectively, to compare with 
the lattice results. In the above definitions of the decay constants, our f^,, for example, 
obtained in Chapter 4 becomes 130.7 (129.8) for the HO (linear) parameters. 



References 


fD 


Id- 




fB 


fB- 






HO 


179.7 


211.6 


218.6 


160.9 


173.0 


207.0 


Li 


icar 


196.9 


238.9 


233.1 


171.4 


185.8 


203.9 


LAT 


ioi] 


200±30 




220±30 


170±35 




195±35 




m 


195 ± ll^giV^ 




213 ± 


159 ± Ullf^o^^ 




175 ± 10+^g+f 


Expt. |4£|] 


< 219 




137 - 304 








References 




^dO 


DJ 






'bo 




HO 


0.182 


-0.309 


0.106 


0.420 


-0.208 


-0.081 


Linear 


0.176 


-0.301 


0.101 


0.438 


-0.217 


-0.083 



processes except heavy-to-light decays, e.g., B ^ tx and B ^ p. We also quantify 
for 0^ 0^ decays the nonvalence contributions from q~^y^O frame as in the case of 
the EM form factor analyses. Summary and discussion of our main results follow in 
Section 7.4. In the Appendix F, we show the derivation of the matrix element of the 
form factors for 0^ 0~ semileptonic decays both in the g^=0 and q'^y^O frames. 



7.1 EM Form Factors of B and D Mesons 



The EM form factor of a pseudoscalar meson is defined by the the relation 

(P2|Je"m|A) = (Pr + P2Wg'), (7.1) 

where Jl^^ is the EM current for the quarks and Pi{P2) is the four- momenta of the 
initial(final) meson. In the calculations of the hadronic matrix elements in the LP 
frame, one usually use q^=P^ — P'^=0 frame. This leads to = q^q~ — q± = —q± = 
—Q^ implying a spacelike momentum transfer. 

In this work, however, we calculate the valence contributions from both g+=0 and 
g+7^0 frames and compare the difference in spacelike region so that we can indirectly 
estimate the nonvalence contributions from q^^O frame. The form factor from the 
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valence contribution in g"''=0 frame is an exact solution since the "Z" -graph does 
not contribute at all in this q~^=0 frame as far as the "good" components of the 
currents are concerned, i.e., or PB| , In the standard g+=0 frame, the 

quark momentum variables are given by In the standard q'^=0 frame, the quark 
momentum variables are given by 

pt = (l-x)P+, p^ = xP^, 

Pl_L = (1 - x)Pi_L + k_L, Pq± = xPi_L - k_L, 

pt = {1-X)P+, P'g+=XP+, 

P2± = (1 - x)P2± + k'x, p'^-^ = xP2± - k'x, (7.2) 
which requires that Pq=p'q and Pg±=p'g_L. 

Taking Pi_l=0, P2±=— q±, we obtain from the matrix element of the " + " component 
of the current J"*": 

Fm = e,£dxJd'k^^^^^<l,;ix,k\)MxM) 
+ k_L ■ k'_L 

X , _ , = + eg(q g of the first term), (7.3) 
^A^ + kl^A^ + k'l 

where k'_L=k_L — xqj_, A=xmq + (1 — x)mq and eg(eg) is the charge factor of the quark 
(antiquark). 

The EM form factor of a pseudoscalar meson obtained from the valence contributions 
in q'^^Q frame is given in appendix F. Because of the deficit of the nonvalence (or 
"Z-graph") contributions in q'^^Q frame, the EM form factors in g+=0 [Eq. (|7.3|) ] and 
g"'"7^0 frames [Eq. (F.IO)] are different with each other for the non-zero momentum 
transfer. Thus, we can estimate quantitatively the effect of nonvalence contributions 
by comparing Eq. (|7.3| ) and Eq. (F.IO). In Figs. 7.2 and 7.3, we show the EM form 
factors of the D and B mesons and compare the results of g^=0 frame with those 
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of Q'^T^O frame, respectively. For comparison, we also calculated the EM form 
factors FT,{q^) and Fx^q^) in the spacelike region using both g"'"=0 and g"^7^0 frames to 
estimate the nonvalence contributions in q'^^^ frame. As shown in Figs. 7.4 and 7.5 
for FT^{q^) and Fj^iq^)-, respectively, our predictions in g+=0 frame are in a very good 
agreement with the available data |7^, ^ while the results for q^^Q frame deviate 
from the data significantly. The deviations represent the nonvalence contributions in 
g"'"7^0 frame. However, the deviations are clearly reduced for Fxiq^) because of the 
large suppression from the energy denominator shown in Fig. 1.1 in Chapter 1 for the 
nonvalence contributions. The suppressions are much bigger for the heavier mesons 
such as D (see Fig. 7.2) and B (see Fig. 7.3). Especially, for the B meson case, the 
nonvalence contributions are almost negligible up to = — ~ 10 GeV^. This 
seems quite natural because it is much harder to create the heavy quark pair than 
the light one in "Z-graph". 

7.2 Weak Form Factors for 0~ 0~ Transitions 

The matrix element of the current J^=q2l^Qi for Q~{Qiq) 0^{q2q) decay is given 
by two weak form factors /+ and /_, viz., 

(P2ig27^giiPi) = uiq'm + P2Y + (7.4) 

where q^={Pi — is the four-momentum transfer to the lepton and ml < q^ < 
(Ml — 1^2)"^. Equation (|7.4|) can also be parametrized into two other different ways 



{P2\q2rQi\Pi) = ^M,M, 



h+{q^){vi+V2r + h^{q^){v^-V2y 



(Pi + ^2)^ - \^ \ ' 



Q 

2 T\J2 



— /If ^ 

+ Uq')^^^q'. (7.5) 
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Figure 7.2: The EM form factor of D meson for the spacehke = —(f region. The 
thick sohd (hght sohd) and thick dotted (hght dotted) hnes are the results from the 
q^—Q (?^7^0) frame for the HO and hnear potential parameters, respectively. 7.3: 
The EM form factor of B meson for the spacelike — —q^ region. The same line 
code as in Fig. 7.2 is used. 
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Figure 7.4: The EM form factor of pion for low = —q^ compared with data [70|. 
7.5: The EM form factor of kaon compared with data [Q. The same hne code as in 
Fig. 7.4 is used. 
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where Vi = Pi /Mi and satisfy the following relations: 



/+(0) = /o(0), Mq') = U{q') + M2_^| /-(g^)- (7.6) 

The second parametrization uses the helicity basis, where /+ is related to the exchange 
of a vector particle and /o to the exchange of a scalar particle. In the heavy-quark 



limit Ml 2 oo, heavy-quark symmetry requires that | 132 | 



h+{q')=av,-V2), h.{q') = 0, (7.7) 

where ^(fi ■ V2) is the universal Isgur-Wise (IW) function[| normalized to unity at the 
point of equal velocities, ^(1)=1. 

Since the lepton mass is small except in the case of the r lepton, one may safely neglect 
the lepton mass in the decay rate calculation of the heavy-to-heavy and heavy-to-light 
transitions. However, for decays, the muon(yu) mass is not negligible, even though 
electron mass can be neglected. Thus, including nonzero lepton mass, the formula 
for the decay rate of Kf^ is given by |^ 



dr a 



dq^ 24773' wm^2i j ^ y ^ 



X 



where Gp is the Fermi constant, Vq^q^ is the element of the CKM mixing matrix 
and the kinematic factor Kf{q^) (=|P2|) in the decaying meson rest frame (Pi=0) is 
given by 

Kf{q^) = ^ V(Mf + Ml - g2)2 _ ^m^mI (7.9) 

Note that the form factor fo{q^) does not contribute to the decay rate in the limit 
of massless leptons. Since our analysis will be performed in the isospin symmetry 



^ From Eqs. ( |7.4| )-( [7/?l) , one can easily find that the form factors f+i'f') f^^d fo{q^) are related 
to the IW function via ^Vi ■ ^2) = ^^^ /+(?')= [i.^./^m.+m.)^] ■ That means the 

dependence of /+ is different from that of /o by an additional factor. 
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{mu=md) but SUf{3) breaking {ms^mu(d)) limit, we do not discriminate between the 
charged and neutral kaon weak decays, i.e., f^"=f^^ . For Ki^ decays, the three form 
factor parameters, i.e., A+, Aq and ^a, have been measured using the following linear 



parametrization 49 



-,2 



f±{q') = f±iq' = m'l) (l + j , (7.10) 

where A±^o is the slope of /±^o evaluated at q^=mf and ^A=f-/ 



In LFQM, the matrix element of the weak vector current can be obtained by the 
convolution of initial and final LP meson wave functions in g'^=0 frame where the 
decaying hadron is at rest: 

{P2\q2l^Qi\Pi) 







xTr 



dx V dx 2(1 -x) Ui=i jMl-{nn-m^Y 



75(^2 + "^2)7''(A + "^1)75(^5 - mq 



(7.11) 



where MlQ={k']_ + m^)/(l — x) + {k\ + m?)/x. Using the matrix element of the "+" 
component of the current, J"*", we obtain from Eqs. ( |7.4| ) and ( [7.11j ) the form factor 
f+{q\) as follows 



f+iq±) = I dxjd kxW^W^02(a;,k^)0i(x,k 



Al + kl^Al + k'l 

(7.12) 



where q\=—q'^, Ai=miX + rriqil — x) and k'_L=kj_ — xq_L. As we discussed in the 
introduction, we need the "±" component of the current, J^, to obtain the form 
factor /_(gi) in Eq. (|7^), viz., 

(P2|g2(q± ■ ii.)Qi\Pi) = ql[f-{ql) - UiqDl (7-13) 

after multiplying on both sides of Eq. ( |7.4| ). The left-hand side (LHS) of Eq. ( [7.13| ) 
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is given by 



{P2\q2{q± ■'y±)Qi\Pi) = - dx d ——\ — , r 

h J \ dx \ dx 2JAl + klJAl + k'l 



X Tr 



75(^2 + "^2)(q± ■ 1±){A + "^1)75(^5 - 



• (7.14) 



Using the quark momentum variables given in Eq. (7^), we obtain the trace term in 
Eq. ( 17.141) as follows: 



Tr 



75(^2 + "^2)(q± ■ li.){A + "^1)75(^5 - mq 

{Al + k'l) 



x{l — x) 
+ [(mi - 1712^ + ql]k± ■ q± \. 



x{l — x) 



k_L ■ q_L 



(7.15) 



The more detailed derivation of Eqs. ( [7.12| ) and ( [r.l4| ) are presented in Appendix 
F. Since both sides of Eq. ( |7.13| ) vanish as 0, one has to be cautious for the 

numerical computation of /_ at q^=0. Thus, for the numerical computation at q^=0, 
we need to find an analytic formula for /-(O). In order to obtain the analytic formula 
for the form factor /-(O), we make a low q'j_ expansion to extract the overall q'j_ from 
Eq. ( [7.14|) . Then, the form factor /_(0) is obtained as follows: 



/_(0) = /+(0) + fdx I d^kn 



Idkin dk2n x(j)l{x,k±)(j)i{x,k±) 
dx V dx /Al + klJAl + kl 



Cti{Cji — Cj2 + Cm + Ciij + Ct2 



k] cos^ d + G 



T3 



(7.16) 



where the angle (p is defined by k_|_ ■ q_|_ =|k_|_| |q_|_| cos and the terms of C's are given 
by 



Cji 
Cm 

Cr 



'l-xm + f3l)Ml 

1 



1 



-1 



i-[K-^i)/M^]^ 4 



Pi 



M20 - {m2 - mqY Mfo - (mi - mqf 



A{l-xm + Pl) 



mi, — m?^ 2 



m| — m 



2\ 2- 
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Cti 

Ct2 



1 



{Al + Al + 2kl) + (mi 

Ct3 



7712) 



xP! Al + kl 



Pi + Pi 



x{l — x) 



(7.17) 



The form factors /+ and can be analytically continued to the timelike q^>0 region 
by replacing q± by iq± in Eqs. ( |7.12|) and (|7.13|) . Since /-(O) in Eq. ( [7.16| ) is exactly 
zero in the SUf{3) symmetry ||131|| , i.e., m„(rf)=ms and Py^^=Pus=Pss, one can get 
/+(?i)=-^vr('?i) for the 7r+ — 7r° weak decay (TTga), where -^^(q'i) is the EM form 
factor of pion, and f-{q^) = because of the isospin symmetry. 



7.3 Weak Form Factors for ^1 Transitions 

The form factors for the semileptonic P{Qiq) —>■ V{q2q) decay of a pseudoscalar to a 
vector meson is given by 



{V{P2,e)\q2Y{l-^,)Q,\P{P,)) 

— Ojp *upppa 



-£*^(Mi + M2)Ai(g2) + (P, + p^Y^e* ■ q)- "^'^^'^ 



M1 + M2 



where 



W{e* ■ q)^-^[A,{q') - A,{q% (7.18) 



, , Ml + Ms , , Ml - M2 , , , 



Alternatively, one can also use another set of the form factors in the following way: 

{V{P2,e)\q2i^{l-^,)Q,\P{P,)) 
= le^.p^e^'P^P^giq^) 
+e*^fiq^) + P^ie* ■ P)a+{q^) + q^^i^e* ■ P)a4q^), (7.20) 



Ill 



where P=Pi + P2- The relationship between the two sets of the form factors is given 
by 



iM, + M2)g{q^ 
1 



2M2 
fiq' 



[f(q^) + (Mt-Mi)a4q') + q'a4q% 



(Mi + M2)a+(g2). 



(7.21) 



The form factors proportional to g^, i.e., ^0,^3 and a_, do not contribute to the 
decay rate in the limit of massless leptons and the corresponding differential decay 
rate for P — > Viu^ is given by 



dV _ Gl 
1^ ~ 967r3Mf 



VQ.Q,?Kf{q')q' 



H4q')\' + \H4q 



2m2 



\Hoiq')\' 



(7.22) 



where Kf{q^) is the kinematic factor given in Eq. ( [7.9|) and the helicity amplitudes 
are given by 



2M2V2 



{M^ - Ml - q^){M, + M2)Ai(g2) - 4 ' ^ ^' A^^q^) 



M1 + M2 



= (M, + M2)A,iq^) T 2 ^'f!^f} v{q% 



(7.23) 



M1 + M2 

Note that the form factor Ai{q^) dominates at large q^ all three of the helicity am- 
plitudes because of Kf{q^)=Q at ql^^^={Mi — M2Y. Thus, to obtain the decay rates 
for 0^ 1^ processes in the massless lepton limit, we need to calculate the form 
factors of V, Ai, and A2 (or g, f, and a+). The form factors g, a+, and / can be cal- 
culated using the '+'-component of the current in g"*" = frame, i.e., the transverse 
decay modes of the vector and axial-vector currents determine g^ and a_|_ and the 
longitudinal mode of the axial- vector current determines /. 



The use of the longitudinal polarization vector gives zero on both sides of Eq. ( 7.18 ). 
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and 



The explicit expressions of the form factors g, a+ and / are given by 

[rrii — 1712 



Idkin \dk'2n x(f)2ix, k'_L)0i(x, kj 



(7.24) 



-k_L ■ q_L + 



M20 + m2 + mq 



■^2 _ (k±-q±^'^ 



cix / (fk ■ ' ■ ' ^<^2(a;, k'x)0i(2;, k^) 



^a; V dx + fclV^i + '^l 



X ^ (1 - 2x)A - ^^^-^^[(1 - 2x) A - vAa] 



xqi 



M20 + m2 + m-q 



k'x-kx + ^ii32; 



(7.25) 



/(g2) = (M|-M2 + g>+(g2 



ii(a;, kj 



-2M2 /rf^fcx dk^^ix^ 



x<^2a;(l-a;)^iM2o + 



1 - 2x)M2o + 1712- rUg 
M20 + m2 + mq 



(7.26) 



[k^-kx+^ii32] 



where Bi=—xmi + (1 — x)mg. Accordingly, the form factors V(5'^), ^i(q'^) and ^2(0'^) 
can be obtained from Eq. ( [7.21|) . As in the case of 0^ 0^ transitions, we now 
analytically continue the form factors ^(q'^), ^i(q'^) and A2{q'^) to the timelike 
region by changing q± to iq± in Eqs. ( |7.24]) - ([7726|) . Of special interest, we also compare 



our analytic solutions with the following form factor parametrization @, |123|| : 

^(,2) = m 

J l^q2/A2 + s2qyAr 
where S2=±l and the parameters Aj are determined by the calculation of the appro- 



(7.27) 



priate derivatives of J^(g2) at q^=0 using Eqs. (|7l2|) and (^!24|)-(^). 
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7.4 Numerical Results 



For the numerical calculations for various — (1 ) semileptonic decays, we used 



the quark model parameters given in Tables [7.1| and |7.2| . Since the smearing effect is 
very small, we do not include the results obtained from the parameter sets of the SF 
in our numerical analyses. 

7.4.1 K^^, decays 

Our predictions of the parameters for K^^^ decays in g''"=0 frame, i.e., /+(0), A4., 
Ao, (r2)^,=6/;(0)//+(0)=6A+/M2+, and U=f-I f+\e='m^, given by Eq. (Ol), are 
summarized in Table [7.5| . We do not distinguish K^.^ from i^^a in the calculation 
of the above parameters since the slopes of f± are almost constant in the range of 
^1 < Q'^ < However, the decay rates should be different due to the phase 

space factors given by Eq. ( |7.8| ) and our numerical results for r(_K'e3) and Vi^K^'^) 
in g^=0 frame are also presented in Table |7.5| . Our results for the form factor /+ 
at zero momentum transfer, /+(0)=0.961 [0.962] for the HO [linear] parameters, are 
consistent with the AdemoUo-Gatto theorem [|133|] and also in an excellent agreement 



with the result of chiral perturbation theory [ [L26|| , /+(0)=0.961±0.008. Our results 
for other observables such as A+, S,Ai and T{Ki'i) are overall in a good agreement with 



the experimental data [^. We have also investigated the sensitivity of our results 
by varying quark masses. For instance, the results^ obtained by changing the strange 

quark mass from ms=0.48 GeV to 0.43 GeV (10% change) for the HO parameters 

^ Even though we show the results only for the HO parameters, we find the similar variations for 
the linear parameters; i.e., the posivive sign of Aq can be obtained when ms/mK<1.8 for both HO 
and linear paramters. In addition to the observables in this work, our predictions for Jk, 'r'^+, and 
r^o in [H are changed to 108 MeV (1% change), 0.385 fm^ (0.3% change), and -0.077 fm^ (15%), 
respectively. 



114 



Table 7.5: Model predictions for the parameters of i^£3 decay form factors obtained 
from g"*" = frame. As a sensitivity check, we include the results in square brackets 
by changing m<j=0.48 to 0.43 GeV for the HO parameters. The CKM matrix used in 
the calculation of the decay width (in units of 10^ s"^) is = 0.2205 ± 0.0018 



Observables 


HO 


Linear 


Other models 


Expt. [i9|] 


/+(0) 


0.961 
[0.974] 


0.962 


0.961±0.008'', 0.952^0.98^0.933 




A+ 


0.025 
[0.029] 


0.026 


0.03^,0.033^0.025'' 
0.028^0.018■'',0.0199 


0.0286±0.0022[X+] 
0.0300±0.0016[X03] 


Ao 


-0.007 
[+0.0027] 


-0.009 


0.017±0.004'' 0.013^0.0" 
0.0026^-0.0024/ -0.0059 


0.004±0.007[is:+3] 
0.025±0.006[is:^3] 


u 


-0.38 
[-0.31] 


-0.41 


-0.164±0.047'' -0.24^-0.28" 
-0.28^-0. 25-'', -0.289 


-0.35±0.15[i^+3] 
-0.11±0.09[i^03] 


(r)^K(fm) 


0.55 
[0.59] 


0.56 


0.61^0.57^ 0.47^0.48f 




r(i^e"3) 


7.30±0.12 
[7.60±0.12] 


7.36±0.12 




7.7±0.5[if^'3] 




4.57±0.07 
[4.84±0.08] 


4.56±0.07 




5.25±0.07[X;]3] 



" Ref. ^2§° Ref. |27[,'= Ref. [[128 



Ref. il29i, " Ref. |130|. ^' Ref. il3l|,g Ref. 18 



are included in Table |7.5| . As one can see in Table |7^, our model predictions are 
quite stable for the variation of except Aq, which changes its sign from —0.007 
to +0.0027. The large variation of Aq is mainly due to the rather large sensitivity of 
/_(0) (18% change) to the variation of m^. Similar observation regarding on the large 
sensitivity for Aq compared to other observables has also been reported in Ref. [|130|1 
for the variation of quark masses. As discussed in Refs. ||131|| and ||134|| , /-(O) is 
sensitive to the nonperturbative enhancement of the SU(3) symmetry breaking mass 
difference — mu(^d) since /-(O) depends on the ratio of and mu{d)- 

Of special interest, we also observed that the nonvalence contributions from g^T^O 
frame are clearly visible for A+, Aq and S,a even though it may not be quite significant 
for the decay rate T{K(^^). Our predictions with only the valence contributions in g^T^O 
frame are /+(0)=0.961 [0.962], A+=0.081 [0.083], Ao=-0.014 [-0.017], ^=-1-12 
[-1.10], r(is:e3) = (8.02[7.83]±0.13)xl06 s^^ and r(is:^3) = (4.49[4.36]± 0.13)xl06 s"^ 
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Figure 7.6: The form factors f+{q^) for the K ^ ir transition in timehke momentum 
transfer > 0. The sohd and dotted hues are the results from the q~^=0 and q~^j^O 
frames for the HO and hnear parameters, respectively. The differences of the results 
between the two frames are the measure of the nonvalence contributions from 
frame. 

for the HO [linear] parameters. Even though the form factor /+(0) in q~^^0 frame 
is free from the nonvalence contributions, its derivative at q^=0, i.e., A+, receives 
the nonvalence contributions. Moreover, the form factor f-{q^) in q'^^Q frame is 
not immune to the nonvalence contributions even at g^=0 |]2B|. Unless one includes 



the nonvalence contributions in the q'^^^ frame, one cannot really obtain reliable 
predictions for the observables such as A+, Aq and for K^^ decays. 

In Fig. 7.6, we show the form factors /+ obtained from both g+=0 and q'^^^ frames 
for 0<g^<(Mx — Mtt)^ region. As one can see in Fig. 7.6, the form factors /+ obtained 
from g+=0 frame (solid lines) for both HO and linear parameters appear to be linear 
functions of q^ justifying Eq. ( [7.10| ) usually employed in the analysis of experimental 



data 491. Note, however, that the curves without the nonvalence contributions in 
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Figure 7.7: The decay rates dV/dq^ of -ft'e3(solid line) and -K'^3(dotted line) for the 
HO parameters in g"'"=0 frame. The dot-dashed line is the contribution from the term 
proportional to /o in Eq. ( [7. 81 ) for K^^, decay. The results for the linear parameters 
are not much different from those for the HO parameters. 

g"'"7^0 frame (dotted lines) do not exhibit the same behavior. In Fig. 7.7, we show 
dV /dq^ spectra for K^.'^, (solid line) and K^^ (dotted line) obtained from g'^=0 frame. 
While the term proportional to /o in Eq. (|7.8| ) is negligible for i^es decay rate, its 
contribution for K^^^ decay rate is quite substantial (dot-dashed line). Also, we show 
in Fig. 7.8 the form factors f+iq^) (solid and dotted lines for the HO and linear 
parameters, respectively) at spacelike momentum transfer region and compare with 
the theoretical prediction from Ref. ||13CI|| (dot-dashed line). The measurement of this 
observable in < region is anticipated from TJNAF [[L3CI|| . 

We have also estimated the zero-mode contribution by calculating the " — " compo- 
nent of the current. Our observation in an exactly solvable scalar field theory was 
presented in Chapter 6 ^6). Using the light-front bad current J~ in g+=0 frame, 
we obtained /_(0)=12.6 (18.6) for the HO (linear) parameters. The huge ratio of 
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Figure 7.8: The form factors f+{q^) for the K ^ ir transition in spacehke momentum 
transfer — < 0. The sohd and dotted hues are the results from the HO and hnear 



parameters, respectively. The dot-dashed line is the result from Ref. ||130 |. 



/_(0)|j-//_(0)|jj^~ —36 (—48) for the HO (linear) parameters is consistent with our 
observation in Chapter 6 [^. We also found that the zero-mode contribution is 
highly suppressed as the quark mass increases. 

7.4.2 D and Ds decays 



(1) D — ^ TT e^Ue'. Our predicted decay rate for D n in q^=0 frame is r{D^ —>■ 
7r-e+z/e)=0.110 (0.113)|KdP ps~^ Using the lifetime r^o = 0.415 ± 0.004 ps and 
I V^rf|=0.224±0.016 E9i, we obtain the following branching ratio 



Br(D'' 



71 e Up 



(2.30 ±0.33) X 10"^ (HO), 

(2.36 ± 0.34) X 10"^ (Linear), (7.28) 

7r-e+z/e)= (3.9 ±tli ±0.4) x 



while the experimental data |^ reported Brexp.(-D° 
10^'^. Even though our predicted decay rates are rather smaller than the central 
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Table 7.6: 

D^o \K* 



Summary of the parameters for the various form factors for D—>-tt [K] and 





D^TT [K] 


D^p 


K*] 


F 


/+ 


V 


Ai 


A2 


HO 


F{0) 


0.593 [0.724] 


0.709 [0.809] 


0.510 [0.640] 


0.342 [0.450] 




1.73 [1.25J 


1.006 [1.029] 


0.572 [0.698] 


0.419 [0.534] 


Ai [GeVJ 


1.95 I2.06J 


1.94 [2.06] 


3.37 [3.41] 


2.45 [2.41] 


S2A2 [GeV] 


2.88 [3.04] 


2.86 [3.04] 


-5.05 [-5.05] 


3.18 [3.43] 




2.71 [1.30J 


1.039 [1.029] 


0.570 [0.698] 


0.428 [0.534] 


^ A.1 A2 (^max) 


1.71 [1.21] 


1.007 [1.042] 


0.571 [0.697] 


0.420 [0.538] 


Linear 


F(0) 


0.618 [0.736] 


0.739 [0.822] 


0.517 [0.638] 


0.311 [0.414] 




1.57 [1.23J 


1.002 [1.024] 


0.563 [0.683] 


0.357 [0.475] 


Ai [GeVJ 


2.06 [2.161 


2.05 [2.15] 


3.99 [3.85] 


2.92 [2.67] 


S2A2 [GcV] 


3.17 [3.291 


3.26 [3.32] 


-4.62 [-4.61] 


3.27 [3.76] 


^ Ai (9max) 


2.10 [1.231 


1.036 [1.033] 


0.559 [0.683] 


0.362 [0.475] 








^KiK2 ("Jmax) 


1.61 [I.ISJ 


1.018 [1.024] 


0.561 [0.681] 


0.357 [0.477] 


LAT 1 

1 
1 


lilf 

Loe 
llf 


1 


m 
m 

F(0) 
F{0) 


0.65(10)[0.73(7)] 
■ ■ ■ [0.78±0.08] 
0.58(9) [0.63(8)] 
0.61(i?)[0.67(«7)] 


1.1(2)[1.2(2)] 
••• [1.08±0.22] 
0.78(12) [0.86(10)] 
0.95(29)[1.01(?«)] 


0.65(7)[0.70(7)] 
••■ [0.67±0.11] 
0.45(4) [0.53(3)] 
0.63(6)[0.70(In)] 


0.55(10)[0.6(1)] 
■ • ■ [0.49±0.34] 
0.02(26) [0.19(21)] 
0.51(}«)[0.66(i«)] 


SR 1 
1 


IV 

Llf 


1 
1 


F{0) 
F{0) 


0.5±0.1 [■ ■ ■] 
■ ■ ■ [0.60±0.15] 


1.0±0.2 [■•■] 
■■■ [1.10±0.25] 


0.5±0.2 [■•■] 
■■■ [0.50±0.15] 


0.4±0.1 [■ ■ ■] 
■ ■ ■ [0.60±0.15] 


QM [ 
1 


11!- 
12f 


1 


m 
m 


0.69 [0.76] 
0.67 [0.78] 


1.23 [1.23] 
0.93 [1.04] 


0.78 [0.88] 
0.58 [0.66] 


0.92 [1.15] 
0.42 [0.43] 


Expt. 1 


L35 


1 


F(0) 


■ ■ ■ [0.76±0.03] 


■ ■ ■ [1.07±0.09] 


■■■ [0.58±0.03] 


■ ■ ■ [0.41±0.05] 



value of the data, one shoud note that the number of events for the D ^ n data is 
currently very small compared to other processes We also found that the decay 
rate obtained from the valence contributions in q^^O frame is about 18 (17)% larger 
for the HO (linear) parameters than that obtained from g+=0 frame. 

In Figs. 7.9(a) and 7.9(b), we present the dependence of the form factors (thick 
lines) f^^{q^) and f^^iq^) obtained from g+=0 frame (thick solid lines for HO 
paramters and dotted lines for linear parameters), respectively. For comparison, we 
also show the results (thin lines) obtained from the valence contributions in q'^^Q 
frame as well as the lattice calculations ||109[ . As on can see in Figs. 7.9(a) and 
7.9(b), while the differences from the choice of the model parameters (i.e., HO or 
linear) are not significant, the ones between the two frames (i.e., g"'"=0 and q'^^Q 
frames), which are the measure of the nonvalence contributions from g"*"7^0 frame, 
are quite substantial, especially for the form factor f^^iq"^) case. Our predictions in 
g"'"=0 frame are also comparable to the the lattice results 
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Table 7.7: Comparison of the form factors f-{q^) at g^=0 between g"'"=0 and q^^^ 
frames. The differences between the two frames are the measure of the nonvalence 
contributions in q^^^ frame. Note that c^=| cos5p| and c^' = | sin5p|. 



Models 




D ^-n 


D^K 


Ds^V 




HO 


g+=0 
g+^0 


-0.428 
-0.600 


-0.400 
-0.723 


-0.416c^ 
-0.764c^ 


-0.416c^. 
-0.431CV 


Linear 


g+=0 
g+^0 


-0.454 
-0.625 


-0.422 
-0.719 


-0.460c^ 
-0.779c^ 


-0.460c^. 
-0.449c^/ 



In Table [7.6| , we compare our form factor f^^{q^) at both g^=0 and q'^=q'^s.x with the 
simple pole parametrization given by Eq. (|7.27|) as well as other theoretical predic- 
tions. It is quite interesting to note that our analytic solution /+(q'max)= l-'''^ (1-57) 
for the HO (linear) parameters is well approximated by Eq. ( |7.27| ), ^AiA2(Q'max)= ^■'^^ 
(1.61), but not fitted by the simple pole approximation, jFAi(g^)=F(0)/(l — g^/A^) = 
2.71 (2.10). We also show in Table [7.7| the form factors f^^{0) obtained from both 
^"^"=0 and q^^O frame. As one can see in Table [77^ , the form factor f^'^{0) obtained 
from q^^O frame for the HO (linear) parameters is about 40 (38)% larger than that 
obtained from g+=0 frame. 

(2) K^i^u: Our predicted decay rate for D K in q'^=0 frame is r{D^ — 

K-e+Ue)= 8.26 (8.36) |V;/ x lO'^ ps"^ and T{D^ K-^i+v^)= 6.40 (6.43) x 
10~^ ps~^ for the HO (linear) parameters. Using | Vcs|=1.04±0.16 [^, our predictions 
for the branching ratio are given by 

Br(L>° /Sf-e+z/,) = (3.71 ± 1.14)% (HO), 

= (3.75 ±1.16)% (Linear), 

Br(D° ^ /r-/i+z/^J = (2.87 ±0.88)% (HO), 

= (2.89 ±0.90)% (Linear), (7.29) 

while the experimental data are BrgxpX-D'^ — > -ft'^e+z/e) = (3.66±0.18)% and Broxp.(-D'^ — 
K^/i+z/^) = (3.23±0.17)%, respectively. Our results for the electron (e) decay mode are 
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in a good agreement with the experimental data for both HO and hnear parameters. 
For the muon (/i) decay mode, our values are rather smaller than the central value 
of the experimental data, neverthless, are quite comparable with the data within the 
given uncertainties of CKM matrix, Vcs- The result of the decay rate obtained from 
the valence contributions in q^^O frame is about 16 (15)% larger for the HO (lin- 
ear) parameters than that obtained from g"^=0 frame. Our value of /:|^^(0) =0.724 
(0.736) for the HO (linear) parameters is very close to the available experimental 
data /^^^*'(0)=0.7±0.1. Unlike the D ^ n transition, our analytic solution 
/+'^(Q'max)= 1-25 (1.23) for the HO (linear) parameters is well approximated by both 
Eq. ( |7.27| ), ^AiA2(Q'max)= 1-21 (1-18), and the simple pole approximation, ^Ai('?max) = 
1.30 (1.23). This implies that the simple pole approximation is enough to fit the 
form factor f^^{q^) and one may neglect the A2 and higher order contributions in 
Eq. (7.27). Moreover, our simple pole masses Ai=2.06 and 2.16 GeV for both HO and 



linear parameters, respectively, are in good agreement with the value of 2.11(=Z}*^) 
GeV expected from the closest resonance with the proper quantum number J^=l^ . 
As one can see in Table |7.7| , the form factor f!^^{0) obtained from q^^O frame for 



the HO (linear) parameters is about 80 (70)% larger than that obtained from ^"^=0 
frame. 

In Figs. 7.10(a) and 7.10(b), we present the dependence of the form factors f^^{q^) 
and fQ^{q^) obtained from both g+=0 (thick lines) and Q'^t^O (thin lines) frames, 
respectively, and compare with the available experimental data (•) as well as the 



lattice QCD calculations | 109 |. Again, the nonvalence contributions from q^^Q frame 
are quite sizeable especially for f^^^q^) case. Our results for both f^^{q^) and 
fo^i^f) ^'"''=0 frame are overall in a good agreement with the lattice calculations 
in Ref . ITOl . 
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Figure 7.10: (a) The form factor f+{q^) for the D ^ K transition compared with the 
experimental data (•) and the lattice calculation ||109|| . The same line code as in 
Fig. 7.2 is used for our results, (b) The form factor fo{q^) for the D ^ K transition. 
The same line code is used as in fa). 
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Figure 7.11: The form factors V{q'^) (solid line), v4i(g^) (dotted line) and A2{q'^) 
(long-dashed line), where the thick (thin) lines represent the HO (linear) parameters, 
ioT D p transition compared with the lattice calculation ||110|| of v4i(g^) (data). 
7.12: The form factors ^(i?^), Ai{q'^) and ^2(0'^) for D K* transition compared 
with the lattice calculations ||105|| (square) and ||11CI|| (circle) of Ai{q'^). The same line 
code as in Fig. 7.11 is used. 
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(3) — > p^i'^uf. Our predicted decay rate for D ^ p in g+=0 frame is r{D^ 
p~i~^i'i)= 6.38 (6.61)|Vcdp X 10^^ ps~^ for the HO (linear) parameters. We also obtain 
the ratio Tl/Tt, where the subscript L represents the longitudinal polarization and T 
transverse polarization, for the HO (linear) parameters as 1.39 (1.48), which is quite 



comparable to the QCD sum rule result: 1.31±0.11 in Ref. [[L14|| . Our predictions for 
the branching ratio given by 

Bt{D'^ p-i+ue) = (1.33 ±0.19)% (HO), 

= (1.38 ±0.19)% (Linear), (7.30) 



are overall in a good agreement with the available experimental data [^: Brgxp. = 
(1.87 ± 0.9)%. We show in Table 0| our values of V, Ai and A2 at both q^=0 and 



q^=qa\a.yL ^^^^ comparc with Eq. ( |7.27|) as well as other theoretical predictions. As one 
can see in Table |7^, our anal5d;ic solutions for the form factors V{q^^-^= 1.006 (1.002), 
Ai(g^^^)=0.572 (0.563) and A2(g^^x)=0.419 (0.357) for the HO (linear) parameters 
are well approximated by Eq. (17:2^) , yA,A,(gLx) = l-007 (1.018), AiA,A,(g^,J =0.571 
(0.561) and A2K^K^{q^^-^=QA2Q (0.357). Our analytic solutions are also well fitted 
by the simple pole formular, ■l/Ai(gLx)=l-039 (1.036), AiAi(gLx)=0-570 (0.559) and 
^2Ai(Q'max)=0-428 (0.362). This again means that the simple pole formular may be a 
good approximation to fit the form factors V , Ai and A2. In Fig. 7.11, we present our 
q^ dependence of the form factors V{q^) (solid line), Ai{q^) (dotted line) and ^2(0'^) 
(long-dashed line) obtained from g^=0 frame, where the thick and thin lines represent 
the HO and linear parameters, respectively. We also compare our Ai{q^) with the 
available lattice result | [110| j (data) and show they are overall in a good agreement 
with each other. 

(4) — s> K*''i^i'£: Our predicted decay rate for D K* in g+=0 frame is 
T{D^ K*-i+ue)=5M (5.85) x 10"^ pg-i for the HO (linear) parameters. We 
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obtain the ratio Tl/Tt as 1.37 (1.42) for the HO (hnear) parameters, while the avail- 
able experimental data reported as 1.18±0.18±0.08 (E653) [|13|] and 1.20±0.13±0.13 
(E687) | |137|| . Our predictions for the branching ratio given by 

Bi{D^ ^ K*-(+ve) = (2.63 ±0.81)% (HO), 

= (2.62 ±0.81)% (Linear), (7.31) 



are quite comparable to the available experimental data [49|: Broxp. = (2.02 ± 0.33)%. 
We summarize in Table |7.6| our results of the form factors V, Ai and A2 at both 
q^=0 and q'^=5'max compare with the available experimental data ||135|| as well 



as Eq. (|7.27|) and other theoretical results. While our form factors yli(0)=0.640 
(0.638) and 742(0)=0.450 (0.414) for the HO (linear) parameters are quite compara- 



ble with the data, Ai{0) = 0.58 ± 0.03 and A2(0)=0.41 ± 0.05 1T35|. our form factor 
y(0)=0.809 (0.822) underestimates the current data, 1/(0) = 1.07 ± 0.09 |^]. As 
one can see in Table [7^ , our analytic solutions for the form factors ^^(g^g^^) =1.029 
(1.024), Ai{ql^^)=0.698 (0.683) and A2(gLx)=0-534 (0.475) for the HO (linear) pa- 
rameters are well approximated not only by Eq. ( [71271) , 1^AiA2(g^g^)=1.029 (1.024), 
^iAiA2(9max) =0-698 (0.683) and A2A,A^iqi^J=0.53A (0.475), but also by the simple 
pole formular, VA,iqiJ=lM2 (1.033), Aia, (g^x) =0-697 (0.681) and A^^iqiJ = 
0.538(0.477). Our simple pole masses Ai=2.06 and 2.15 GeV (see Table |7l6|) for 
both HO and linear parameters, respectively, are also in good agreement with the 
value of D*^=2A1 GeV. We also present in Fig. 7.12 our dependence of the form 
factors V(g^), Ai(g^) and ^2(9^) and compare our Ai(g^) with the available lattice 
result [ p,10| | (data). Our results of Ailq"^) for both HO and linear parameter cases are 
in a good agreement with that of Ref. ||110| . 



(5) Dg r]{ri')i^i'£: These two semileptonic decays are very interesting processes to 
check our LFQM predictions of rj-i]' mixing angle. In our previous analysis of quark 
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potential model [|T3| in Chapter 4, we predicted the rj-T]' mixing angle for the HO 
(linear) potential model as 6'^^'^3-,=— 19.3°(— 19.6°). 

Our predicted decay rates for Dg —>■ 1] and Dg Vj in g^=0 frame are V(Ds — > 
1])= 0.100 (0.104)cos2(5|KsP ps"^ and Y{p, i)= 0.026 (0.028)sin2 pg-i, 
respectively, for the HO (linear) parameters. Using the lifetime rD^=0.467±0.017 ps, 
we obtain the branching ratio for 6su{2,)=—^'^° as follows:^ 

Br(D, ^ 7]tvi) = (1.7 ± 0.5)%, 

Bt{Ds r(tve) = (8.7 ±2.7) x 10"^ (HO), 

Br{Ds Titui) = (1.8 ± 0.6)%, 

Bt{Ds itvt) = (9.3 ± 2.9) X 10"^ (Linear). (7.32) 



While the experimental data are Bt^^pX^s ^ ri£^iy£) = (2.5±0.7)% and BicxpXDg 
7]'i^i'e) = {8.8±3A) X 10~^, respectively. Our results for the rj' decay mode are in a good 
agreement with the experimental data for both HO and linear parameters. For the 1] 
decay mode, our values are rather smaller than the central value of the experimental 
data, neverthless, are quite comparable with the data within the given uncertainties 
of CKM matrix, V^s- 

In Table |7]^, we compare our form factors f+''^{q'^) and f^"^ (g^) at both g^=0 and 
l^^lmax with Eq. (|7.27|) and the lattice calculation in Ref. ||106|| . Our predictions of 
/^-^(0)=0.692 (0.721)c^ and /+^'''(0)=0.692 (0.721)c^/ are in good agreement with 
the lattice results, /+=''(0) = (0.67± 0.07)c^ and /+^'''(0) = (0.67± 0.07)c^ |T06[, respec- 



tively. As one can see in Table [7l8| , our analytic solutions /+''^(5'max)= 1-31 (1.28)c^ 
and f^"^ (^7^)= 0.93 (0.94)c^/ for the HO (linear) parameters are well approximated 

by Eq. (^), ■^a(gLx) = l-27 (1.24)c, and -^i^^ (gLx)=0.92 (0.93)c,,, respec- 

^In numerical calculations, we use the common 77-77' mixing angle for both the HO and linear 
potentials. 
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Figure 7.13: The form factor f+iq"^) for the Dg ^ 7] transition with the 77-77' mixing 
angle ^5^/(3) =—19°. The same hne code as in Fig. 7.2 is used for our resuhs. 7.14: The 
form factor f+{q^) for the Ds — > rj' transition with the rj-r)' mixing angle 6*3(7(3)=— 19°. 
The same hne code as in Fig. 7.2 is used for our results. 
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Table 7.8: Summary of the parameters for the Ds viv') ^"^^ 
Note that c^=| cos5p| and Cr,' = | sin5p|. 



transitions. 







Ds^(j){ss) 


F 


u 


V 


Ai 


A2 


HO 


F{0) 


0.692c,, [0.692c,,'] 


0.830 


0.597 


0.409 




1.31c„ [0.93c,,'] 


1.070 


0.672 


0.492 


Ai [GeVJ 


1.94 [1.94] 


1.94 


2.83 


2.25 


S2A2 [GeVJ 


2.68 [2.68] 


2.69 


4.80 


3.00 


■^Ai (Qmax) 


1.49c,, [0.95c,,'] 


1.091 


0.673 


0.498 


^AiA^ ('Zmax) 


1.27c,, [0.92c,,'] 


1.069 


0.672 


0.491 


Linear 


m 


0.721c,, [0.721c,,'] 


0.861 


0.608 


0.385 




1.28c,, [0.94c,,'] 


1.080 


0.668 


0.445 


Ai [GeV] 


2.06 [2.06] 


2.05 


3.22 


2.53 


S2A2 [GeV] 


2.94 [2.94] 


2.94 


-12.91 


3.31 


^ Ai ('Zmax) 


1.37c„ [0.95c,,'] 


1.095 


0.666 


0.448 


^ AlA2 ('Zmax) 


1.24c,, [0.93c,,'J 


1.079 


0.666 


0.445 


LAT 106] 

|iog| 


m 
m 


0.67(7)c,, [0.67(7)c,,'] 
... [...] 


0.86±0.10 
1.30±0.32(43) 


0.52±0.03 
0.73±0.12(24) 


0.17±0.17 
0.55±0. 10(24) 



tively. From the resuhs of the simple pole approximation, jF;^^°^(g^)=1.49 (1.37)c,, 
and J^Xi^ (g^)=0.95 (0.95)c,,', one can find that Aj(z > 2) can be neglected for the 
Ds — > ^' transition but not for the Dg ^ rj case. For comparison, we also present in 
Table |7.9| the results obtained from 6su{3)=~^0° and —23° in g~'"=0 frame. Although 



it is not easy to conclude which mixing angle is the best from Table 7.9, we can see at 



least that our predictions with the mixing angle of ^5(7^3) ~~ 19° are overall in a good 
agreement with the available experimental data. The difference of the decay rates 
between the g"'"=0 and q'^^O frames is larger for the t] decay mode (19% for HO and 
18% for linear) than for the rj' decay mode (12% for HO and 11% for linear). Likewise, 



comparing the form factors /_ (0) for the two decay modes given by Table \7.1\ , one 
can see that the difference between the g"'"=0 and q~^^0 frames is much larger for the 
7] decay mode (80% for HO and 70% for linear) than for the 77' decay mode (4% for 
HO and 2% for linear). In Figs. 7.13 and 7.14, we present the dependence of the 
form factors f+{q^) for Dg —>■ 7] and Dg v( with the mixing angle 6'5[/(3)=— 19°, 
respectively. 



129 



Table 7.9: Branching ratio of the r],r]' transitions for various rj — rj' mixing 

angles, Osuj'i) = i^p + 35.26°. 



Processes 




^st/(3) — ^10° 


^5(7(3) - -19° 


(^su{3) - -23° 


Expt. 




HO 
Linear 


(2.5 ±0.8)% 
(2.6 ±0.8)% 


(1.7 ±0.5)% 
(1.8 ±0.6)% 


(1.4 ±0.4)% 
(1.5 ±0.5)% 


(2.5 ±0.7)% 




HO 
Linear 


(6.7 ±2.1) • IQ-^ 
(7.1 ±2.2) • 10-3 


(8.7 ±2.7) • 10-3 
(9.3 ±2.9) • 10-3 


(9.6 ±3.0) • 10-3 
(10.0 ±3.1) • 10-3 


(8.8 ±3.4) • 10-3 



(6) Dg — > (pi'^uf. Our predicted decay rate for Dg 0(ss) in q~^=0 frame is T{Ds — > 
(j)i~^i'e)= 4.99 (5.14)|VcsP X 10^^ ps-^ for the HO (linear) parameters. We obtain the 
ratio Tl/Tj' as 1.33 (1.37) for the HO (linear) parameters, while the lattice QCD 
calculation in Ref. [|106|] reported 1.49±0.19. Our predictions for the branching ratio 
are given by 



Br(D, ^ (Ptue) 



(2.51 ±0.77)% (HO), 
(2.59 ±0.80)% (Linear) 



(7.33) 



where we used 6v = —3.3° for u-cp mixing angle. Our results in Eq. ( |7.33|) are quite 



comparable to the available experimental data Brexp. = (2.0 ± 0.5)%, within the 
given uncertainties of the CKM matrix, Vcs- We summarize in Table ^]8| our re- 
sults of the form factors V, A\ and at both g^=0 and q'^=q'max ^^'^ compare with 
Eq. ( |7.27| ) and the lattice results | |106| , [109| ] . Our results of the form factors V , and 
A2 at g^=0 are quite comparable with the lattice calculations in Refs. [|106| , |109|| . As 
one can see in Table |778| , our analytic solutions for the form factors l^(g^3^^) =1.070 
(1.080), Ai(g^^^)=0.672 (0.666) and A2(gLx) =0-492 (0.445) for the HO (linear) pa- 
rameters are well approximated not only by Eq. (|7.27|) , VA.iA2(f7max)=l-069 (1.079), 
^iAiA2(gmax) =0-672 (0.666) and A2AiA2(gLx)=0-491 (0.445), but also by the simple 
pole formular, \/A,(gLx)=l-091 (1-095), Aia^ (gLx) =0-673 (0.666) and A2Ai(gLx) = 
0.498(0.448). The simple pole approximation is already enough good approximation 



for 



transition. We also present in Fig. 7.15 our dependence of the form 
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Figure 7.15: The form factors V^(g^), ^1(5'^) and A2{q^) for 
The same hne code as in Fig. 7.11 is used. 

factors V{q'^), Ai{q'^) and ^2(9^). 



[ss] transition. 



7.4.3 B and Bg decays 



(1) ^ -K i^uf. Our predicted decay rate for i? ^ tt in g+ = frame is r(i?° — >■ 
TT^i^ug) = 7.06(8.16)1 Vafep ps~^ for the HO (hnear) parameters. Using the hfetime 
tbo = 1.56 ± 0.04 ps and |Kfe| = (3.3 ± 0.4 ± 0.7) x lO'^O, our predictions for the 
branching ratio given by 

Br(5° ^ n-e-^ue) = (1.20 ± 0.29) x 10"^ (HO), 

= (1.40 ± 0.34) X 10"^ (Linear), (7.34) 

are quite comparable wiht the recent experimental data ||138|| , Br(i? 7r^£+z/^) = (1.8± 
0.6) X 10~^, within the given error range. We also performed the calculations of the 
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Table 7.10: Summary of the parameters for the B 



APE Wm and ELC 



vr and B p transitions. The 
results are from their method 'b', which uses the heavy 



quark scaling laws to extrapolate from D- to 5- mesons at fixed to. 





B ^71 


B^p 


F 


u 


V 


Ai 


A2 


HO 


m 


0.234 


0.273 


0.216 


0.196 


-^(O'max) 


3.34 


0.703 


0.449 


0.857 


Ai [GeV] 


4.37 


4.35 


6.91 


4.70 


S2A2 [GeV] 


4.74 


5.84 


11.01 


6.40 


^ A.1 (^max) 


-0.61 


-3.66 


0.376 


2.43 


^ A.iA2(?max) 


1.60 


0.972 


0.359 


0.599 


Linear 


m 


0.273 


0.324 


0.249 


0.220 


'^(^max) 


2.80 


0.688 


0.461 


0.851 


Ai [GeV] 


4.59 


4.57 


8.06 


5.09 


S2A2 [GeV] 


6.53 


7.11 


14.62 


7.02 


A_i (O'max) 


-1.07 


13.00 


0.362 


1.026 


^ A.iA2(Q'max) 


2.16 


1.736 


0.358 


0.573 


LAT |104 




F(0) 

'^(^'max) 


0.27±0.11 


qpi+U-(Jti 

U-'JO-O.OS 


n 97+U.U6 

U--^' -0.04 


0.26:^^^ 

0.88lHi 




105 
111 


m 
m 


0.35±0.08 
0.30±0.14(5) 


0.53±0.31 
0.37±0.11 


0.24±0.12 
0.22±0.05 


0.27±0.80 
0.49±0.21(5) 


SR 


112 
113 


] 
J 


m 
m 




0.34±0.05 
0.35±0.07 


0.26±0.04 
0.27±0.05 


0.22±0.03 
0.28±0.05 



decay rate using the valence contributions from q^^O frame and found that the decay 
rate obtained from q~^y^O frame is about 9 (7)% smaller for the HO (linear) parameters 
than the result obtained from q~^=0 frame. 



We summarize in Table [7.10| our analytic solutions of f^^ obtained at both g^=0 
and 5'^=5'max compare with Eq. (|7.27| ) as well as other theoretical results. As 
one can see in Table 7.10 , our results of f^'"{0) for both HO and linear parameters 
are quite comparable with the lattice QCD calculations ||104| , |105| , |1 1 1|| . Gomparing 
with the pole dominace formular given by Eq. (|7.27|) , we found that our solution of 
/+''(Q'max)=3.34 (2.80) for the HO (linear) parameters are neither fitted by Eq. ( |7.27|) , 



•^AiA2(?max) = l-60 (2.16), uor by the simple pole approximation, J^Xi{qma.x)=~^-^^ 
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Table 7.11: Comparison of the form factors f-{q^) at g^=0 between g"'"=0 and 
frames. The differences between the two frames are the measure of the nonvalence 
contributions in g^T^O frame. 



Models 




B ^71 


B ^ D 


Bs^D, 


HO 


g+=0 


-0.203 
-0.235 


-0.318 
-0.478 


-0.305 
-0.518 


Linear 


g+=0 


-0.238 
-0.273 


-0.328 
-0.499 


-0.325 
-0.533 



(—1.07). Even though we do not know yet how many terms of Ai{i > 3) are needed 
to obtain the well approximated solutions motivated by the pole dominace model, it 
seems quite clear that the pole approximation up to A2 given by Eq. ( [7.27| ) does not 



work well for the heavy-to-light transitions. In Table |7.11| , we also summarize the 
results of the form factor /-(O) at q^=0 obtained from both q'^=0 and q^^O frames. 
The difference, i.e., nonvalence contribution from q^y^O frame, of f^'"{0) between the 
two frames is about 16 (15)% for the HO (linear) parameters. 

In Figs. 7.16(a) and 7.16(b), we present the dependence of the form factors f^'^{q^) 
and fo^iq"^) obtained from both g+=0 (thick lines) and q^^O (thin lines) frames, re- 



spectively, and compare with the available lattice QCD calculations ||108|| . Again, the 
nonvalence contributions from q~^y^O frame are quite sizeable especially for f^^{q'^) 
case. Our analytic solutions of f^^{q^) for both HO and linear parameters are overall 



in a good agreement with the lattice calculations in Ref. |108|. 



(2) B^ —>■ p^i^uf Our predicted decay rate for i? ^ p in g+=0 frame is r{B^ 
p~i^i'i)= 11.44 (14.25)|Vufep ps~^ for the HO (linear) parameters. We obtain the ratio 
^l/^t as 1.07 (1.19) for the HO (linear) parameters, while other QM calculations 
predicted as 0.3 |81[, 1.34 ITT8[ and 1.13 |123[. The lattice QCD of Ref. and 



the QCD sum rule of Ref. ||113|| also calculated the ratio as 0.80lo;o3 and 0.52±0.08 
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Figure 7.17: (a) The form factor Vlq"^) for the ^ p transition. The thick (thin) 
sohd hne corresponds to our HO (hnear) parameters. For comparison, we include 
the resuhs of J^K^K.^i.'f) (dashed hne) given by Eq. ( |7.271 ) and lattice QCD from 
UKQCD (•), APE [|mg (square) and ELC fTTTI (o). (b) The form factor Ai{q^) 
for the ^ p transition. The same line code as in (a) is used. 
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1.6 




Figure 7.17: (c) The form factor A2{q'^) for the ^ p transition. The same hne 
code as in (a) is used. 

respectively. Our predictions for the branching ratio are given by 



:i.94l}:^^) X 10-^ (HO), 



i2.36tlit) X 10-^ (Linear) 



(7.35) 



While the central value of linear parameters is quite comparable to the available 
experimental data Brexp. = (2.5^5 o) ^ 10~^, the central value of the HO parameters 



underestimates the data. 



We summarize in Table [7?T0| our results of the form factors V, Ai and A2 at both g^=0 
and 5'^=5'max compare with Eq. ( [7.27|) and other theoretical results. Our results 
of the form factors V, Ai and A2 at g^=0 are in good agreement with the lattice 
calculations in Refs. ||104| , |105| , |lllj| as well as the QCD sum rule in Refs. ||112| , |113 |. 
Moreover, our results of the axial vector form factors Ai{q'^^^) =0.44:9 (0.461) and 
^2(Q'max)=0-857 (0.851) for the HO (linear) parameters are in excellent agreement 
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with the lattice results in Ref. [|T0|]: Ai(g^^^)=0.46t[}:^? and A2(gLx)=0-88+°-°^ 



0.035 

2 



respectively. Our prediction of the vector form factor V{ql^^^)=0.703 (0.688) for 
the HO (linear) parameters, however, disagree with the lattice result in Ref. ||104|| , 
^(?max)=2.07to;o6- ouc cau see in Table [7.10| , our analytic solutions for the form 
factors V(q'max) ^i^nd A2{q'^^^) are neither approximated by Eq. ( [7.27|) nor the simple 
pole formular, i.e., ^Ai(5'max)- Only the form factor of Ai is reasonably approximated 
(but still ~ 20% difference from our analytic solution) by the simple pole dominace, 
i.e., AiAi (gmax)=0-376 (0.362) for the HO (linear) parameters with the pole mass Ai = 
6.91 (8.06) GeV. It is interesting to note that UKQCD ||103|| in their analysis showed 



the simple pole behavior for Ai(g^) is preferred with mpoie = 7^^ GeV, which is quite 
compatible with our values of Ai. 

In Figs. 7.17(a), 17(b) and 17(c), we show the dependence of the form factors V, 
Ai and A2, respectively, and compare with the available lattice results ||103 , 105 , |111|| . 



For comparison, we also include in Fig. 7.17(a) the result (dashed line) of Eq. (|7.27|) , 
i.e., ^AiA2(q'^) obtained from the linear parameters. As one can see in Figs. 7.17(a) 
and 17(b), our analytic solutions for the Ailq"^) and A2{q'^) are in excellent agreement 
with the lattice calculations ||103| , |105| , |111|| . For the vector form factor V{q^) case. 



our analytic solution gives reasonable agreement with the ELC result ||111|| (o) but 
underestimates the result of UKQCD ||103|| (•). 



(3) B{Bs) D{Ds)£i'e: Our predicted decay rates for B ^ D and Bg —>■ Dg in 
q+=0 frame are 1(5° ^ D-i+ue)= 8.93 (9.39)|KfeP ps"^ and 1(5, ^ D;i+ue)= 8.80 
(9.30)|Vcfep ps~^, respectively, for the HO (linear) parameters. Using rB^ = (1.54±0.07) 
ps and I Vbc|=0.0395±0.003 [HUl, we obtain the branching ratio as follows: 



Br(B° D-tui) = (2.17 ± 0.19)%, 
Bt{Bs D-tvt) = (2.11 ± 0.18)% (HO), 
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Table 7.12: Summary of the parameters for the B D{Bs —>■ Dg) and B —>■ D*{Bs — > 
D*) transitions. The experimental data for B ^ D* transition was obtained by as- 
suming an exponential dependence of the form factors on according to the formal- 
ism outline in Ref. 





B -> D{Bs ^ D,) 


B D*{Bs Dl) 


p 


J+ 


T/ 
V 


^1 


■^2 


HO 


m 


0.686 (0.664) 


0.711 (0.692) 


0.652 (0.619) 


0.581 (0.541) 




1.11 (1.13) 


1.10 (1.11) 


0.80 (0.80) 


0.86 (0.84) 


Ai [GeV] 


5.22 (4.95) 


5.25 (4.94) 


7.74 (6.77) 


5.49 (5.15) 


S2A2 [GcV] 


7.35 (6.64) 


7.45 (6.62) 


-19.17 (10.07) 


7.81 (6.88) 


^ Ai (<Zmax) 


1.20 (1.26) 


1.16 (1.22) 


0.79 (0.81) 


0.90 (0.90) 


^ A.1A2 ilmax) 


1.11 (1.11) 


1.10 (1.11) 


0.79 (0.79) 


0.86 (0.83) 


Linear 


m 


0.709 (0.689) 


0.727 (0.717) 


0.667 (0.638) 


0.580 (0.549) 




1.12 (1.14) 


1.10 (1.13) 


0.80 (0.80) 


0.83 (0.82) 


Ai [GeV] 


5.38 (5.09) 


5.39 (5.08) 


8.24 (7.16) 


5.70 (5.35) 


S2A2 [GeV] 


7.77 (6.93) 


7.82 (6.91) 


-13.68 (11.46) 


8.29 (7.24) 


^ Ai (^max) 


1.19 (1.24) 


1.15 (1.22) 


0.79 (0.80) 


0.86 (0.87) 


^ A.1A2 (9max) 


1.12 (1.13) 


1.10 (1.12) 


0.79 (0.80) 


0.83 (0.82) 


ISGW |BC 
BSW |118 


1 


"^(^max) 

F{ql,x} 




1.20 (•••) 
0.97 (• • •) 


0.94 (•••) 
0.85 (•••) 


1.08 (•••) 
0.90 (•••) 


EXP. |148 




^(9max) 




0.91±0.49±0.12 


0.85±0.07±0.11 


0.87±0.22±0.10 



Br(S° ^ D-tu^) = (2.28 ± 0.20)%, 

Br(5, ^ D-tue) = (2.23 ± 0.20)% (Linear). 



(7.36) 



Our predictions of the B ^ D transition for both HO and linear cases are in a good 
agreement with the experimental data [|9|, Br(5° L>-£+z/f) = (2.00±0.25)%. The 
decay rates of B ^ D and Bs Dg obtained from the valence contributions in 
7^ frame for the HO (linear) parameters are about 6 (6)% and 8 (8)% larger 
than those obtained from g"''=0 frame, respectively. 



We summarize in Table |7.12| our analytic solutions of f^^ and f:^"^" obtained at 
both g^=0 and q'^=q'^ax compare with Eq. ( [7.27|) . As one can see in Table [7.12| , 
our results of /:f^(g^^J=l.ll (1.12) and =^=(g^^^)=1.13 (1.14) for the HO (lin- 
ear) parameters are well approximated by Eq. (|7.27| ), jF^^^^(g^^^)=l.ll (1.12) and 



^AiM ilmsix) — ^-^^ (1.13). The results of the monopole approximations, ^yfi^(i?max) = l-20 
(1.19) and jF£^^'^°(g^^^)=1.26 (1.24) are also quite comparable to our analytic solu- 
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tions. The form factors f^^{0) and /^"^"{O) obtained from g+=0 frame are also 
summarized in Table [7.11| and compared with those obtained from q~^y^O frame. As 
one can see m Table the results of /:^^(0) and f^'^' (0) in frame are about 
50 (52)% and 70 (64)% larger for the HO (linear) parameters than those in q~^=0 
frame, respectively. In Figs. 7.18 and 7.19, we show the dependence behaviors of 
the form factors /+ and /q for B ^ D and Bs transitions, respectively. For 

comparison, we also include the results of the valence contributions in g^T^O frame. 
For these heavy-to-heavy transition cases, we can easily see the nonvalence contribu- 
tions (i.e., the difference of the results between the two frames) are much suppressed 
compared to the previous light-to-light and heavy-to-light transition cases. 

These heavy-to-heavy transitions are also used to investigate model reliability by 
checking the universal IW function given by Eq. (^3) heavy-quark symmetry. The 
slope of the IW function at the zero-recoil point is defined as 

^{w) = l-p\w-l\ (7.37) 

where w = vi ■ V2 = (Mf + M| — g^)/(2MiM2). Our predictions of the slope = 
0.80 < < 0.92 for possible combination^ of (3d and (3b are quite comparable with 
the current world average p^^g =0.66±0.19 extracted from exclusive semileptonic 
B D£u decay as well as other theoretical estimates, 0.7 < p^ < 0.88 in 



and 0.87 in ||139|| . In Fig. 7.20, our predictions of the IW function ^^^{vi ■ V2) for 
the HO (sohd) with the common gaussian /3 = /3b=0.496 and linear (dotted line) 
parameters with the common (3 = /?s=0.5266, respectively, are compared with the 



available experimental data ||140| , |141|| and show good agreement with the data. Other 
combinations of the gaussian (3 paramters give very similar behavior to those as 



^ The slopes of the IW function for the HO (Unear) parameters are p^=0.85 (0.80) for (3b = 
/3r,=0.496 (0.5266) GeV, 0.92 (0.86) for /3b = /3d=0.4216 (0.4679) GeV and 0.90 (0.85) for /?b 7^ /3d 
case as usual, respectively. 
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Figure 7.18: (a) The form factor f+{q^) for B^D transition. The same hne code 
as in Fig. 7.2 is used, (b) The form factor /o(g^) for B^D transition. The same hne 
code as in Fig. 7.2 is used. 
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Figure 7.19: (a) The form factor f+{q^) for Bg^Dg transition. The same hne code 
as in Fig. 7.2 is used, (b) The form factor fo{q^) for Bg^Dg transition. The same 
hne code as in Fig. 7.2 is used. 
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Figure 7.20: The IW function ^^^(fi • ^2) for the HO (sohd) with the com- 
mon gaussian f] = /3b =0.496 and hnear (dotted hne) parameters with the common 
P = /?_B=0.5266, respectively. We compare our results with the experimental data of 
ARGUS [|Tg (square) and CLEO (circle). 



shown in Fig. 7.20. 



(2) S°(S°) ^ D*+{D*+y-vf. Our predicted decay rates for B ^ D* and B., Ds in 
g+=0 frame are 1(5° D*+(i-Ui)= 22.13 (22.90)|VbcP ps^^ and 1(5° D*+i~u^) = 
21.34 (22.25)|VfecP ps~^ for the HO (linear) parameters, respectively. The ratio rj;,/rT 
foiB ^ D and 5, D, decays are obtained as 1.20 (1.23) and 1.19 (1.21) for the HO 
(linear) parameters, respectively. While the experimental data for B D* reported 
as 0.85±0.45 [ig and 0.82±0.36 |43 . 



Our predictions for the branching ratio are given by 



Br(5° ^ D*+rue) = (5.39 ± 0.8 



Br(B° ^ D*+rue) = (5.13 ± 0.75)% (HO), 
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Figure 7.21: The dV/dq^ distribution for B D*ii' decays using the hnear 
parameters. The thick sohd, dotted, and long-dashed hues are the results from 
the central value, upper limit, and the lower limit of the CKM matrix element. 



I Vbc|=0.0395±0.003, respectively. The experimental data are taken from Ref. [|147 



Br(5° ^ D*+rue) = (5.58 ± 0.85)%, 

Br(5° ^ D*+e~ue) = (5.35 ± 0.78)% (Linear). (7.38) 

Even though the current average value for —>■ D* decays, Brexp.(-B'^ —>■ D*^i^i?i) 
= (4.60±0.27)%, is smaller than our predictions, it is interesting to note that the recent 
experiments ||144| , |145| , |146|| reported results close to our predictions, i.e., (5.08±0.21± 



0.66)% m, (5.53±0.26±0.52)% [|l45l, and (5.52±0.17±0.68)% 1T4|, respectively 



In Fig. 7.21, we present the differential rate for B D* decays using the linear 
parameters and compare with the available experimental data ||147|| . Within the 
uncertainties of the CKM matrix element | Vf,c|=0.0395±0.003, our results are overall 
in good agreement with the data. The results for the HO parameters are not much 
different from those for the linear parameters. 
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We summarize in Table |7.12| our analytic solutions of V, Ai and A2 obtained at both 
g^=0 and q^=q^^^ and compare with Eq. (|7.27|) as well as the available experimental 



data | 148 | and other theoretical results |^U|, 118 |. As one can see in Table [7.12| , our 
results of the three form factors for B ^ D* transition, i.e., V{q^^^=l.lQ (1-10), 
Ai(g^^^)=0.80 (0.80) and A2(gLx)=0-86 (0.83) for the HO (linear) parameters are 
in good agreement with the experimental data [ p.48|| , \4xp.(?max)=0-91='=0.49±0.12, 
v4iexp.(gLx)=0-85±0.07±0.11 and A2exp.(9max)=0-87±0.22±0.10, respectively. Our 



analytic solutions are also well approximated by not only Eq. ( |7.27| ) but also the 
monopole form factor JF^^ (see Table |7.12| ). In Figs. 7.22 and 7.23, we show the 
dependence behaviors of the three form factors for B D* and B^ D* decays, 
respectively. 



Experimentally, two form-factor ratios for B ^ D* decays defined by ||132| , [149 



Ri{q') 
R2{q') 



1 - 



q 

q' 



A,{q^y 
Aiiq^y 



(7.39) 



have been measured by CLEO [|14|] as follows: i?i(gLx) = 1-24 ± 0.26 ± 0.12 
and i?2(gLx) = 0-72 ± 0.18 ± 0.07. From Table ^1% we obtain the ratios as 
/?i(g^^^) = 1.10 (1.11) and /2i(g^^x)=0.86 (0.84) for the HO (linear) parameters, which 
are in good agreement with the data. Our results are also agree with the predictions 
of Neubert [|13|], i?i(g^^x)=1.35 and i?2(gLx)=0-79, and the ISGW2 model [|ll, 
i?i(g^j^^) = 1.27 and i?2(5'max)=l-01- ■'■^ -^^S- 7.24, we present the two form-factor ra- 
tios Ri and R2 as a function of w and compare with the data ||149|| at q^ = ^^^^^(or 
w=l) anticipating more data in the entire physical q^ range. Our results show almost 
w independent behavior as indicated by Ref. ||132|| . The form factors for 0^ 1^ 
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Figure 7.22: The form factors V^(g^), and ^2(5^) for B ^ D* transition. 

The same hne code as in Fig. 7.11 is used. 7.23: The form factors V{q'^), ^i(g^) and 
A2{q'^) for B ^ D* transition. The same hne code as in Fig. 7.11 is used. 
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Figure 7.24: The two form-factor ratios, Ri{w) (solid line) and R2{w) (dotted line) 
for the HO parameters compared with the experimental data ||149|| of -Ri(l) (circle) 
and -R2(l) (square). 7.25: Our predictions of the ratio of C,^^{w)/^^^* (w) for both 
HO (solid line) and linear (dotted line) parameters compared with the experimental 
data m. 
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transitions in the heavy-quark hmit are also related to the IW function via 

Mp + Mn[l- qV{Mp + Mnf] ' ^ ' ' 

The slope of the IW function at zero-recoil for the B —>■ D* decays is obtained as 
0.85 < < 0.97 for all possible combinations of the gaussian parameters jSp and jSo*- 
Our results are in good agreement with the data: p^=0.91±0.15±0.06 (CLEO) ||149 



and p2^g =0.71±0.11 (PDG) ||. We present in Fig. 7.25 the ratio of the two IW 

functions for B ^ D andB ^ D*, i.e., ^^^{w)/^^^' (w), where ^^^{w) and 

is obtained by the form factors f+{q^) and V{q^)^ respectively, and compare with the 



data [^. Our results for both HO (solid line) and linear (dotted line) parameters 
are almost equal to 1, the value expected from the heavy quark limit (m;,, rric oo) 
and thus show a good agreement with the data |^ . 



Moreover, the decay B^ D*~^i~ue allows one to measure the chirality of the weak 
b c transition | |132| , |140| , |149| , |15U| | , which is only possible if the daughter meson has 
spin J > 0. This b c chirality can be tested by measuring the forward-backward 
asymmetry [|13|, |148|, 



^ 3 JKfiq^)q^i\H_\'-\H_\')dq^dpe 

4JKf{q^)q^{\H+\'^ + \H_\^ + \Ho\^)dq^dpe' 

3r_ -r. 



4 r 

The helicity alignment of the iy~ is given by ||140 



(7.41) 



^poi = 2 ^^^°^ -1, (7.42) 

which describes the D*"*" polarization extracted from the D*~^ decay angle distribu- 
tion (see for more detail in Ref. ||132| , |150|| ). From Eqs. (|7.41|) and (p.lO|) , we ob- 



tain Af^*= 0.195 (0.193) and A^^*= 1.396 (1.510) for the HO (linear) parameters. 



147 



Table 7.13: Decay rate F (in unit of ps ^) and branching ratio Br(0 —>■ iu^) 
for various semileptonic decays. We use ^'^^'^3-)=— 19° to obtain the branching ra- 
tio for Ds^i]{ri') decays. The used CKM matrix elements are | K;s|=1.04±0.16, 
|Kd|=0.224±0.016, |Kf,| = (3.3±0.4±0.7)xlO-3, and |Hc|=0.0395±0.003 ^ 



Mode 




r 


Br. 


Expt.(Br.) 


D ^71 


HO 


O.llOlKdP 


(2.30 ±0.33) • 10-3 


(3.7 ±0.6) • 10-3 




Linear 


0.113|KdP 


(2.36 ±0.34) ■ 10-3 




D ^ K 


HO 


8.26-10~2|y^^|2 


(3.71±1.14)%[e+z/e] 


(3.66 ±0.18)% 






6.40-10-2 2 


(2.87±0.88)%[/i+z/^] 


for ^ K-e+Ue 




Linear 


8.36-10-2|V;,|2 


(3.75± 1.16)%[e+z/e] 


(3.23 ±0.17)% 






6.43-10-2|V;,|2 


(2.89±0.90)%[/i+z/^] 


for L)0 ^ K-fx+Uf, 




HO 


O.lOOcos^^lK^p 
0.104cos2 5 


(1.7 ±0.5)% 


(2.5 ±0.7)% 




Linear 


(1.8 ±0.6)% 




Ds^i 


HO 


0.026sin^5|Ks|^ 
0.028sin2 5 


(8.7 ±2.7) ■ 10-3 


(8.8 ±3.4) • 10-3 




Linear 


(9.3 ±2.9) ■ 10-3 




B 


HO 


7.06|K6|' 


(1.20 ±0.29) ■ 10-^ 


(1.8 ±0.6) ■ 10-^ 




Linear 


8.16|Kb|' 


(1.40 ±0.34) ■ 10-^ 




B^ D 


HO 
Linear 


8.93|KfeP 
9.39|KfeP 


(2.17 ±0.19)% 
(2.28 ± 0.20)% 


(2.00 ± 0.25)% 


Bs~^D, 


HO 


8.80|Kfe|' 


(2.11 ±0.18)% 






Linear 


9.30|Kfe|' 


(2.23 ±0.20)% 





Our results for B ^ D* are in excellent agreement with the available experimental 
data: Afb=0.20±0.08±0.06, Apoi=l.l±0.4±0.2 HO] and Afb=0.197±0.033±0.016, 



/lpoi=1.55±0.26±0.13 |14|. Similarly, we also obtain Ap^^''= 0.194 (0.194) and 



A 



BsD', 
pol 



1.376 (1.422) for the HO (linear) parameters anticipating future experimen- 



tal data for Bg ^ D*. Finally, we summarize in Table |7.13| and 7.14 our results of 
the decay rates for 0- 0^ and 0- — > 1- transitions, respectively. 



7.5 Summary and Discussion 



In conclusion, in this Chapter, we analyzed the exclusive 0-— >0- and 0-^1- semilep- 
tonic heavy meson decays extending our LFQM constrained by the variational prin- 
ciple for the QCD-motivated effective Hamiltonian disscussed in Chapter 4. Our 
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Table 7.14: Decay rate F (in unit of jV^igaP ps~^) and branching ratio Br(0^ —>■ l^iue) 
for various semileptonic decays. We use |5^/|=3.3° to obtain the branching ratio for 
Dg (f) decays. The Tl and Tt are the decay rates for the longitudinal and transverse 
helicity contributions, respectively. 



Mode 




r 




Br. 


Expt.(Br.) 


K* 


HO 
Linear 


5.856-10-2 
5.846-10-2 


1.37 
1.42 


(2.63 ±0.81)% 
(2.62 ±0.81)% 


(2.02 ±0.33)% 




HO 
Linear 


6.382-10-2 
6.610-10-2 


1.39 
1.47 


(1.33 ±0.19) - 10-=^ 
(1.38 ±0.19) - 10-3 


(1.87 ±0.9) - 10-=^ 




HO 
Linear 


4.994-10-2 
5.137-10-2 


1.33 
1.37 


(2.51 ±0.77)% 
(2.59 ±0.80)% 


(2.0 ±0.5)% 


B^-^p- 


HO 
Linear 


11.445 
14.254 


1.07 
1.19 


(2.361^1^) - 10-^ 


(2.51^:^) • 10^^ 


5° D* 


HO 

Linear 


22.133 
22.904 


1.20 
1.23 


(5.39 ±0.82)% 
(5.58 ±0.85)% 


(4.60 ±0.27)% 


B's ^ d: 


HO 
Linear 


21.341 
22.253 


1.19 
1.21 


(5.13 ±0.75)% 
(5.35 ±0.78)% 





model not only provided overall a good agreement with the available experimental 
data and the lattice QCD results for the weak transition form factors and branching 
ratios of the light-to-light {K^^), heavy-to-light and heavy-to-heavy meson decays but 
also rendered a large number of predictions to the heavy meson mass spectra and 
decay constants. Our predicted meson mass spectra given by Table fO were overall 



in a good agreement with the data [^. Our values of the decay constants given by 
Table were also in a good agreement with the results of lattice QCD ||102| , |107|| 
anticipating future accurate experimental data. 



As an application of our model, we first calculated the spacelike EM form factors of D 
and B mesons as well as vr and K mesons in q~^=0 frame and estimated the nonvalence 
contributions of the EM form factors from the frame by comparing the valence 

contributions obtained from the q~^^0 frame with those obtained from the q'^=0 
frame. The nonvalence contributions from the q'^^0 frame are highly suppressed as 
the quark mass increases, i.e., in the order ofTT>K>D>B,a.s one can see from 
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Figs. 7.2-7.5. 

We have overcome the difficulty associated with the nonvalence Z-graph contribution 
in timehke region by the analytic continuation of weak form factors from the spacelike 
region. In other words, the form factors f± for 0^ 0" decays and the form factors 
V, Ai and A2 were obtained in the g"'"=0 frame and then analytically continued to 
the timelike region by changing qj_ to iq± in the form factors. The matrix element of 
the " _L " component of the current was used to obtain the form factor /_, since 
the J~ is not immune to the zero mode contributions even in the q~^=0 frame 



Our numerical computation conffimed the equivalence of our analytic continuation 
method and the dispersion relation method ||123|| p|. The nonvalence contributions 
from the q~^^0 were also quantified for the 0~ —>■ 0~ decays by comparing the valence 
contributions between the q~^=0 and q^y^O frames. Another interesting thing is that 
our predicted value of the rj-rj' mixing angle, i.e., Osu(3) ~ —19°, works fairly well for 
the predicitions of the decay rates for Ds viv') processes. 

Furthermore, our analytic solutions of the form factors for both 0~ 0~ and 0^ 1~ 
decays are well approximated by monopole-type form factors except heavy-to-light 
processes such as B ^ n and B ^ p. For these two processes, our analytic solutions 
of the form factors are neither approximated by monopole-type form factors nor the 
Eq. (|7.27| ). In order to explain the discrepancies between our analytic solutions and 
the approximations of the pole dominace model, more accurate experimental data 
or the lattice results should be accomodated. Nonetheless, if we assume the pole 
dominance model is a good approximation for the heavy-to-light decays, this may 
indicate that our LQFM with a simple picture of qq contents for the description of 



^If we were to use the model parameters given in ]32| , we obtain, for example, the values of 
/+(0) - /+(9max) as follows: 0.783 (0.781) - 1.2 (1.2) for D ^ K, 0.682 (0.681) - 1.61 (1.63) for 
D ^TT, 0.682 (0.684) - 1.12 (1.12) for B D, and 0.293 (0.293) - 2.7 (2.3) for B tt, where the 
values in parentheses are the results obtained by the auther in Il23| . 
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a meson is not valid for really high (f region. Anyway, the comparison between our 

analytic solutions and the results of the pole dominance model provides a good testing 
ground for the validity of our model calculations. 
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Chapter 8 
Conclusion 



In this thesis, we have investigated the electroweak form factors and the semileptonic 
decays of pseudoscalar and vector mesons within the framework of hght-front con- 
stituent quark model. In Chapters 2 and 3 of this thesis, we have discussed the two 
different schemes in treating the meson masses, e.g., SM and IM schemes. Regardless 
of the difference between the two schemes, once the best fit parameters were used, 
both schemes provided the predictions that were not only pretty similar with each 
other but also remarkably good compared to the available experimental data for form 
factors, decay constants, chargi radii, etc., of various light pseudoscalar and vector 
mesons as well as their radiative decay widths. Similarly, once the best fit parameters 
were chosen, the difference from the Jacobi factor in the IM model wave functions is 
substantially reduced in the numerical predictions for the physical observables. How- 
ever, the difference in the choice of the radial wave function, e.g., harmonic-oscillator 
(HO) wave function versus power-law (PL) wave function, was appreciable no matter 
what parameters were used. For example, in the phenomenology of various meson 
radiative decays at low Q^, we observed that the Gaussian type wave function was 
clearly better than the PL wave function in comparison with the available experimen- 
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tal data. 

In order to justify the model wave function as a solution of the QCD-inspired dy- 
namic equation, we applied the variational principle in Chapter 4 of this thesis to 
the QCD-motivated potential, which includes not only the Coulomb plus confining 
potential but also the hyperfine interaction, to obtain the correct p-vr splitting. For 
the confining potential, we took both HO (~ r^) and linear (~ r) type potentials 
and compared the numerical results for these two cases. The variational principle 
for the effective Hamiltonian turns out to be crucial to find the optimum values of 
our model parameters as shown in Chapters 4 and 7. We adopted the IM scheme 
to assure the orthogonality of model wave functions. As shown in Figs. 4.1(a) and 
4.1(b), our central potentials for both HO and linear potentials in Eq. (^]^) are not 
only very similar to each other but also quite close to the Isgur-Scora-Grinstein-Wise 
model 2(ISGW2) ||81|] potentials. Using our quark potential model discussed in Chap- 
ter 4, we predicted the ground state meson spectra (see Table [773|) . Our predictions 
of ground state meson mass spectra agreed quite well with the experimental data |^ 
(within 6% error). Furthermore, our model predicted the two unmeasured mass spec- 
tra of ^So{bb) and ^Si{bs) systems as 71/^^=9295 (9657) MeV and Mb5=5471 (5424) 
MeV for the HO (linear) potential, respectively. 

We have also predicted the uj-(j) and 7]-r]' mixing angles using the parametrization to 
incorporate the quark annihilation diagrams P7| , [75| , [75[ | mediated by gluon exchanges 



and the SU(3) symmetry breaking, i.e., m„((i) ^ nis. Our results of the uj-(j) and rj-i]' 
mixing angles are \6v\ ~ 4.2°(7.8°) and Osuis) ~ -19.3°(-19.6°) for the HO (linear) 
potential model, respectively. The sensitivity of the {uj,(f)) mass spectra for ~ 1°(5°) 
variation of 6v, i.e., from Sy = 4.2°(7.8°) to 3.3° for the HO (linear) potential case, is 
within the 1%(5%) level. As we discussed in the appendix E, the sign of mixing angle 



153 



depends on the sign of SU(3) breaking parameter X (see Eqs. (E.3)-(E.5)). While 
Xp > is well supported by the Particle Data Group [||] (-23° < ^^^^3) < -10°), 
the sign of Xy is not yet definite at the present stage of phenomenology. Regarding 
the sign of Xy, it is interesting to note that 5v ~ — 3.3°(= ^^5(7(3) ~ 35.26°) (i.e., 
Xy < 0) is favored in Refs. [l^, ^ |8^, while the conventional Gell-Mann-Okubo 
mass formula for the exact SU(3) limit (X 1) predicts 6y ~ 0° in the linear mass 



scheme and 6y ~ +3.3° (i.e., Xy > 0) in the quadratic mass scheme Even 
though it is not yet clear which sign of uj-(j) mixing angle should be taken, the overall 
agreement between our HO potential model with the positive sign, i.e., 6y ~ +3.3°, 
and the available experimental data seem to be quite good. If we were to choose the 
sign of X as X > in Eq. (E.4), then the fact that the mass difference — rup 
is positive is correlated with the sign of the uj-(j) mixing angle [Q. In other words, 
rriuj > mp implies 6v > from Eqs. (E.3)-(E.5). Perhaps, the precision measurement 
of — s> ?7'7 envisioned in the future at TJNAF experiments might be helpful to 
give a more stringent test of 5y. We also predicted the decay constants for various 
pseudoscalar and vector mesons as summarized in Tables [4.3| and [r.4| . Our predictions 
for the decay constants are quite consistent with the available experimental data |^ 



as well as the lattice QCD ||102| , |107|| . For the heavy meson decay constants, however. 



the current experimental data have large error-bars and the more accurate data are 
needed. 

We further applied our quark potential model to compute the radiative (see Chapter 4) 
and the semileptonic (see Chapter?) decays of pseudoscalar and vector mesons. One of 
the distinguished features in our LF approach is the ability to compute the timelike 
form factors as shown in 0^ 0^ and 0^ 1~ semileptonic decays. While we 
calculated the exclusive processes in the Drell-Yan-West {q~^=0) frame, we analytically 
continued the amplitudes to the timelike region. Our analytic continuation method 
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(i.e. changing q±_ to iq\_ in the form factors) avoided encountering the nonvalence 
diagram (black blob in Fig. 1.1(b)) and yielded the outcome of timelike form factors. 
As summarized in Tables |7.13| and |7.14 for the 0~ 0~ and 0~ — > 1~ semileptonic 



decays, respectively, our numerical results for both HO and linear parameters were 
overall not only quite similar to each other but also in good agreement with the 



available experimental data To assure our method of analytic continuation, we 
used the exactly solvable model of scalar field theory interacting with gauge fields 
and analyzed the full information for the nonvalence diagram. As shown in Chapter 
5, we were able to check if the analytic continuation of the results in the g^=0 frame 
(without the black blob) indeed reproduces the exact results in timelike region. We 
established that this method can actually broaden the applicability of the standard 
LF frame a la Drell-Yan-West to the timelike form factor calculation. Exploring the 
timelike region, we have also found the importance of zero-mode contributions. The 
estimation of the zero-mode contribution was presented in Chapter 6. To the extent 



that the zero-modes have a significant contribution to some physical observables 

it seems conceivable that the condensation of zero-modes could lead to the nontrivial 

realization of chiral symmetry breaking in the LF quantization approach. 

Throughout the thesis, we attempted to fill the gap between the model wave func- 
tion and the QCD-motivated potential and explore covering as many observables as 
possible. We think that the success of our model hinges on the advantage of LF quan- 
tization realized by the rational energy-momentum dispersion relation. It is crucial to 
calculate the "good" components of the current in the reference frame which deletes 
the complication from the nonvalence Z-graph contribution. We anticipate further 
stringent tests of our model with more accurate data from future experiments and 
lattice QCD calculations. 
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Appendix A 



Spin-Orbit Wave Functions 

JJ; 

AA 



The 4x4 Dirac matrices 7^ are given by 



with 7 hermitian and 7* antihermitian. In the chiral representation flO 
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where / is the 2x2 unit matrix and a* are Pauh matrices defined as 
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(A.l) 



(A.2) 



(A.3) 



We frequently need 7"*" = 7 + 7 and 7 = -^7 7 7 7 , which are given as 



7 



The constituent quarks are spin ^ particles and can be described by Dirac 
u{k, A) and v{k, A) satisfying the Dirac equation 



(A.4) 



spmors 



m)u{k,\) = 0, {^ + m)v{k,\) = 0, 



(A.5) 
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where ^ = k/^j'^. Using the 7-matrices given by Eqs. (A.1)-(A.4), we now introduce 
the spinor basis: 



u{k, A) 



{^ + m)u{X), v{k,X) 



(^-mMA), 



u 1 I 



1 ^ \ 






















V ) 




1 1 / 



(A.6) 



(A.7) 



where v{X) — u{—X). We then can derive the spinors u{k,X) and v{k,X) as well as 
the antispinors u{k, X) and v{k, A), where we use u — u'^^^ and v — v'^^'^, as 



( m \ 


k+ 
V k^ J 



> u[k,--]^ 



k+ 


\ -m J 



v(k,~]^ 



k+ 


\ m / 

( —m \ 


k+ 
\ k^ J 



(A.8) 



k+ 

^ (0,m,-k^,k 
1-+ \ 



(A.9) 
(A.IO) 

(A.ii; 



u ( k, - \ = , — f /c^, /c^, m, 0) , u ( k, — 

V (k, 1)=^ (0, -m, -fc^, , V (k, = {k+, k\ -m, o) . (A.ll) 

Here, k^ and k^ are defined as A;^ = A;^ + i/c^ and k^ = k^ — ik"^, respectivly. The 
Melosh transformation TZM{.x,'k^,m) transforms the usual instant frame spin basis 
into the LF frame spin basis and thus enables one to assign proper total angular 
momentum quantum numbers to hadronic states. In the basis given Eq. (A.8)-(A.ll), 
the spin-orbit wave functions of pseudoscalar and vector mesons can be represented 
in the following covariant way: 



-1 



v^[M2-(mi-m2)2]V2 



u{pi,X)lbv{p2,X), 



(A.12) 
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and 



V2[M^ - (mi -ms) 2] 1/2 

£ • (Pi - P2) 



X u{pi,X) 



v(p2, A). 



(A. 13) 



Mo + mi + m2. 

The four- vectors pi,p2 are given in terms of the LF relative momentum variables 
{x, k±) as follows 



Pi 



= xP+, p+ = (1 - x)P+, 



Pi_L = xP_L + k_L, P2± = (1 -a;)P± -k_L, 



(A.14) 



and satisfied hj p'^ = mj{i = 1,2). The polarization vectors are given by 



e±{± 



^ ' Mo r ' p+ ' 

The explicit form of Eqs. (A. 12) and (A. 13) are given by 
^00 ^0 / -k^ xm2 + (: 



V2\ -xm2 - (1 - x)mi 



xm2 + (1 — a;)mi 



^11 ^ I xm2 + (I - x)mi + ^r-^^-— Jf'^T 

7? 11 — ^ ' \ / i ' Mo+mi+m2 Mo+mi+m2 



^^-R (l-a:)Mo+m2 
M0+TO1+TO2 



(fc^) 
Mo+mi+m2 



(A.15) 

(A.16) 
(A. 17) 



10 _ "^0 



7^^"- - 



J L (l-2x)Mo+m2-mi , x , 2^^ \ 

_ n _ ^^m, -U 2fci j^R {l-2x)Mo+m2-mi j ' 



'^^ ^/2 \ Tm -K'1 - -^ 2fci J R {l-2x)Mo+m2- 



Mo+Tni+m2 Mo+mi+m2 

-k'j^^S^ ^^^2 + (l-a;)mi + ,,-L 



j„L {i-x)Mo+m2 



where 



and 



7^o = 



'° ^a;(l -a;)[Mo2 - (mi - ma)^] V2 



k\ + m\ ^ k\ + ml 



X 
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1-x 



(A.18) 
(A.19) 

(A.20) 
(A.21) 
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Appendix B 

Spin- Averaged Meson Masses of 
7/, 77^, uo and (j) 



Wc list our flavor wave functions of the neutral pseudoscalar (77, rj') and vector (a;, 0) 
nonet states to show explicitly how we obtained the values of spin-averaged mass used 
in Chapter 2: 

r}^ ^{uu^ dd-2ss), (B.l) 
V6 

77' = ^(-u-u + fiJ+ss), (B.2) 
v3 

uj = ^={uu -\- dd) , (B.3) 
v2 

= ss. (B.4) 
For ideally mixed isocalar and isovector mesons, we take the flavor wave function as 

= 71 (^^ + xf ^"^^ = s-s, (B.5) 

where PiV) denotes pseudoscalar(vector) meson states. In terms of this basis, the 
neutral meson nonet states are given by 
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where t] = rjs, rj' = rji, uj = Xns a-iid (j) = Xs ■ We define the spin- averaged masses M 
of Xns and Xs as follows: 

M„, = ^M„^, + ^M^,, M, = ^Mf + ^Mj. (B.7) 

Using Eqs. (B.3)-(B.6), we assign 

Then, we obtain the following spin averaged masses of rj and rj': 

M„ = ^Mn. + ^M. = 842 MeV, (B.9) 

2 - 1 - 

M^, = -M„, + -M, = 885 MeV. (B.IO) 

3 3 

Likewise, from Eq. (B.6), we obtain 

Mu., = + ^M0, (B.ll) 
= + ^M^. (B.12) 

Using = M^^g and M,,/ = M^^, we can then evaluate the spin-averaged masses of 
u and as M,^ = 928 MeV and = 799 MeV, respectively. 

To calculate the spin-averaged meson masses depending on the schemes of fiavor 
mixing, let's consider the "perfect mixing" {9 su (3) = —10°) fj and rj' states defined 



hyMM 



V = ^{xL-X^), V' = ^{XL + Xf). (B.13) 



Using Eqs. (B.13) and Eq. (B.6), we obtain 
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Thus, the spin-averaged masses fh of the fj and r]' for the "perfect mixing" scheme 
are related to and M^i calculated in our scheme (see Eqs. (B.9) and (B.IO)), 
respectively, as follows: 

- (a/2 -1)2 (v^+l)2 
M,, = ^ + ^ 

= ^ ' + ^ > m^,. (B.15) 

Prom Eq. (B.15), the spin-averaged masses for the "perfect mixing" states are given 
by = 843 MeV and m;^, = 884 MeV, respectively. Using the same method as 
above, for the 9su{2,) — —23° mixing scheme, we obtain the following spin averaged 
meson masses: = 838 MeV and m^/ = 873 MeV, respectively. 
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Appendix C 

The Electromagnetic Decay Width 

r{A^B + 7) 

In this appendix, wc derive the decay widths oi A ^ B + 'j [A = p, K* ,LV,(j), and 
B — TTjK, T], T]') and Ai ^ tt + 7. 

Consider the electromagnetic decay process A — > S + 7. If the masses of A and B 
mesons are given by and me, respectively, then the decay rate V{A ^ B -\- ^) m. 
the rest frame of A is given by = c = 1) 

X 5{mA-EB-E^), (C.l) 

where S = for each group of j identical particles in the final states. Here, A4 is 
the transition matrix element defined by 

M^el{X,){B{P')\J'^\A{P,X)). (C.2) 

To allow decays to all possible spin configurations, we consider the replacement 
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where Sa is the spin of particle A and Xj corresponds to the state of transverse 
polarization of the emitted photon. The above replacement means the average over 
the initial spin and the sum over the helicities of the emitted photon. Then, from Eq. 
(C.l), we get 



X 



S{mA - Jml + pl - \pb\ 



E \M\'. (C.4) 



If we let \pb\ — P and E = m% + + p, then 



Therefore, we get 



A-y=±l 



where 



5 = 1, if B^7, 

S^\, if B = 7. (C.7) 

Here, po is the value of p when E — ttia, i.e., po = {'itT'A ~ fny)/'^fnA- In the rest 
frame of A, i.e.. Pa = ijriA,^) ,Pb = (^^s^Ps) and = (|ps|,— Pb), the invariant 
amplitude square is given by 



|2 



E lA^r = E \eGAB{Q'y"'^e;iX,)e^{P,X)P:,Pfsi 

= e'\GABiQ')\'2\e''-^PLPp\' = e'\GAB{Q')\'2mlpl 



{ml - m|)2 



2 

Therefore, we get the following decay width T{A — > 5 + 7) 



e'\GAB{Q')\'. (C.8) 
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where a{— e^/Air) is the fine-structure constant. 

The decay width of Ai — > ttj can be calculated in the same manner using our defini- 
tion of the transition matrix element given by Eq. (2.26) as follows: 



E \M\'=\e;{\,){n{P')\J^\A,{P,\))\' = 2 



= 2 



e 



iPA,+P',)-P, 

(2mAiP^)^ 



(C.IO) 



and the result is given by 



r(Ai ^ 7r7) 



4a 

T 



rriAi 



m 



2 \ 3 



2mAi 



(C.ll) 
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Appendix D 

Fixation of the Model Parameters 
Using the Variational Principle 

In this appendix, we discuss how to fix the parameters of our model, i.e., quark 
masses (m„, rris), Gaussian parameters ( Pun = Pud^ Pus, Pss), and potential parameters 
{a,b,K) in Vgg given by Eq. (^). In our potential model, the p-vr mass splitting is 
obtained by the hyperfine interaction Vhyp. 

Our variational method first evaluates (\I'|[ifo + with a trial function 0io(A;^) 

that depends on the parameters (m, (3) and varies these parameters until the expecta- 
tion value of Hq + Vq is a minimum. Once these model parameters are fixed, then, the 
mass eigenvalue of each meson is obtained by Mgg = (\Ef | [ifo+^o] |^I/) + (^E'|-ffhyp|^) Q- In 
this approach, we do not discriminate the Gaussian parameter set (3 = {Puu, Pus, Pss) 
by the spin structure of mesons. 

Let us now illustrate our detailed procedures of finding the optimized values of the 

-'^As we will see later, in our fitting of the p-n splitting, the rather big mass shift due to the 
hyperfine interaction is attributed to the large QCD coupling constant, k = 0.3 — 0.6. 
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model parameters using the variational principle: 

dmHo + Vo]\^) 



0. (D.i: 



d/3 

From Eqs. ([j:.lD-(^?^ and Eq. ( p.lD , we obtain the following equations for the HO 
and linear potentials: 

Lmear : k = ^[^-^ ^j, (D.3) 

where the subscript h (/) represents the HO (linear) potential parameters. Equa- 
tions (D.2) and (D.3) imply that the variational principle reduces a degree of freedom 
in the parameter space. Thus, we have now four parameters, i.e., {m„,/3„j, a, k (or 
b)}. However, in order to determine these four parameters from the two experimen- 
tal values of p and vr masses, we need to choose two input parameters. These two 
parameters should be carefully chosen. Otherwise, even though the other two pa- 
rameters are fixed by fitting the p and vr masses, our predictions would be poor for 
other observables such as the ones in Section 4.2 as well as other mass spectra. From 
our trial and error type of analyses, we find that mu=0.25 (0.22) GeV is the best 
input quark mass parameter for the HO (linear) potential among the widely used 
u{d) quark mass, m„=0.22 GeV H, 0.25 GeV 0, and 0.33 GeV ||, |80|, 111]. For 
the linear potential, the string tension bi = 0.18 GeV^ is well known from other quark 



model analyses [^8|, |80|, ^ commensurate with Regge phenomenology. Thus, we take 
m„=0.22 GeV and 6^=0.18 GeV^ as our input parameters for the linear potential case. 
However, for the HO potential, there is no well-known quantity corresponding to the 
string tension and thus we use the parameters of mu(d)=0.25 GeV and /5„j=0.3194 
GeV as our input parameters which turn out to be good values to describe various 



observables of both the vr and p mesons for the Gaussian radial wave function |10 
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Using Eqs. (D.2) and (D.3) with the input value sets of (1) (m„=0.25 GeV, 

/?„rf-=0.3194 GeV) for the HO potential and (2) (m„=0.22 GeV, 6;=0.18 GeV^) for the 
linear potential, we obtain the following parameters from the p and tt masses, viz., 
{^\H^^^\^!) = M^j^^ = nipi^^-^ {P= pseudoscalar and V= vector): 

HO : ah = -0.144 GeV, = 0.010 GeV^ kh = 0.607, (D.4) 
Linear : a/ = -0.724 GeV, = 0.3659 GeV, ki = 0.313. (D.5) 

As shown in Fig. 4.1(a), it is interesting to note that our two central potentials. 
Coulomb plus HO (solid line) and Coulomb plus linear (dotted line) potentials, are 
not much different from each other and furthermore quite comparable to the Coulomb 
plus linear quark potential model suggested by Scora and Isgur (ISGW2) |^ (long- 
dashed line for /t=0.3 and dot-dashed line for k,=O.Q) up to the range of r < 2 
fm. Those four potentials (HO, linear, and ISGW2) are also compared with the GI 
potential model |^8| (short-dashed line) in Fig. 4.1(a). The corresponding string 
tensions, i.e., fo{r) = —dVo{r)/dr, are also shown in Fig. 4.1(b). 

Next, among various sets of {ms,Pus} satisfying Eqs. (D.2) and (D.3), we find 
m,=0.48 [0.45] GeV and /3„,-=0.3419 [0.3886] GeV for HO [linear] potential by fitting 
optimally the masses of K* and K, i.e., M^j^^ = TnK*{K)- Once the set of {ms,/?„s} 
is fixed, then the parameters /3s5=0.3681 [0.4128] GeV for the HO [linear] potential 
can be obtained from Eq. (D.2) [(D.3)]. Subsequently, M]^ and are predicted 
as 996 [952] MeV and 732 [734] MeV for the HO [linear] potential, respectively. As 
shown in Fig. 4.2(a) [4.2(b)], the solid, dotted, and dot-dashed lines are fixed by 
the HO[linear] potential parameter sets of {rriu, (^ud\^ {"^s? Pus}^ and [3ss, respectively, 
and these three lines cross the same point in the space of h and if the parameters 
in Table 4.1 are used. 

We have also examined the sensitivity of our variational parameters and the corre- 
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sponding mass spectra using a Gaussian smearing function to weaken the singularity 
of 6^{r) in hyperfine interaction, viz., 

^'(0--4^exp(-aV). (D.6) 



By adopting the well-known cutoff value of a=1.8 and repeating the same 

optimization procedure as the contact term [i.e., 5^(r)] case, we obtain the following 
parameters^] for each potential: 

HO : ah = -0.123 GeV, bh = 9.89 x 10"^ GeV^ = 0.636, (D.7) 
Linear : a/ = -0.7 GeV, k = 0.176 GeV^ ki = 0.332. (D.8) 

The changes of other model parameters and mass spectra are given in Tables 4.1 and 
4.2. As one can see in Eqs. (D.7), (D.8) and Tables 4.1, 4.2, the effects of smearing out 
6^{r) are quite small and the smearing effects are in fact negligible for our numerical 
analysis in Section 4.2. 



^For the sensitivity check of smearing out S'^ir) [Eq. ( p.6| )], we kept = 0.3659 GeV for the 
Hnear potential case given by Eq. (D.5) as an input value and checked how much bi changed. 
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Appendix E 

Mixing Angles of (77,77^) and (uj^cj)) 



In this appendix, we illustrate the mixing angles of (77,77') and (co',0) by adopting the 
formulation to incorporate the quark-annihilation diagrams and the effect of SU(3) 
symmetry breaking in the meson mixing angles. 



Equations ( [4 .51 ) satisfy the (mass) eigenvalue equation 

M'\f.) = Ml\f,) (2 = 1,2). (E.l) 

Taking into account SU(3) symmetry breaking, we use the following parametrization 
for TVI^ suggested by Scadron|76|: 



2_(M^^ + 2X V2XX \ 

The parameter A characterizes the strength of the quark-annihilation graph which 
couples the J=0 uu state to /=0 uu, dd, ss states with equal strength in the exact 
SU(3) limit. The parameter X, however, pertains to SU(3) symmetry breaking such 
that the quark-annihilation graph factors into its flavor parts, with A, AX, and XX"^ 
for the uu uu{dd), uu — > ss(or ss — > uu), and ss ss processes, respectively. Of 
course, X — 1 in the SU(3) exact limit. Also, in Eq. (E.2), M^^^ and describe 
the masses of the corresponding mesons in the absence of mixing. 
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Solving Eqs. (E.l), and (E.2), we obtain Eq. (gj) and 

^ 2{Ml-Ml){Ml-Ml) 

(Mi-M4)(M|^-M4)' ^ ■ ^ 

^ (Mi-Mf+AX^-2A)- ^^-^^ 

Equations (E.3) and (E.4) are identical to the two constraints Tr(A^^)= Tr(MU and 
det(A^^)= det(M|.). The sign of 5 is fixed by the signs of the A and X from Eq. (E.5). 
Also, since Eq. (E.2) is decoupled from the subspace of {uu — dd)/\/2^ the physical 
masses of and rUp are confirmed to be the masses of and M^, respectively, 
as we used in Section 4.1 to fix the parameters (a, 6, k). 

Given the fixed physical masses of M^—uIt, and M^—nip together with Mf^{i = 1, 2), 
the magnitudes of mixing angles for rj-ri' and uj-(f) now depend only on the masses of 
and M]^, respectively, from Eqs. ([4.5|). However, from Eqs. (E.3)-(E.5), one 
can see that the sign of mixing angle depends on the sign of parameter X. While 
Xp > is well supported by the Particle Data Group [|| (-23° < e^.-gf^-^ < -10°), 
the sign of Xy is not yet definite at the present stage of phenomenology. Regarding 
the sign of Xy, it is interesting to note that 6v ~ — 3.3°(= 9su{3) — 35.26°) (i.e., 
Xv < 0) is favored in Refs. 0, ^, |S|, |§, while the conventional Gell-Mann-Okubo 
mass formula for the exact SU(3) limit (X 1) predicts 6v ~ 0° in the linear mass 
scheme and 6v ~ +3.3° (i.e., Xy > 0) in the quadratic mass scheme [|^. Our 
predictions for the uj-(j) and rj-r]' mixing angles are given in Section 4.1. 

The corresponding results of the mixing parameters Xv{p) and Xy(p) in Eqs. (E.3) 
and (E.4) are obtained for the HO (linear) potential as follows: 

Xv = 0.57 (0.73)m^ GeV^ Xy = ±2.10 (±3.08), 
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\p = 13.5 (13.3)m^ GeV^ Xp = 0.84 (0.85). (E.6) 

Our values of Xy and Xp for both HO and linear potential cases are not much different 
from the predictions of Ref. [|7^]. The reason why Ay is much smaller than Xp, i.e., 



Ap ~ 23Ay(18Ai/) in our HO (linear) case and Ap ~ 18Av in Ref. may be 
attributed to the fact that in the quark-annihilation graph, the 1 annihilation 
graph involves one more gluon compared to the ^ annihilation graph. This also 
indicates the strong departure of rj-rj' from the ideal mixing. 
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Appendix F 

Derivation of the Matrix Element 
of the Weak Vector Current 



F.l D r ell- Yan- West (g+=0) Frame 

In this appendix, we show the derivation of the matrix element of the weak vector 
current < -P2|^27^<5i|-Pi > for 0^ — > 0^ decays given in Eq. ([7!^ ). 



In the LFQM, the matrix element of the weak vector current can be calculated by 
the convolution of initial and final LF wave function of a meson as follows: 

Ai,A2,A 

X ,J^t(.,k'J%Mv%illK^,(.,kJ, (F.l) 



'Pi yPi 

where the spin-orbit wave function 7l'-''^''{x,k±) for pseudoscalar meson 
obtained from Melosh transformation is given by 

1 

and 



K% = ^7tt^=— ^^fe' >^^hMP,, A), (F.2) 



Ml = + — ^. F.3 

^ 1 - a; X ^ 
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Subsituting Eq. (F.2) into Eq. (F.l) and using the quark momentum variables given 
in Eq. (|7.2| ), one can easily obtain 

k'_L)</'i(x, k_L) 



< P2\q2l''Qi\Pi > 



2{1 -x)UiJMl-{m,-m, 



X Tr 



75(^^(2 + "^2)7^(^1 + "^1)75(^5 - rrig 



(F.4) 



where we used the following completeness relations of the Dirac spinors 
u{p, \)u{p, A) = ^ + m, ^ v{p, X)v{p, X) = — m. 

^1,2 Ai 2 



(F.5) 



In the standard q~^=0 frame where the decaying hadron is at rest, the trace terms in 
Eq. (F.4) for the "+" and "_L" components of the vector current J^=q2l^qi-, respec- 
tively, are obtaind as follows 



Tr 



75(^2 + m2h^{^i + mi)75(^g- - TUg 

-4 p'^{p2 ■ Pq + m2mq) + p'^{pi ■ pq + niiniq) + p^g{-pi ■ P2 + mim2) 
4P+ 



X 



A1A2 + k^ ■ k'_ 



x{l — x) 



(kx - q±) + 



, for /i = + 

(Al + k'l) 



(F.6) 



x{l — x) 



k± + [(mi - 1712^ + gi]kj 



for fi 



(F.7) 



where Ai=miX + mq{l —x) and k'^=k_L — xq^. Our convention of the scalar product. 
Pi 'P2={PiP2 +PrP^)/2 — Pi± ■ P2± were used to derive Eqs. (F.6) and (F.7) from the 
second line of the above equation. Substituting Eqs. (F.6) and (F.7) into Eq. (F.4), 
we now obtain the matrix element of the weak vector current < P2\q2l'^Qi\Pi > for 
fi = + [see Eq. ( [7.12| )] and _L [see Eq. ( |7.14| )] in q^=0 frame, respectively. 



F.2 Valence Contributions in q^^O Frame 
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F.2.1 Electromagnetic Form Factors 



The analysis in q^^Q frame is carried out in a reference frame where the momentum 
transfer is purely longitudinal, i.e., q±=0 and q^—q^q~. In this case, we have for the 
meson rest frame Pi±=P2±=0; 



-< 1 j^i ~ q 



+ Ml M| 



-,o±). 



(F.8) 



Pi Pi - 9+ ' 

where Mi ( M2) is the mass of the initial(final) meson. The momentum transfer square 
q^ can be rewritten in terms of the fraction r—P^/Pi—'^ — q'^ /Pi as follows 

M2 



(F.9) 



where M—M1—M2. Note from the constraint < r < l(or equivalently, < q'^ < 
Pi+) that the square of the momentum transfer is always spacelike, i.e., q^ < 0. 
Constraining r as < r < 1, we obtain the solution for r from Eq. (F.9) as follows 



2M2, 



(F.IO) 



where 
are 



-q^ > 0. For given Pi and P2, the relevant quark momentum variables 



Pi± = 



Pt 



P2± = 



(1 - x)P+, p+ = xP+, 

(1 -x)Pi_L +k_L, = xPi_L - k_L, 

(1 - x')P2± + k'x, p',-^ = x'P2± - k\, 



(F.ll) 



where x{x' — x/r) is the momentum fraction carried by the spectator q in the ini- 
tial(final) state. Note that the spectator model satisfies p^=p'^ and Vqi.=v' q±- 



185 



Taking a Lorentz frame given in Eq. (F.8), the matrix element of the "good" current 
J"*" in LHS of Eq. (|7. 1|) has the form 



Ai ,Xo ,A 



X flS(-M'x)5<J^^7-%#fl;;,(..kx). (F.12) 



Using the covariant form of the spin-orbit wave function given in Eq. ( [A. 121) , Eq. (F.12) 
is transformed into 



Tr 



X 



2MioM'2oa/(1-x)(1-x') 



[F.13) 



where Mio=y^M|) - ( nii — niqY. After some manipulations, the trace term in Eq. (F.13) 
can be reduced to 



Tr 



75(^2 + "^2)7"^(^i + "^1)75(^5- - 



X' 



where 

Ai = mqX + mq{l — x), A2 = mqx' + 777.5(1 — x'). 
One can also easily show the following identity 



X 



Al + klJAl + kl. 



(F.14) 



(F.15) 



[F.16) 



Finally, using Eqs. ( [7.1D and (F.13)-(F.16), we obtain the EM form factor resulted 
from the valence contribution in q^^O frame as follows 

2 



1 + r Jo 



A1A2 + k] 



lAl + kl^Al + kl 
which is normalized to unity at Q^=Q. 



dx j (i^kj_02(a;', k_L)0i(x, k^) 



^qil ^ ? of the first term), (F.17) 
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F.2.2 ^0 Semileptonic Decays 

While the momentum transfer square for the EM transitions of a pseudoscalar 
meson in g'^T^O frame is always negative [see Eq. (F.9)], the for the P —>■ P 
semileptonic decays is positive, i.e., mf < < (Mi — M2)^. One can easily obtain 
in terms of the fraction r as follows 

g2 = (l-r)(M2-^). (F.18) 



r 



Accordingly, one gets the two solutions for r: 

M2 



^ Ml 



Vi-V2± \J (f 1 ■ f 2)^ - 1 



[F.19) 



where f 1(1^2) being the four velocity of the parent (daughter) meson and 



The +(— ) signs in Eq. (F.19) correspond to the daughter meson recoiling in the 
positive(negative) z-direction relative to the parent meson. At zero Yecoi\{q^=q^^^) 
and maximum recoil(g^=0), r± are given by 

r40) = l, r40)= . (F.21) 



VMi 

In order to obtain the form factors f±{q^) which are independent of r±, by putting|[119 
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{P2\Q2l^Qi\Pi)\r=r^ = 2P^H{r±), (F.22) 
one obtains from Eq. ( |7.4| ) 

= ^(lTr-)H(r,)-(lTr,)g(r-)_ 

r 1 — r_ 
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Using a Lorentz frame given in Eq. (F.8) and the quark momentum variables in 
Eq. (F.ll), we obtain H{r) in Eq. (F.22) as follows 

H{r) = dx ciV02(x',k^)</.i(a;,k^)^^^^:^=, (F.24) 
where Bi=miX + mg{l — x) and B'i=mix' + mq(l — x'). 
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